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INTRODUCTION 


1.1  General  Objectives 

Research  on  combustion  of  solid  rocket  propellants  is  typically  mo¬ 
tivated  by  practical  needs  of  the  solid  rocket  business,  as  perceived  by 
propellant  chemists,  propulsion  engineers,  and  by  combustion  scientists. 
Direction  and  level  of  effort  are  usually  determined  by  current  applied 
problem  areas  such  as  need  for  large  changes  in  burning  rate  or  specific 
impulse,  encounters  with  unstable  combustion  in  development  programs,  need 
for  smokeless  exhaust  plumes,  etc.  However,  there  has  always  been  a  modest 
level  of  more  fundamental  research  aimed  at  understanding  the  underlying 
details  of  the  combustion  process.  It  is  this  kind  of  research  that  pro¬ 
vides  the  physical  insight  and  factual  base  on  which  more  applied  efforts 
are  built  and  evaluated.  Typical  of  such  research,  it  is  fundamental  in 
its  objective  of  seeking  understanding,  but  applied  in  the  sense  that 
physically  and  chemically  complex  problems  are  accepted  in  order  to  pre¬ 
serve  relevance  to  the  practical  problem.  The  present  program  has  been 
designed  to  help  serve  this  need  for  "fundamental"  research.  The  scope 
has  been  limited  for  practical  reasons  to  combustion  of  heterogeneous 
systems,  without  powdered  metal  fuel  ingredients.  Objectives  are  to  de¬ 
termine  the  structure  of  the  combustion  zone,  the  behavior  of  the  indi¬ 
vidual  ingredients  in  the  combustion  zone,  and  the  relative  importance 
of  competing  processes.  As  a  means  of  controlling  variables,  combustion 
experiments  have  until  recently  used  ordered  microstructure  (e.g.,  oxidizer- 
binder  sandwiches).  From  a  practical  viewpoint,  the  objectives  are  to 
assist  in  selection  and  control  of  ingredient  characteristics,  design  of 
propellants,  prediction  of  steady  and  unsteady  combustion,  and  in  trouble¬ 
shooting  for  systems  that  are  not  meeting  specifications. 

1.2  Historical  Background 

Combustion  of  solid  rocket  propellants  is  an  applied  science  that 
has  its  origins  in  the  science  of  internal  ballistics  of  guns.  Indeed, 
the  propellant  charges  used  in  the  early  days  of  the  modern  era  of  rocketry 
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(1937  -  1944)  were  made  from  gun  propellants  and  trench  mortar  sheet  pro¬ 
pellant.  These  propellants  consisted  primarily  of  nitrocellulose,  with 
various  amounts  of  plasticizer  (often  including  nitroglycerin).  The  pro¬ 
pellants  were  colloidal  in  nature,  virtually  homogeneous.  Rocket  applica¬ 
tions  called  for  predictable,  smooth  surface  burning  of  relatively  large 
pieces  of  propellant,  at  pressures  much  lower  than  in  guns.  Because  of 
urgent  need,  development  of  such  propellant  charges  proceeded  on  a  sub¬ 
stantially  trial-and-error  basis,  while  a  science  of  internal  ballistics 
of  rockets  was  developed  as  the  needs  of  the  immediate  situation  dictated. 
Progress  of  the  1940's  is  typified  by  the  book  "Internal  Ballistics  of 
Solid-Fueled  Rockets"  by  Wimpress  (Reference  1),  and  by  a  collection  of 
papers  on  theory  of  combustion  of  solid  propellants,  in  the  Journal  of 
Physical  and  Colloid  Chemistry,  Volume  54,  of  1950  (Reference  2).  It  is  worthy 
note  that  the  science  of  propellant  combustion  remained  substantially 
empirical,  with  heavy  reliance  on  laboratory  tests  to  determine  burning 
rates  of  modified  propellants;  and  with  little  or  no  capability  to  describe 
or  anticipate  such  all-too-common,  unwanted  behavior  as  spontaneous  os¬ 
cillatory  combustion. 

By  1950,  significant  progress  had  been  made  in  development  of  pro¬ 
pellants  made  from  hydrocarbon  fuel  and  inorganic  crystaline  oxidizers. 

In  such  propellants,  the  granular  oxidizer  was  contained  in  a  matrix  of 
rubber-like  fuel.  The  theoretical  science  of  propellant  combustion  was 
hardly  prepared  to  deal  with  the  new  complexity  of  a  heterogeneous  pro¬ 
pellant,  and  the  progress  of  the  1950' s  was  largely  dependent  on  the  em¬ 
pirical  methods  and  such  physical  insight  as  could  be  developed  from  ex¬ 
perimental  research.  It  didn’t  take  too  much  physical  insight  to  recognize 
that  a  propellant  that  consisted  of  a  near  stochiometric  ratio  of  oxidizer 
and  fuel  would  burn  faster  than  the  usual  fuel-rich  mixture,  or  that  fine 
oxidizer  size  yielded  higher  burning  rate.  It  took  more  effort  to  deter¬ 
mine  the  effects  of  such  propellant  variables  on  the  sensitivity  of  burning 
rate  to  pressure,  temperature  and  flow  environment  in  the  rocket  motor. 

But  laboratory  tests  could  be  run,  and  motor  firings  could  be  made  until 
the  propellant  developers  got  what  they  wanted.  Financial  support  for 
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research  aimed  at  principles  of  combustion  languished.  In  an  American 
Rocket  Society  Meeting  in  the  late  1950's,  a  panel  of  propellant  developers 
was  asked  what  help  research  had  been  to  them,  and  they  said,  "none"  (pro¬ 
bably  not  entirely  accurate,  but  compatible  with  the  relative  amount  of 
past  effort  devoted  to  trial  and  error  development  vis-a-vis  research). 

In  1959,  the  role  of  research  was  drastically  changed  by  the  increas¬ 
ing  cost  of  motor  testing  associated  with  development  of  larger  and  more 
sophisticated  motors  and  new  propellants.  A  Defense  Department  Ad  Hoc 
Committee  (Reference  3)  studied  the  risks  to  large  motor  programs  posed 
by  the  problem  of  combustion  instability,  and  recommended  an  aggressive 
combustion  research  effort,  to  provide  a  more  fundamental  basis  for  anti¬ 
cipating  and  correcting  problems  in  development  programs.  Combustion 
instability  was  used  in  Congressional  briefing  by  liquid  propellant  pro¬ 
ponents  as  an  argument  against  use  of  solid  rockets  in  ballistic  missile 
programs  (even  though  combustion  instability  was  a  proven  problem  in  liquid 
rockets,  an  unknown  in  large  solid  rockets). 

By  this  time  a  much  wider  range  of  composite  propellants  came  into 
consideration  or  use,  further  complicating  the  problem  of  predicting  com¬ 
bustion  behavior.  The  years  1958  -  1970  were  extremely  productive  years 
from  the  standpoint  of  increased  understanding  of  propellant  combustion. 

The  role  of  research  was  particularly  evident  in  the  problem  of  combustion 
instability,  primarily  because  it  involved  a  complex  transient  interaction 
of  combustion,  gas  dynamics  of  the  combustor,  and  geometry  of  the  combus¬ 
tor.  Propellant  developers  were  not  equipped  to  deal  with  such  a  problem, 
and  motor  designers  could  not  solve  the  problems  by  trial-and-error  firings 
of  motors  because  of  the  costly  motor  firings  involved.  While  laboratory 
scale  burners  could  be  used  to  advantage  they  had  to  be  developed  by  re¬ 
search  teams.  Only  limited  aspects  of  the  problem  could  be  tested  this 
way,  and  considerable  advance  in  understanding  was  required  even  to  under¬ 
stand  the  laboratory  burner  behavior  (Reference  4).  This  situation  has 
led  to  a  much  more  conspicuous  role  for  research,  the  benefits  of  which 
were  not  limited  to  solving  combustor  stability  problems.  A  panel  of  pro¬ 
pellant  developers  might  today  grudgingly  concede  that  combustion  research 


has  provided  improved  laboratory  test  methods,  better  interpretation  of 
test  results,  better  advance  identification  of  key  problems,  and  even  help¬ 
ed  bring  desirable  new  propellant  formulations  to  successful  application. 

Recent  trends  in  propellants  and  applications  continue  to  produce 
new  needs  for  understanding.  Interest  has  grown  in  propellants  that  pro¬ 
duce  no  optically  visible  exhaust  trail  from  the  missile,  a  trend  that 
poses  combustion  stability  problems  associated  with  exclusion  of  aluminum 
from  the  formulation.  In  these  applications,  even  the  traditional  oxidizer 
(ammonium  perchlorate)  is  undesirable  because  of  its  tendency  to  give  a 
water  condensation  trail  due  to  the  HC1  in  the  reaction  products.  Efforts 
to  change  to  nitramines  (e.g.,  HMX)  have  posed  difficulties  over  the  whole 
spectrum  of  combustion  problems.  Another  general  problem  in  modern  pro¬ 
pellants  comes  from  the  ever-present  desire  to  increase  the  energy  content 
of  the  propellant  (particularly  important  in  upper  stage  and  space  motors, 
and  in  volume-limited  missile  systems).  The  ingredients  of  such  propel¬ 
lants  (e.g.,  nitrocellulose,  nitroglycerin,  HMX)  are  individually  and  col¬ 
lectively  detonable.  The  processing  industry  has  worked  the  miracle  of 
learning  how  to  produce  motors  with  these  dangerous  ingredients  (which 
are  considerably  less  delicate  once  they  are  manufactured  into  a  propel¬ 
lant),  but  they  still  represent  a  serious  hazard  if  the  combustion  does 
not  proceed  as  planned.  As  a  result  there  has  been  increased  research 
on  combustion  as  it  relates  to  transition  to  detonation. 

In  the  face  of  changing  military  and  space  program  needs,  the  problems 
emphasized  in  propellant  combustion  research  have  shifted  on  a  time  scale 
of  from  one  to  five  years,  with  the  most  fundamental  efforts  traditionally 
facing  a  struggle  for  survival  because  of  lesser  apparent  relevance  to 
solution  of  "today's"  problem.  No  matter  that  the  fundamental  knowledge 
contributes  to  all  problems,  it  is  not  identified  easily  enough  with  today's 
problem  to  compete  for  support.  Unlike  low  support  in  the  1950's,  there 
is  considerable  support  for  scientific  work  on  propellant  combustion  in 
the  1970's  (albeit  less  than  in  the  1960's),  a  result  of  increased  complex¬ 
ity  of  problems  and  of  increased  cost  of  "solving"  them  by  empirical  means. 
But  the  problem-oriented  nature  of  the  motivation  for  support  tends  to 
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obscure  Che  fact  Chat  science  and  scientists  are  not  produced  by  solving 
practical  problems,  but  rather  by  study  of  the  underlying  principles  and 
by  advancement  of  understanding  and  methods.  For  this  reason,  a  modest 
level  of  research  on  combustion  mechanisms  has  continued  for  the  last  30 
years,  as  an  essential  part  of  capability  to  deal  with  "today's"  problems. 
Such  "fundamental"  research  is  only  fundamental  in  the  context  of  an  ap¬ 
plied  science  such  as  propellant  combustion,  and  in  the  context  of  its 
motivation,  which  is  to  understand.  In  this  sense,  the  work  described 
in  this  report  is  fundamental. 

1.3  Current  Situation 

One  may  examine  the  current  situation  in  two  ways  (a)  What  kind  of 
practical  combustion  problems  do  we  have?,  or  (b)  What  aspects  of  the  com¬ 
bustion  mechanisms  do  we  understand,  what  do  we  not  understand?  These 
questions  are  addressed  here. 

Of  course,  one  would  ideally  like  to  be  able  to  predict  in  advance 
all  combustion  characteristics  of  a  given  propellant,  and  even  calculate 
modifications  in  propellant  to  produce  desired  effects  in  combustion  char¬ 
acteristics.  However,  such  a  capability  is  not  likely  to  be  achieved, 
primarily  because  of  the  complexity  of  propellants  and  of  their  combustion. 
It  has  been  suggested  that  such  goals  might  be  achieved  by  a  massive  sys¬ 
tematic  experimental  investigation  of  the  effect  of  propellant  and  envi¬ 
ronmental  variables  and  suitable  correlation  of  test  results  and  test  var¬ 
iables.  However,  the  scope  of  such  an  effort  is  almost  boundless,  and 
the  cost  prohibitive  unless  restricted  to  a  very  limited  set  of  variables. 
Thus  it  appears  to  be  preferable  to  seek  the  capability  to  predict  trends, 
and  leave  the  "fine  tuning"  of  propellant  formulation  to  the  time  when 
the  application  has  established  constraints  such  as  the  domains  of  acces¬ 
sible  propellant  variables,  expected  environmental  variables,  and  accept¬ 
able  combustion  behavior.  This  is,  in  fact,  the  strategy  that  has  been 
followed,  although  probably  more  out  of  expediency  than  logical  choice. 

The  classical  list  of  practical  combustion  problems  and  research  ob¬ 
jectives  will  probably  always  be  relevant  to  some  degree.  The  objectives 


include  assurance  of  the  following: 

1.  Desired  burning  rates,  uniformity  of  rate. 

2.  Desired  (usually  minimal)  dependence  of  burning  rate  on  variables 
of  the  combustor  environment  (pressure,  temperature,  etc.). 

3.  High  combustion  efficiency. 

A.  Stable  combustion  (non-oscillatory) . 

5.  Adequate  ignition  properties,  quenching  properties  in  some  ap¬ 
plications  . 

6.  Low  smoke  in  combustion  products  (some  applications). 

7.  Low  susceptibility  to  detonation. 

8.  Physical  and  chemical  stability  in  storage. 

9.  Minimum  toxicity  of  combustion  products. 

Unfortunately,  resolution  of  combustion  problems  is  never  absolute,  because 
the  problems  are  not  independent  of  each  other;  resolution  of  one  problem 
often  worsens  others.  This  is  the  very  reason  that  purely  empirical  fixes 
of  problems  have  such  limited  applicability.  New  applications  demand  new 
trade-offs,  and  there  is  no  substitute  for  "mechanistic  understanding" 
in  design  of  the  trade-off. 

An  illustration  of  strengths  and  weaknesses  in  our  present  knowledge, 
and  of  the  impact  of  demands  by  new  applications,  is  provided  by  the  trend 
to  low  smoke  motors  in  recent  years.  This  trend  has  had  two  major  impacts 
on  propellants,  with  companion  combustion  problems.  Powdered  metal  ingre¬ 
dients  are  eliminated,  and  efforts  are  being  made  to  avoid  use  of  ammonium 
perchlorate.  The  metals  yield  condensed  oxide  smoke  in  the  exhaust;  and 
AP  yields  HC1 ,  which  causes  water  condensation  in  the  exhaust  trail.  But 
the  oxide  smoke  from  metals  had  virtually  eliminated  oscillatory  combustion 
in  earlier  motors,  and  now  the  problem  is  back  "full  force".  The  proposed 
abandonment  of  AP  is  causing  consideration  of  a  new  spectrum  of  ingredients 
such  as  nitramines,  which  have  very  different  combustion  characteristics 
that  are  only  poorly  understood.  These  are  driving  forces  in  a  large  part 
of  today's  propellant  combustion  research,  some  of  which  was  abandoned 
years  before  because  of  apparent  decline  of  need  at  that  time. 
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Another  need  that  has  heavily  affected  recent  research  is  the  hazard 
of  very  high  energy  propellants  selected  for  those  applications  where  vi¬ 
ability  of  the  weapon  system  is  dependent  on  attaining  high  specific  im¬ 
pulse.  This  has  necessitated  better  understanding  of  detonation  hazards, 
a  situation  that  has  led  to  research  on  both  mechanical  and  combustion 
behavior  of  heterogeneous  solids . 

Both  the  oscillatory  combustion  and  the  detonation  hazards  problems 
are  made  relatively  intractable  by  the  necessity  to  embrace  the  complex¬ 
ities  of  transient  nonlinear  processes,  in  heterogeneous  systems. 

In  science,  problems  involving  complex  systems  are  necessarily  ap¬ 
proached  by  a  great  deal  of  observation,  ordering  of  observations,  recog¬ 
nition  of  critical  features  of  behavior,  construction  of  hypotheses  and 
analytical  models  to  explain  and  exhibit  the  observed  features  of  behavior, 
followed  by  cycles  of  refinement  through  experiments  and  analysis  (a  some¬ 
what  grand  example  is  the  entire  science  of  cosmology).  Success,  however 
viewed,  quite  often  depends  on  understanding  of  processes  originating  at 
the  microscopic  level  (in  the  case  of  cosmology,  understanding  at  the  atomic 
and  nuclear  level),  and  this  is  true  of  propellant  combustion,  where  re¬ 
action  zones  and  thermal  layers  that  control  the  macroscopic  burning  are 
of  the  same  dimensional  order  as  the  granularity  of  the  propellant  ingre¬ 
dients.  Unfortunately,  observational  methods  generally  do  not  have  ade¬ 
quate  spatial  resolution  at  those  dimensions.  In  spite  of  this  difficulty, 
it  is  necessary  to  obtain  knowledge  regarding  the  combustion  zone  structure 
in  order  to  recognize  critical  features  of  behavior  and  construct  relevant 
models . 

The  present  state  of  theory  of  propellant  combustion  is  somewhat  like 
that  of  cosmology  in  1938.  Many  observations  of  the  macroscopic  aspects 
of  combustion  and  a  bountiful  and  growing  array  of  hypotheses  and  analy¬ 
tical  models  are  in  hand,  but  there  is  only  minimal  information  on  the 
microscopic  processes  that  drive  the  overall  phenomenon  (nuclear  processes 
in  cosmology,  molecular  and  microflame  processes  in  propellant  combustion). 
This  difficulty  has  not  been  too  serious  in  developing  an  understanding 
of  steady  state  combustion,  because  only  a  time-wise  average  of  the  micro- 
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processes  is  important.  In  transient  processes,  the  time  constants  of 
the  micro  processes  (mixing  times,  flame  transit  times,  collision  times) 
become  important  aspects  of  the  macroscopic  behavior,  and  the  knowledge 
to  deal  with  this  fact  is  largely  lacking.  A  classic  example  is  the  in¬ 
ability  to  explain  the  dynamic  response  of  combustion  to  pressure  waves 
in  high  frequency  oscillatory  burning  (i.e.,  explain  the  effect  of  mean 
pressure,  oxidizer  particle  size,  and  ballistic  modifiers  on  dynamic  com¬ 
bustion  response).  Another  example  is  the  response  of  a  burning  porous 
bed  to  the  penetration  of  hot  gases  ahead  of  the  reaction  wave  (an  element 
of  the  def lagration-to-detonation  transition  process).  There  is  no  way 
these  processes  can  be  correctly  explained  without  looking  at  the  situation 


on  a  three-dimensional  microscopic  scale  and  describing  the  transient  pro¬ 


cesses  on  that  basis. 

One  might  summarize  the  present  situation  in  combustion  science  of 
propellants  by  the  following: 

1 •  The  response  of  propellant  ingredients  to  controlled  heating  has  been 
studied  by  a  variety  of  processes  in  the  hopes  that  response  in  the  pro¬ 
pellant  combustion  zone  can  be  learned  under  less  difficult  conditions. 
Tests  have  never  been  standardized,  results  have  not  been  catalogued,  and 
the  various  studies  have  not  used  the  same  ranges  or  increments  of  control 
variables.  Many  ingredients  have  been  tested  only  superficially.  For 
some  ingredients  (e.g.,  ammonium  perchlorate)  the  test  temperatures  have 
remained  too  low  to  simulate  combustion  zone  behavior.  Very  little  of 
this  work  is  being  done  in  the  U.S.A.  anymore.  The  results  remain  frag¬ 
mentary. 

2.  Deflagration  and  combustion  of  ingredients  under  controlled  conditions 
is  often  done  as  a  means  of  overcoming  the  limitations  of  relevance  of 
thermal  response  tests,  while  avoiding  the  difficulties  of  observation 
with  propellants.  These  studies  are  often  a  source  of  embarassment  to 
modelers  because  they  reveal  complex  behavior  not  included  in  analytical 
models.  They  eventually  lead  to  better  models,  and  more  realistic  raech- 
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anistic  insight  (deflagration  of  ammonium  perchlorate  crystals  is  the  most 
conspicuous  example  of  this  kind  of  study).  At  present  almost  no  work 
of  this  kind  is  going  on  in  the  U.S.A.  except  with  HMX. 

3.  The  steady  state,  one  dimensional  representation  of  propellant  burning 
can  be  coded  for  computers  and  solutions  in  principle  can  be  obtained  if 
suitable  reaction  kinetics  can  be  established.  However  such  models  have 
not  to  date  provided  for  the  phenomenon  of  selective  concentration  of  stable 
ingredients  or  reaction  intermediates  on  the  burning  surface,  a  process 
usually  observed  in  experiments.  Further,  the  reaction  kinetics  are  still 
undetermined  for  most  ingredients.  Although  propellant  combustion  calcula¬ 
tions  can  sometimes  be  technically  justified  from  one-dimensional  models, 
such  calculations  are  not  used  much  because  they  are  not  relevant  to  hetero¬ 
geneous  propellants.  Of  course  burning  rates  can  be  measured  by  a  variety 
of  methods. 

4.  The  present  modeling  of  transient  combustion  (one-dimensional)  is 
essentially  the  same  as  it  was  ten  years  ago.  Improvements  are  probably 
needed  in  modeling  the  high  frequency  response,  as  problems  exist  with 
combustion  instability  in  low-smoke  tactical  rockets.  However  improvements 
have  been  slow  becuse  of  lack  of  kinetic  data  on  surface  and  gas  phase 
reactions,  inapplicability  of  1-D  models  to  heterogeneous  propellants, 

and  mistrust  of  experimental  data  available  to  evaluate  models.  On  the 
other  hand,  the  experimental  results  from  T- Burner  tests  are  widely  used 
as  a  basis  for  combustor  stability  calculations . 

5.  Steady  state  burning  of  heterogeneous  propellants  has  been  a  continu¬ 
ing  topic  of  study,  both  experimental  and  analytical  approaches  being  used. 
Most  of  the  experimental  work  has  been  aimed  at  determining  burning  rates , 
and  the  effects  of  propellant  and  environmental  variables  on  burning  rates. 
Analytical  models  are  based  on  the  calculation  of  heat  transfer  rates  from 
reaction  zones  back  to  the  propellant  surface  (a  one-dimensional  represen¬ 
tation).  The  propellant  heterogeneity  enters  in  through  the  calculation 


of  the  effective  distances  of  various  exothermic  reaction  steps  from  the 
surface,  a  calculation  based  on  assumptions  and  calculations  of  the  micro¬ 
structure  of  the  combustion  zone.  Analytical  modeling  work  has  been  at 
a  relatively  high  level  for  several  years.  Recent  advances  in  theory  have 
mostly  been  concerned  with  "liberalizing"  the  assumptions  about  microstruc¬ 
ture,  primarily  by  inclusion  of  more  versatile  descriptions  of  particle 
size  distribution.  However  these  changes  do  not  include  much  improvement 
in  description  of  the  state  or  configuration  of  the  propellant  burning 
surface  and  the  real  three-dimensionally  complex  interaction  of  gas  phase 
flames  and  heat  transfer  with  the  surface.  Further,  experimental  studies 
of  these  processes  are  also  minimal.  The  vigorous  modeling  activities, 
lacking  new  experimental  inputs,  tend  to  become  preoccupied  with  details 
and  speculations. 

6.  Erosive  burning  is  a  phenomenon  that  has  resisted  decisive  measurement 
and  modeling,  probably  partly  because  engineering  schemes  for  prediction 
have  been  good  enough  to  invalidate  the  need  for  a  more  decisive  attack 
on  the  problem.  However  the  low  level  of  progress  in  this  area  is  also 
attributable  to  lack  of  sufficient  knowledge  of  combustion  zone  microstruc¬ 
ture,  difficulty  in  laboratory  scale  simulation  of  the  motor  flow  environ¬ 
ment,  difficulty  in  measuring  burning  rate  under  such  conditions,  and  in¬ 
tractability  of  relevant  analytical  models  of  the  combined  combustion-flow 
process.  Recently  increased  efforts  have  made  significant  progress  in 
modeling  the  effect  of  flow  field  on  combustion  zone  and  microstructure. 

The  models  are  necessarily  less  rigorous  representations  of  combustion 
zone  microstructures  than  the  no-flow  models  because  of  the  necessity  to 
incoporate  the  complication  of  the  flow  field.  Such  experimental  work 
as  has  been  conducted  is  aimed  at  determining  erosive  burning  rate  itself, 
rather  than  at  obtaining  the  information  on  combustion  zone  structure  needed 
to  improve  the  model.  As  always,  experiments  in  laboratory  scale  apparatus 
raise  doubts  about  adequacy  of  simulation  of  the  rocket  motor  flow  field. 
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7.  Oscillatory  combustion  of  heterogeneous  propellants  is  a  matter  of 
great  practical  concern,  and  a  logical  extension  of  topic  #4,  with  the 
added  complication  of  a  three-dimensionally  complex  combustion  zone.  It 
is  helpful  to  discuss  the  topic  in  2  parts,  because  the  status  differs 
markedly  relative  to  them.  The  first  part  is  called  pressure-coupling, 
i.e.,  the  dynamic  response  of  the  combustion  to  pressure  oscillations. 

The  other  part  is  called  velocity  coupling,  i.e.,  the  dynamic  response 
of  the  combustion  to  gases  oscillating  parallel  to  the  burning  surface. 

As  one  might  anticipate  by  analogy  with  steady  state  burning,  the  state 
of  knowledge  of  pressure  coupling  (analogous  to  pressure  dependence  of 
steady  state  rate)  is  much  more  advanced  than  the  state  of  knowledge  of 
velocity  coupling  (analogous  to  "erosive"  burning).  If  for  no  other  rea¬ 
son,  this  is  due  to  the  macroscopically  one  dimensional  nature  of  the  pro¬ 
cess  in  pressure  coupling,  as  compared  to  the  necessarily  two  dimensional 
nature  of  the  process  in  velocity  coupling. 

Experimental  methods  and  analytical  models  of  pressure  coupled  oscil¬ 
latory  combustion  have  evolved  continuously  since  1956,  with  modern  de¬ 
velopments  mainly  extending  the  range  of  frequency  and  types  of  propellants 
whose  behavior  can  be  described.  Measurements  are  costly  and  only  semi- 
quantitative,  being  better  for  nonaluminized  propellants  than  for  alumi¬ 
nized  ones.  Analytical  models  have  until  recently  been  purely  one-dimen¬ 
sional.  The  most  recent  efforts  using  statistical  description  of  the 
combustion  zone  rely  on  one-dimensionalizations  and  steady  state  calcula¬ 
tions  at  various  points  in  the  computational  sequence,  while  failing  to 
describe  the  transient  three-dimensional  behavior  at  the  microscopic  level. 
This  latter  aspect  of  the  problem  has  thus  far  escaped  experimental  study 
as  well.  In  general,  experimental  results  at  all  levels  have  been  from 
programs  of  limited  scope  in  many  laboratories  without  benefit  of  any  stan¬ 
dardized  increments  of  propellant  or  test  variables,  so  few  generalizations 
can  be  made  from  experimental  results  (a  JANNAF  committee  to  establish 
standards  was  unable  to  gain  acceptance).  However  results  have  been  useful 
in  selection  of  propellants  for  development  motors  that  were  expected  or 
known  to  have  combustor  stability  problems,  and  the  test  methods  and  services 
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were  the  product  of  combustion  research  teams. 

The  state  of  knowledge  of  velocity  coupled  combustion  response  is 
very  poor  indeed,  although  a  matter  of  particular  concern  due  to  experiences 
in  development  programs.  Valid  analytical  modeling  is  even  more  difficult 
than  in  the  case  of  erosive  burning,  and  the  results  of  a  valid  model  are 
in  a  form  that  is  difficult  to  apply  to  motors  (in  effect,  the  model  calcu¬ 
lation  has  to  be  repeated  at  time  intervals  during  burning  in  every  motor, 
in  combination  with  the  charge  geometry  and  internal  ballistic  equations 
of  the  motor)  (see  Appendix  A).  Models  to  date  are  few,  complex,  unevaluated, 
but  clearly  naive  as  far  as  relevance  to  behavior  of  ingredients  of  hetero¬ 
geneous  propellants  and  combustion  zone  microstructure  are  concerned. 

Further,  velocity-coupled  combustor  instability  has  certain  nonlinear  char¬ 
acteristics  that  have  proven  to  be  very  much  a  part  of  the  difficulty  in 
development  programs,  but  singularly  intractable  in  analytical  modeling 
and  in  design  of  laboratory-scale  simulation  of  motors.  As  in  the  case 
of  erosive  burning,  the  problem  is  to  produce  a  flow  field  adjacent  to 
the  burning  surface  that  simulates  that  in  the  rocket  motor,  except  that 
in  the  velocity  coupled  instability  experiment  one  must  also  simulate  cor¬ 
rectly  the  gas  oscillations ,  and  measure  the  effect  of  the  combustion  response 
on  the  oscillations.  The  so-called  velocity-coupled  T-Burner  apparently 
does  exhibit  velocity  coupled  response,  and  can  probably  distinguish  between 
propellants  with  high  and  low  velocity-coupled  response.  However  simulations 
of  motor  flow  fields  is  very  superficial,  interpretation  of  test  results 
is  speculative,  results  are  not  very  reproducible,  and  the  method  is  current¬ 
ly  usable  only  over  a  limited  range  of  test  variables.  No  serious  effort 
has  been  made  to  explore  experimentally  the  primary  source  of  nonlinear 
behavior,  which  relates  to  the  superposition  of  mean  flow  field  and  parallel- 
acoustic  oscillations.  An  independent  approach  to  measurement  of  velocity- 
coupled  response,  using  a  system  with  driven  oscillations,  is  under  develop¬ 
ment  at  one  laboratory,  and  has  some  potential  for  better  measurements 
in  the  frequency  range  where  the  T-Burner  is  least  satisfactory  (low  frequen¬ 
cy)  . 
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8.  Transition  from  Def iagration-to-Detonation  is  the  combustion  topic 
under  most  intensive  study,  and  the  most  complex  one.  The  phenomenon  in¬ 
volves  dynamic  mechanical  failure  of  the  propellant  as  much  as  it  involves 
combustion.  In  addition,  it  involves  extremely  high  pressures  and  short 
times,  so  that  it  is  not  clear  what  chemical  steps  actually  have  time  to 
participate  in  transition.  Analytical  modelers  are  confronted  with  analysis 
that  bridges  the  domain  from  classical  diffusion-controlled  combustion 
to  shock  propagated  detonation;  an  analysis  that  is  almost  prohibitively 
difficult  for  a  pure  detonable  liquid  or  pre-mixed  gas.  Given  the  added 
complications  of  a  heterogeneous,  chemically  complex  solid  with  possible 
porosity,  it  seems  unlikely  that  any  comprehensive  solution  of  the  problem 
will  be  achieved.  Under  these  conditions  the  analytical  modelers  can  run 
wild  solving  dozens  of  simplified  models,  while  critics  can  show  that  the 
models  are  patently  naive.  Ultimately,  an  integrated  program  of  analysis 
and  experiments  must  emerge,  that  establishes  and  models  the  key  aspects 
of  the  phenonemon  and  identifies  what  properties  of  practical  propellants 
prevent  an  accelerating  reaction  front.  Progress  in  research  is  limited 
not  only  by  the  difficulty  of  identifying  the  key  processes  to  study  and 
model,  but  also  by  the  difficulty  of  designing  and  conducting  corresponding 
experiments  (complicated  further  by  hazardous  materials  and  the  risk  of, 
or  need  for  detonations). 

It  is  not  intended  to  summarize  here  the  status  of  research  and  know¬ 
ledge  on  the  topic  of  DDT,  which  is  still  somewhat  primitive.  However 
some  points  merit  particular  note: 

a)  To  the  extent  that  mechanical  degradation  of  the  solid  propellant 
under  dynamic  loading  is  involved,  the  theory  is  in  very  bad  shape  for 
heterogeneous  materials  (put  more  charitably,  a  relevant  theory  is  incred¬ 
ibly  complex). 

b)  The  high  temperatures  and  rapid  changes  involved  during  the  accel¬ 
erating  and  "shock  building"  phase  of  transition  probably  invalidate  the 
use  of  conventional  combustion  concepts  such  as  "ignition"  and  "burning 
rate",  which  are  still  used  by  modelers. 

c)  It  is  not  clear  what  chemical  reactions  actually  have  time  to 
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contribute  to  the  acceleration  of  reaction  (a  more  fundamental  statement 
of  (b)). 

d)  It  is  not  yet  clear  whether  transition  can  occur  in  a  practical 
propellant  without  prior  mechanical  degradation,  leaving  unclear  the  issue 
of  the  kind  of  processes  to  include  in  design  of  experiments  and  analyses. 

1.4  Approach  in  Present  Investigations 

In  the  present  program,  propellant  combustion  has  been  studied  by 
examination  of: 

1.  Physical,  chemical  and  combustion  behavior  of  propellant  ingre¬ 
dients  at  elevated  temperatures . 

2.  Similar  behavioral  studies  of  ingredient  combinations. 

3.  Suitable  means  to  describe  the  microstructure  of  a  heterogeneous 
propellant  and  its  burning  surface. 

4.  Structure  of  the  combustion  zone  of  heterogeneous  systems. 

5.  Analytical  modeling  of  ideal  systems. 

6.  Combination  of  the  results  in  models  most  appropriate  for  their 
application  (e.g.,  representation  of  perturbation  response  to 
pressure,  velocity;  erosive  burning;  role  of  ballistic  modifiers; 

a  context  for  description  of  surface  accumulation  as  with  aluminum). 
There  is  nothing  particularly  novel  about  this  strategy,  but  in  many  critical 
respects  past  work  has  been  indecisive,  leaving  the  objectives  of  understand¬ 
ing  and  valid  description  of  propellant  combustion  in  a  speculative  state. 

The  investigations  reported  here  were  designed  to  answer  outstanding  questions 
of  understanding  or  of  validity  of  modeling.  While  gravely  constrained 
by  the  difficulties  of  designing  relevant  and  tractable  experiments  and 
analytical  models,  some  progress  in  development  of  experimental  methods 
is  reported,  and  results  of  experiments  and  analysis  show  some  microscopic 
and  transient  aspects  of  the  combustion  process  that  are  important  in  pro¬ 
pellant  combustion  but  neglected  in  past  work. 
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THERMAL  DEGRADATION  OF  PROPELLANT  BINDERS 


In  most  models  of  the  combustion  of  heterogeneous  propellants,  it 
is  assumed  that  the  binder  is  degraded  from  a  solid  to  a  gas  in  place, 
with  the  geometry  of  its  exposed  surface  regressing  in  accordance  with 
the  local  heat  transfer  to  produce  gasification  temperatures  (Figure  la). 

The  heat  to  the  complex  array  of  binder  is  supplied  by  a  coupled,  complex 
combination  of  exothermic  reactions  in  the  oxidizer  and/or  gas  phase. 

In  analytical  models,  the  geometrical  details  of  this  coupled  process  are 
ignored,  and  the  binder  surface  is  assumed  to  be  flat,  with  a  regression 
rate  of  the  propellant  (Figure  lb)  (without  any  argument  as  to  why  the 
surface  would  be  flat).  In  practice,  burned  surfaces  of  propellants  (quench¬ 
ed  by  various  means)  exhibit  complicated  surface  configurations,  with  the 
details  depending  on  the  types  of  binder  and  oxidizer,  combustion  pressure, 
nature  of  gas  flow  adjoining  the  burning  surface  (e.g..  Figure  lc).  These 
details  of  the  combustion  process  have  not  been  determined  systematically, 
and  their  importance  is  unknown  (omission  from  combustion  models  is  more 
a  matter  of  simplification  of  mathematics  and  ignorance  of  details  than 
physical  justification). 

In  order  to  better  anticipate  or  explain  the  role  of  binder  in  the 
combustion  zone  process,  it  is  natural  to  explore  the  response  of  the  pure 
binder  material  to  heating  in  the  temperature  range  characteristic  of  the 
propellant  burning  surface  (up  to  several  hundred  degrees  C).  Such  exper¬ 
iments  have  been  done  in  the  past,  using  combustion  of  samples,  pyrolysis 
by  radiant  energy,  and  electrical  heating  of  samples.  Observations  includ¬ 
ed  photography,  thermal-gravemetric  measurements,  differential  scanning 
calorimetry,  etc.  Prior  experience  in  the  present  program  had  indicated 
that  "melting"  and  flow  of  the  binder  might  play  an  important  role  in  pro¬ 
pellant  combustion  (the  term  is  used  to  refer  to  a  state  of  the  binder 
and/or  intermediate  decomposition  products  that  will  flow).  Accordingly 
it  was  decided  to  use  a  hot  stage  microscope  to  observe  the  physical  state 
and  activity  of  binder  samples  during  heat-up  to  1200°  C.  Tests  were  done 
at  atmospheric  pressure  in  air,  oxygen  and  argon  atmospheres  with  vapor 


a.  b. 


(c) 


Figure  1.  Burning  surface  profiles. 

(a)  Typical  profile  at  pressures  above  6  MPa  (AP-HC  propellants) 

(b)  Surface  profile  typical  of  analytical  models 

(c)  Burning  surface  after  quench  by  rapid  depressurization 
(UTP  3001  propellant,  AP-HC-Al) 
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products  flushed  by  the  control  gas  flow.  Visual  observations  were  re¬ 
corded.  Five  different  hydrocarbon  binders  were  used. 

A  typical  test  required  about  60  seconds  to  reach  600°  C  and  was  con¬ 
tinued  to  1000  to  1400°  C.  All  binders  showed  some  degree  of  internal 
gasification  with  a  boiling-like  behavior  indicative  of  fluidity.  With 
increasing  temperature,  this  evolves  into  a  tar-like  (usually  discolored) 
substance  of  less  volume  than  the  original  sample,  which  vaporizes  as  tem¬ 
perature  increases  further,  leaving  a  modest  residue  of  char-like  material. 
In  an  oxygen  atmosphere  the  gasifying  tar  ordinarily  ignited  at  about  600° 
C,  while  the  one  binder  that  gasified  at  lower  temperature  did  not  ignite. 
In  the  fluid  stage,  the  different  binders  showed  significantly  different 
degrees  of  fluidity.  Two  binders  yielded  a  small  amount  of  thermally 
stable,  clear  liquid  at  1000°  C,  but  only  when  heated  in  oxygen.  Details 
of  the  results  were  reported  in  Reference  5,  and  are  summarized  in  Figure 
2.  The  binders  used  in  these  tests  were  used  also  in  other  experiments 
sunmarized  in  later  parts  of  this  report,  and  imply  that  surface  temper¬ 
atures  of  the  rapidly  heated  binder  in  the  propellant  combustion  zone  may 
be  as  high  as  600+  degrees  C,  with  gas  evolution  and  melt  properties  at 
much  lower  temperatures.  Polysulfide  binder  decomposed  at  the  lowest 
temperature  with  the  least  melt;  polyurethane  produced  the  most  melt;  PBAN 
and  HTPB  binders  produced  the  most  extensive  melts,  with  PBAN  melt  being 
the  most  fluid.  The  clear  liquid  residue  formed  by  HTPB  and  CTPB  in  oxygen 
atmospheres  still  persisted  (although  in  diminished  amount)  at  1400°C. 
Collectively,  the  results  suggest  that  the  fluid  state  of  the  binder  or 
its  intermediate  products  will  ordinarily  be  present  on  the  propellant 
burning  surface  to  high  temperatures,  with  polysulfide  being  an  exception. 


SANDWICH  BURNING 


3.1  Introduction 

If  one  is  seeking  understanding  of  combustion  or  a  comprehensive  basis 
for  its  control,  a  composite  propellant  is  a  remarkably  difficult  object 
for  study.  Not  only  is  it  very  costly  to  make  up  samples  with  systematic 
changes  in  relevant  variables,  but  the  resulting  samples  are  of  chaotic 
character  on  a  microscopic  (  1  mm)  scale.  Given  the  unavoidable  difficul¬ 

ties  of  producing,  and  making  time-resolved  observations  of  a  microscopic 
subject  in  a  hostile  environment,  the  added  problems  of  sample  preparation 
and  characterization  have  severely  limited  progress  in  research,  and  it 
is  not  surprising  that  investigators  have  resorted  to  other,  more  control¬ 
lable  and/or  inexpensive  schemes  for  preparing  test  samples.  One  of  these 
is  the  ingredient  "sandwich"  illustrated  in  Figure  3.  Alternate  sheets 
of  propellant  ingredients  are  laminated  together,  and  burned  edgewise. 

Such  a  configuration  has  the  advantages  of 

1.  Low  cost 

2.  Control  over  microstructure  and  placement  of  ingredients 

3.  Relative  ease  of  observation  during  burning 

4.  Ordered  structure  that  qualitatively  establishes  the  structure 
of  the  combustion  zone 

5.  Relative  amenability  of  processes  to  analytical  modeling  (2-dimen¬ 
sional  steady  state  behavior) 

It  is  these  attributes  that  have  motivated  the  use  of  sandwiches  in  past 
and  present  investigations  of  propellant  combustion  mechanisms. 

3.2  Combustion  Zone  of  AP/HC  Sandwiches 

The  details  of  sandwich  combustion  are  suggested  by  Figure  4.  Sand¬ 
wiches  are  ordinarily  ignited  on  the  top  edge,  and  side  burning  is  avoided 
by  a  modest  flow  of  nitrogen,  which  also  minimizes  recirculation  of  com¬ 
bustion  products  in  the  line  of  sight  of  combustion  photography.  Above 
its  self-deflagration  pressure  limit,  AP  has  an  exothermic  decomposition 
flame  ( (2) ,  in  Figure  4),  and  a  thick  AP  lamina  will  burn  independently 
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Figure  3.  Sketch  of  AP-binder-AP  sandwiches 

(a)  Conventional 

(b)  Tapered 
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Figure  4.  The  combustion  zone  structure  of  an  AP-HC  sandwich 
(profile  typical  of  pressure  5-10  MPa). 


20 


me*-..  ., 


of  the  binder  (Figure  5a).  Near  the  binder-AP  interface,  a  diffusion  flame 
occurs  between  oxidizer  and  binder  vapors  ((2),  in  Figure  4).  The 
binder  lamina,  which  decomposes  endothermally  under  the  influence  of  the 
diffusion  flame,  tends  to  protrude  and  may  even  divide  the  combustion  region 
into  two  halves  (especially  if  conditions  are  unfavorable  for  pyrolysis 
of  the  binder).  Anything  that  selectively  enhances  the  rate  of  the  AP- 
binder  reaction  tends  to  accelerate  the  surface  regression  in  the  vicinity 
of  the  interface,  giving  a  profile  like  that  in  Figure  5b.  At  pressure 
below  the  AP  deflagration  limit,  the  AP  pyrolysis,  like  the  binder,  is 
supported  by  the  diffusion  flame.  In  this  case  the  AP  burning  surface 
lags  the  interfacial  region  just  as  the  binder  does  (Figure  5c).  Depending 
on  details  of  the  situation,  the  burning  may  simply  proceed  down  and  leave 
a  slot,  or  the  surface  will  develop  into  a  deep  V-shaped  contour  (see  later). 

This  simplistic  description  of  sandwich  burning  is  consistent  with 
current  understanding  of  the  decomposition-deflagration  characteristics 
of  the  binder  and  AP,  and  with  most  of  the  current  literature.  However, 
there  are  difficulties  on  both  the  microscopic  and  macroscopic  scale  that 
must  be  resolved  to  fully  understand  sandwich  burning  and  its  relation 
to  combustion  of  solid  propellants.  These  difficulties  are  addressed  in 
the  following. 

3.3  Classification  of  Sandwiches  and  Expected  Behavior 

As  one  may  observe  in  Figure  4,  sandwich  burning  is  characterized 
by  several  exothermic  reaction  steps  that  "drive"  the  combustion.  AP  oxi¬ 
dizer  burns  as  a  monopropellant  (above  about  20  atm),  with  exothermic 
reaction  at  the  surface,  and  in  a  decomposition  flame  in  the  gas  phase. 

These  heat  sources  are  augmented  by  the  diffusion  flame  between  binder 
and  oxidizer  vapors.  The  locations  of  these  heat  sources  were  denoted 
in  Figure  4  by  the  number  inserts  1-4.  The  way  the  sample  burns  is  de¬ 
termined  by  the  return  of  heat  to  the  solid  surface  to  support  continued 
pyrolysis.  This  in  turn  depends  upon  the  location  of  the  exothermic  re¬ 
action  sites,  the  relative  amount  of  heat  release  in  each,  and  the  overall 
stoichiometry  of  the  sample  (amount  of  heat  release  compared  to  amount 


21 


Figure  5.  Surface  profiles  of  burning  sandwiches  of  AP  and  HC  binders. 

(a)  AP  deflagration  dominating  deflagration  (typical  of  elevated 
pressure,  e.g.  10  MPa) 

(b)  AP  deflagrating,  but  regression  rate  enhanced  by  diffusion 
flame  contribution  (typical  of  medium  pressures,  e.g.  4  MPa) 

(c)  AP  won't  deflagrate  without  heat  input  from  the  diffusion 
flame  (typical  of  low  pressure,  e.g.  <  2  MPa) 

of  heat  required  for  heat-up  and  pyrolysis  of  the  solid). 

The  configuration  of  the  burning  surface  of  sandwiches  has  been  some¬ 
what  stereotyped  in  the  past  according  to  Figure  5a  or  Figure  5c,  where 
5a  is  typical  for  self-deflagrating  oxidizers  and  Figure  5c  is  for  oxidizers 
that  require  heat  from  the  diffusion  flame  to  maintain  gasification  (this 
includes  AP  at  low  pressures).  In  terras  of  propellant  combustion,  these 
classifications  are  superfluous,  because  they  refer  to  samples  with  thick 
binder  and  oxidizer  layers  that  are  not  typical  of  propellants,  and  they 
describe  the  overall  appearance  of  those  samples.  At  the  very  least,  con¬ 
sideration  should  be  directed  to  the  details  of  surface  geometry  near  the 
oxidizer-binder  interface,  and  in  that  region,  consideration  needs  to  be 
given  to  whether  a  thick  binder  layer  is  causing  flooding  by  molten  binder, 
and  to  the  localized  effect  of  the  diffusion  flame  near  the  interface. 

In  the  interest  of  more  objective  observation  and  reporting  of  present 
and  future  results,  a  scheme  of  classification  of  sandwiches  is  proposed 
in  Figure  6.  In  this  categorization,  all  samples  are  first  divided  into 
two  groups  according  to  whether  the  binder  is  thick  or  thin,  where  "thin” 


tH 

CM 

CO 

H 

CM 

CO 

i— 1 

CM 

co 

H 

CM 

CO 

i— 1 

CM 

CO 

H 

CM 

CO 

i—l 

CM 

CO 

i-4 

CM 

CO 

H 

H 

i— 1 

CM 

CM 

CM 

H 

»— 1 

H 

CM 

CM 

CM 

H 

H 

H 

CM 

CM 

CM 

i— 1 

»-4 

i—l 

CM 

CM 

CM 

1-4 

H 

1-4 

H 

r-4 

H 

CM 

CM 

CM 

CM 

CM 

CM 

H 

i—4 

r-4 

H 

H 

i—l 

CM 

CM 

CM 

CM 

CM 

CM 

.-4 

i— 1 

H 

1-4 

H 

H 

i—l 

H 

H 

1—* 

H 

r4 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

n. 0.0.0  Binder  Thickness 
O.n.O.O  AP  Thickness 

O.O.n.O  Pressure  (relative  to  AP  deflagration  limit) 
O.O.O.n  Binder  Pyrolysis  Characteristics 


0.0.0. 1  Heat  resistant,  nonmelting  binder 
0.0.0. 2  Binder  with  thermally  stable  melt 
0.0. 0.3  Easily  gasified  binder,  minimal  melt 


Figure  6. 


Classification  of  sandwich  burning  experiments. 


Chart  designed  for  self-deflagrating  oxidizer,  with  "low  pressure"  refer¬ 
ring  to  below  the  self-deflagration  pressure  and  "high  pressure"  to  con¬ 
ditions  above  the  limit.  For  the  time  being,  "thick"  binder  refers  to 
thickness  compared  to  that  region  of  the  surface  in  which  presence  of  the 
binder  affects  the  regression  front  of  the  oxidizer  at  pressures  above 
the  deflagration  limit  of  the  oxidizer. 


means  comparable  Co  structural  elements  of  binder  in  propellants  (e.g., 

2-20  p,m).  These  categories  1  and  2  are  each  divided  into  two  more  depending 
on  whether  the  oxidizer  layer  is  thick  (1.1  and  2.1)  or  thin  (1.2  and  2.2). 
In  the  present  work,  the  terms  "thick"  and  "thin"  oxidizer  refer  to  overall 
sample  stoichiometry,  with  thick  oxidizer  meaning  oxidizer-rich.  The  third 
categorization  (third  digit  of  the  classification  number)  refers  to  whether 
the  oxidizer  self-def lagrates  or  not  (In  the  classification  diagram  this 
is  referred  to  as  high  and  low  pressure,  because  we  were  concerned  with 
AP  oxidizer).  The  fourth  categorization  (see  diagram  and  4th  digit)  refers 
to  the  melt-gasification  characteristics  of  the  binder.  Using  this  cate¬ 
gorization  of  sandwiches,  their  burning  behavior  can  be  more  comprehensive¬ 
ly  categorized  as  in  the  following. 

Looking  first  at  sandwiches  with  thick  binder  (categories  1.0. 0.0, 
Figure  7,  top),  the  binder  tends  to  protrude  and  divide  the  primary  part 
of  the  diffusion  flame  into  two  flames  that  proceed  relatively  independently 
of  each  other.  These  flames  will  extend  beyond  the  binder  and  merge  if 
the  overall  stoichiometry  of  the  sample  is  not  too  fuel-rich  (thick  AP, 
case  1.1. 0.0).  If  the  stoichiometry  is  fuel-rich,  (case  1.2. 0.0),  the 
flames  will  not  merge  and  the  fuel  protrusion  will  increase  progressively 
during  burning.  If  the  pressure  is  high  (1.1. 0.0,  1.2. 1.0),  the  oxidizer 
will  recede  ahead  of  the  fuel  protrusion  (Figure  7),  while  low  pressure 
(1.1. 2.0,  1.2. 2.0)  will  lead  to  protrusion  of  the  oxidizer  as  well  as  the 
fuel.  The  extent  of  the  AP  protrusion  will  be  less  with  thin  AP  (fuel- 
rich  stoichiometry,  1.2. 2.0). 

When  the  binder  is  thin  (lower  row,  Figure  7),  it  does  not  protrude 
because  it  is  so  easily  pyrolized.  Accordingly  the  two  diffusion  flames 
are  not  distinguishable  (most  work  reported  in  previous  literature  does 
not  correspond  to  this  thin  binder  situation,  although  propellant  combustion 
does  along  most  of  the  exposed  oxidizer-binder  interfaces  on  the  burning 
surface).  At  high  pressure  (2. 1.1.0  and  2. 2.1.0)  the  AP  recedes  with  the 
binder,  but  at  low  pressure  the  AP  recedes  less  rapidly  and  protrudes 
(2. 1.2.0  and  2. 2. 2.0).  In  the  Chick  oxidizer  case  (2. 1.2.0),  the  oxidizer- 
binder  interface  region  burns  down  to  form  a  groove,  the  fuel  being  insuf¬ 
ficient  to  react  with  the  oxidizer  excess  (samples  of  this  type  often  will 
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Figure  7.  Burning  surface  profiles  corresponding  to  the  classes  of 
sandwiches  in  Figure  6. 


not  sustain  burning).  With  thinner  oxidizer  and  low  pressure  (2. 2. 2.0), 
the  oxidizer  protrusion  is  present  but  not  so  conspicuous  because  the  sam¬ 
ple  is  nearer  to  stoichiometric  overall. 

Much  of  the  results  of  early  work  on  sandwich  burning  (of  AP/HC  sam¬ 
ples)  can  be  anticipated  on  the  basis  of  the  foregoing  classification  scheme, 
and  some  apparent  contradictions  can  be  rationalized  in  that  framework. 

As  experimental  methods  are  improved,  many  more  detailed  aspects  of  combus¬ 
tion  behavior  become  evident,  such  as  the  observation  of  characteristic 
surface  patterns  on  the  AP  with  evidence  of  reacting  surface  melts;  flow 
of  binder  melt  onto  the  AP  surface;  accumulation  of  ballistic  modifiers 
on  the  AP  surface;  and  non- stationary  (intermittent)  behavior  of  the  diffu¬ 
sion  flames.  Perhaps  most  significant,  the  behavior  in  the  vicinity  of 
thin  binder  laminae  (<50M,m)  has  proven  to  be  very  complex.  Since  this 
is  the  dimensional  domain  relevant  to  propellant  combustion,  the  present 
research  has  gradually  turned  to  that  domain.  This  has  proven  to  be  a 
much  more  difficult  domain  to  study. 

3.4  Scope  of  Current  Investigations 

The  objective  of  sandwich  burning  in  the  present  research  was  to  deter¬ 
mine  the  effect  of  binder  thickness,  type  of  binder,  and  pressure  on  the 
details  of  combustion  zone  behavior.  In  addition  to  determining  burning 
surface  profiles,  the  investigations  sought  to  determine  the  effect  of 
binder  melt  flow,  and  the  unsteady  nature  of  the  diffusion  flame.  Experi¬ 
mental  methods  consisted  primarily  of  high  speed  motion  picture  photography, 
and  interrupted  burns  followed  by  microscopic  examination.  Binder  thick¬ 
ness  effects  were  investigated  by  burning  the  tapered  edge  of  "tapered 
sandwiches,"  in  which  the  thickness  of  the  binder  (as  viewed  on  the  burning 
edge)  was  graduated  from  about  10  ^.m  at  one  edge  to  about  130  y,m  on  the 
opposite  edge  (all  samples  were  of  this  form  unless  otherwise  specified). 

The  combinations  of  kind  of  binder  and  combustion  pressure  are  shown  by 
the  test  matrix  in  Table  1.  The  shaded  boxes  in  the  matrix  indicate  quench 
tests  in  which  results  were  not  obtained  or  were  compromised  by  experimental 
difficulties.  The  cross-hatched  boxes  indicate  conditions  that  were  not 
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Table  1.  Tests  on  Tapered  Sandwiches 


tested.  For  each  pressure  there  is  a  column  for  indicating  quench  tests 
and  a  column  for  indicating  photographic  tests.  From  the  test  results, 
observations  were  made  on  a  number  of  basic  aspects  of  the  combustion  zone. 
The  results  are  described  in  the  context  of  these  details. 

3.5  Interface  Profile  and  Melt  Flows 

Sandwich  burning  studies  reported  in  the  early  1960's  stated  that 
the  leading  edge  of  the  burning  front  was  at  the  binder-oxidizer  interface 
(Figure  8a),  a  result  that  was  used  by  some  investigators  to  support  an 
argument  that  burning  rate  was  dominated  by  heterogeneous  reactions  at 
the  interface.  However  more  careful  experiments  showed  (Reference  6)  that 
the  leading  edge  was  never  at  the  interface  (AP  oxidizer),  but  was  usually 
on  the  oxidizer  surface  not  far  from  the  interface  (e.g.  50-500  p,m)  (Figure 
8b).  It  was  assumed  that  the  binder,  being  a  heat  sink,  lowered  the  temp¬ 
erature  of  the  AP  near  the  interface  enough  to  retard  its  deflagration. 
However  the  diffusion  flame  preferentially  heats  this  same  region,  so  that 
no  presumption  can  be  made  about  dominance  of  these  two  opposing  effects. 

An  analytical  model  of  combustion  in  the  interface  region  (Reference 
7)  indicated  that  the  oxidizer  surface  should  regress  uniformly  to  within 
a  few  microns  of  the  interface,  at  which  point  the  AP  surface  would  lead 
slightly  due  to  heat  transfer  from  the  diffusion  flame  via  the  binder  (Fig¬ 
ure  8c).  Since  this  was  contrary  to  the  observed  results,  it  was  specula¬ 
ted  that  the  binder  melt  flow  (observed  in  many  quench  tests)  was  retarding 
the  AP  deflagration  at  the  interface.  Accordingly,  samples  were  made  in 
which  the  oxidizer  lamina  was  formed  by  dry-pressing  AP  powder  on  sheets 
of  concrete  and  of  mica.  The  samples  were  burned  and  quenched  (Table  1) 
at  burning  pressures  high  enough  to  burn  without  fuel  (6.9  MPa,  1000  psi). 
The  quenched  samples  with  concrete  sheets  in  place  of  binder  crumbled  dur¬ 
ing  burning;  the  samples  with  mica  sheets  burned  normally  and  gave  good 
quench  samples  (except  for  some  delamination  within  the  mica).  The  quench¬ 
ed  samples  exhibited  surface  profiles  that  were  flat  all  the  way  to  the 
interface,  except  for  the  usual  formation  of  surface  patterns  over  the  en¬ 
tire  AP  surface  (Figure  9).  Since  the  typical  dimension  of  irregularity 


Figure  8.  Different  leading  edge-interface  profiles  postulated  in 
earlier  work. 

(a)  Early  descriptions  of  profile,  probably  relevant  with 
endothermic  oxidizers  (e.g.,  KP) 

(b)  Generally  accepted  (1979)  profile  for  AP 

(c)  Prediction  of  2-D  model  without  binder  encroachment 
on  AP  surface 
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Figure  9.  Surface  of  quenched  AP-mica  sandwich  burned  at  6.9  MPa. 

Note  AP  surface  detail  ’-eriains  typical  all  the  way  to  the 
mica  interface  (a  delaminated  flake  of  mica  obstructs  the 
interface  at  upper  end). 
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Figure  10.  Profiles  of  quenched  surfaces  of  sandwiches,  illustrating 

details  near  the  interface. 

(a)  Typical  of  AP-HTPB  at  4.1  MPa  (binder  thickness  >  70  tun) 

(b)  Typical  of  AP-PBAN  binder  at  13.8  MPa  (binder  thickness 
>  70  urn) 

(c)  Typical  of  AP-HTPB  binder  at  4.1  MPa  (binder  thickness 
~  50  gjn) 

(d)  Typical  of  AP-all  binders  at  1.4  MPa  (binder  thickness 
<  25  p,ra) 

of  the  surface  patterns  is  comparable  to  the  dip  in  the  AP  surface  at  the 
interface  predicted  by  the  two-dimensional  model,  it  was  not  possible  to 
verify  this  prediction  of  the  model.  However  both  this  result  and  the 
model  imply  that  the  typical  retarded  regression  near  the  interface  in 
AP-binder  sandwiches  is  due  to  binder  melt  flow  or  other  influence  of  the 
binder.  It  should  be  noted  that  melt  flow  (as  judged  from  quenched  samples) 
occurs  on  a  dimensional  scale  of  the  same  order  as  oxidizer  particle  sizes 
in  propellants,  raising  the  question  of  the  possible  importance  of  melt 
floj®  in  propellant  combustion.  Accordingly,  it  also  suggests  that  the 
melt  characteristics  of  the  binder  are  important.  Investigation  of  these 
questions  is  continuing. 

Recent  results  with  tapered,  thin-binder  sandwiches  have  revealed 
a  variety  of  profiles  adjoining  the  interface  (illustrated  by  sketch  in 
Figure  10).  These  will  be  discussed  in  appropriate  context  in  following 


30 


sections.  However,  the  following  generalizations  can  be  made: 

1.  With  the  relatively  thick  AP  laminae  used  in  the  present  work, 
most  of  the  quenched  AP  surface  looks  like  that  of  AP  samples 
tested  without  binder  laminae  (i.e.,  above  the  AP  low  pressure 
deflagration  limit).  The  surface  exhibits  characteristic  undu¬ 
lations,  and  froth-like  residue  (reported  in  many  references, 
e.g.  Reference  8). 

2.  There  is  always  a  region  of  the  AP  surface  adjoining  the  inter¬ 
face  where  the  surface  is  relatively  smooth  (Figure  11,  12;  see 
also  later  figures). 

3.  The  smooth  surface  often  exhibits  obvious  flow  patterns,  espe¬ 
cially  near  and  across  the  interface  (See  later  figures.). 

4.  The  surface  effects  extend  for  10-200  p,m  from  the  interface, 
sometimes  further. 

5.  Regression  of  the  AP  surface  is  retarded  in  the  smooth- surf ace 
region  (Figure  10,  11).  As  a  corollary,  the  surface  regression 
in  the  AP  is  higher  in  the  region  just  beyond  the  smooth  region. 
The  regression  is  usually  highest  there,  suggesting  the  impor¬ 
tance  of  the  AP-binder  flame  on  burning  rate  (below  about  6.9 
MPa,  1000  psi). 

6.  The  surface  effects  2-5  are  absent  when  the  binder  is  replaced 
by  mica  (Figure  9). 

It  is  premature  to  say  whether  the  smooth-surface  band  adjoining  the 
interface  is  due  simply  to  melt  flow,  although  the  evidence  is  very  strong 
that  part  of  the  band  nearest  the  binder  lamina  originates  as  flow  from 
the  binder.  The  outer  portion  of  the  band  may  reflect  the  (as  yet  undeter¬ 
mined)  process  by  which  binder-coated  AP  decomposes,  or  may  reflect  AP 
decomposition  modified  by  diffusion  of  gaseous  species  from  the  binder 
or  the  leading  edge  of  the  diffusion  flame.  Given  the  fact  that  the  lead¬ 
ing  edge  of  the  AP  deflagration  front  is  apparently  caused  by  the  added 
heating  resulting  from  thermal  diffusion  from  the  oxidizer-binder  gas 
flame,  it  is  reasonable  to  postulate  that  chemical  species  from  that  flame 
can  also  diffuse  a  comparable  distance  and  modify  surface  processes. 

*  The  meaning  of  "highest"  becomes  unclear  in  profiles  such  as  Figure 

10  c,  d. 


Surface  of  quenched  sandwich  showing  effect  of  low  pressure  and  binder  thickness 
on  the  burning  surface  profile.  This  is  an  example  of  an  AP-PS  sandwich  partial 
burned  at  1.4  MPa. 


Perhaps  the  most  important  point  merits  repetition,  i.e.,  that  relative 
to  real  rocket  propellants,  the  surface  distances  involved  in  the  smooth- 
band  phenomenon  are  of  the  same  order  as  oxidizer  particle  size  in  pro¬ 
pellants.  This  implies  that  deflagration  of  oxidizer  particles  may  be 
materially  affected  by  binder  behavior  and  diffusion  flames,  and  in  a  way 
dependent  on  distance  from  the  particle  periphery  and  possibly  in  a  fluc¬ 
tuating,  time  dependent  way.  The  implication  of  these  observations  re¬ 
lative  to  effects  of  particle  size  on  steady  and  oscillatory  burning  should 
be  noted.  It  is  also  worthy  of  note  that  the  retardation  effect  on  burning 
rate  can  be  quite  large  at  high  pressure  (Section  3.7,  with  profiles  as 
in  Figure  10b). 

3.6  Binder  Thickness 

As  suggested  in  the  earlier  classification  system  for  sandwiches 
(Figure  6),  binder  thickness  is  an  important  factor  in  how  sandwiches  burn. 
Most  investigators  have  used  thicknesses  that  are  an  order  of  magnitude 
larger  than  the  thicknesses  in  high-solids  heterogeneous  propellants. 

This  has  the  effect  of  placing  most  of  the  diffusion  flame  at  a  relatively 
large  distance  from  the  burning  surface,  excepting  that  region  close  to 
the  oxidizer-binder  interface.  Accordingly,  it  is  probably  worth  repeating 
that  only  that  region  of  the  surface  located  near  the  interface  experiences 
propellant-like  conditions.  Further,  even  that  region  may  be  influenced 
by  behavior  further  out  from  the  interface  that  is  not  typical  of  propel¬ 
lants  (for  example,  by  abnormally  large  melt  flow  from  a  thick  binder 
layer) . 

In  the  present  work,  the  effect  of  binder  thickness  was  investigated 
by  testing  "tapered"  sandwiches  (Figure  3)  under  the  conditions  in  the 
test  matrix  of  Table  1.  The  general  trend  of  behavior  is  in  accordance 
with  class  1.1. 1.0  and  1.1. 2.0  samples  (Figure  6,  7:  thick  binder-thick 
AP-high  and  low  pressure),  with  a  transition  to  class  2. 1.1.0  and  2. 1.2.0 
samples  (thin  binder-thick  AP-high  and  low  pressure)  towards  the  thin- 
binder  end  of  the  sample. 

The  thick  binder  (100-150  y,  ra)  end  of  the  samples  yielded  nearly  flat 
AP  with  binder  protruding  in  varying  degree  with  evidence  of  binder  melt 
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SECTION  B-B 


SECTION  C-C 


Figure  13.  General  trend  of  surface  profile  with  binder  thickness. 

flow  (high  pressure  tests,  e.g.,  Figure  lie,  category  as  in  Figure  7  - 
l.l.l.O);  or  burning  into  a  deep,  relatively  flat  bottom  groove  in  the 
AP  with  minimal  binder  projection  (low  pressure  tests,  e.g.,  Figure  12, 
category  as  in  Figures  6  and  7,  1.1.2. 2,  1.1. 2. 3). 

The  thin  binder  (10-15nm)  end  of  the  samples  yielded  nearly  flat 
AP  with  binder  flat  or  recessed  (high  pressure  tests,  e.g.,  Figure  11, 
category  as  in  Figures  6,  7  -  2. 1.1.0);  or  burning  into  a  deep  groove  in 
the  AP  with  binder  recessed  sharply  below  the  AP  (low  pressure  tests,  e.g., 
Figure  12,  category  as  in  Figures  6,  7  -  2. 1.2.0,  and  Figure  lOd).  Burning 
would  not  sustain  at  low  pressure  with  very  thin  binder. 

To  the  extent  that  generalizations  can  be  made,  the  trends  with  binder 
thickness  are  summarized  by  the  sketch  in  Figure  13  (which  is  drawn  to 
correspond  to  a  pressure  of  4-5  MPa).  "Thick"  binder  tends  to  protrude 


and  give  a  melt  flow  across  the  binder/  AP  interface,  the  extent  of  these 
features  depending  on  binder  characteristics  and  pressure.  As  one  proceeds 
to  thinner  binder,  the  protrusion  and  flow  are  less  conspicuous,  and  re¬ 
gions  of  the  binder  profile  sometimes  extend  below  the  interface  AP.  With 
binder  thickness  less  than  about  30  (im  (or  greater  thickness  at  low  pres¬ 
sures),  the  binder  is  completely  recessed  below  the  adjoining  AP.  As  noted 
in  Section  3.5,  the  AP  profile  is  inclined  upward  toward  the  interface 
plane,  and  this  was  true  at  all  binder  thicknesses  and  pressures  tested. 

This  up-sloping  surface  always  exhibited  relative  freedom  from  the  frothy 
characteristic  of  the  rest  of  the  AP  surface  (Figures  11,  12,  14-17),  often 
showing  initial  markings  suggestive  of  melt  flow  from  the  binder. 

The  general  trends  with  lamina  dimensions  and  pressure  were  in  accord 
with  the  arguments  leading  to  Figure  7,  provided  one  makes  allowances  for 
the  binder  properties  (the  fourth  digit  in  the  sandwich  categorization 
scheme  in  Figure  6).  However,  there  are  some  features  of  the  behavior 
that  go  beyond  the  basis  of  Figures  6  and  7.  Most  notable  was  the  demon¬ 
stration  of  a  binder  thickness-dependent  low  pressure  deflagration  limit 
(Figures  12,  14).  At  1.4  MPa,  sandwiches  were  difficult  to  ignite,  and 
did  not  burn  the  thin-binder  end  of  the  sample  (that  portion  with  binder 
thickness  less  than  25  p,m) .  This  pressure  is  below  the  deflagration  limit 
of  the  AP  alone,  and  the  diffusion  flame  apparently  cannot  supply  enough  heat 
to  sustain  losses  and  pyrolysis  of  the  AP  and  binder.  Interpreted  in  the  con¬ 
text  of  a  propellant  surface,  this  result  provides  some  understanding  of 
the  typically  high  dependence  of  burning  rate  on  pressure  and  particle  size 
at  low  pressure,  and  the  high  dynamic  response  during  oscillatory  flow. 

The  flamelets  are  not  only  unstable  individually,  but  also  collectively. 

3.7  Comparison  of  Binders 

As  noted  in  Section  2,  the  propellant  binders  used  in  this  program 
have  observable  differences  in  thermal  decomposition  behavior,  with  the 
polysulfide  binder  being  most  distinctive  because  of  its  low  decomposition 
temperature  (approximately  200°  C  lower  than  the  others,  Figure  2). 

As  noted  in  the  test  matrix,  tapered  sandwich  tests  were  run  with  these 
binders  at  several  pressures.  The  collected  test  results  show  significant 
differences. 
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Surface  of  quenched  sandwich  illustrating  limitation  of  flammabili 
The  thin-binder  end  has  not  sustained  burning  (HTPB  binder,  test  a 


Figure  15a  and  b  compare  the  quench  surface  of  a  sandwich  with  PS 
binder  and  one  with  HTPB  binder  (both  were  burned  at  4.1  MPa,  600  psi). 

The  figure  illustrates  the  pervasive  effect  of  the  binder  on  the  surface 
of  the  adjoining  AP,  an  effect  that  is  present  wi th  al 1  binders  and  all 
binder  thicknesses  tested  (although  less  pronounced  with  thin  binder). 

There  is  less  appearance  of  melt  flow  in  the  case  of  PS  binder,  less  total 
protrusion  of  the  binder  at  the  thick  end  of  the  lamina,  and  the  binder 
is  more  deeply  recessed  in  the  groove  at  the  thin  end. 

A  singular  behavior  observed  with  PBAN  binder  (and  superficially  with 
HTPB)  only  at  high  pressure  (13.8  MPa,  2000  psi)  is  shown  in  Figure  16 
and  corresponds  to  the  sketch  in  Figure  10b.  The  binder  appears  to  se¬ 
verely  retard  or  arrest  the  AP  deflagration  in  the  "smooth"  band  out  to 
about  100  m  from  the  interface.  At  this  pressure  the  balance  of  the  AP 
surface  regresses  ahead  of  the  retarded  surface,  leaving  a  plateau-like 
raised  section  of  AP  and  binder.  This  raised  section  apparently  becomes 
undercut  by  lateral  AP  deflagration,  as  the  plateau  is  discontinuous  along 
the  interface  in  an  irregular  manner.  The  occurrence  of  this  behavior 
was  not  conspicuously  dependent  on  binder  thickness,  except  that  with  thin 
binder  the  regarded  regions  of  the  AP  adjoining  the  binder  simply  slope 
away  from  the  interfaces,  with  the  resulting  ridge  of  AP  being  topped  by 
a  binder  recess  (profile  category  2.2.1  of  Figure  7).  One  may  speculate 
that  binder  melt  flow  formed  a  stable  residue  when  it  decomposed,  one  that 
protects  the  underlying  AP  from  heat  and  chemical  attack.  It's  not  clear 
why  this  occurs  only  in  high  pressure  tests,  and  only  with  selected  binders 
(generalizations  made  hesitantly  because  of  limited  testing),  but  it  may 
be  significant  that  AP  itself  experiences  a  transition  in  mode  of  deflagra¬ 
tion  in  the  12.4-15.1  MPa  (1800-2200  psi)  pressure  range  (Reference  9,10). 

The  retarded  region,  as  at  all  pressures,  is  also  coincident  with  the  smooth 
band. 

Differences  in  binder  behavior  were  most  conspicuous  in  the  low  pres¬ 
sure  tests  (1.4  MPa;  200  psi).  In  all  cases  there  appeared  to  be  a  thresh¬ 
old  binder  thickness  for  the  binder  below  which  burning  would  not  self- 
sustain.  This  was  illustrated  by  Figure  13,  which  shows  the  quench  surface  of 
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a  tapered  sandwich  that  burned  normally  on  the  thick  end.  In  that  portion 
of  the  sample  that  burned,  the  appearance  of  the  various  binders  differed. 

PS  binder  was  most  deeply  recessed  and  showed  a  residue  of  unidentified 
flakey  material  (Figure  17a).  HTPB  and  CTPB  binders  formed  miniscus-like 
profiles  in  the  recess,  with  some  visible  flow  marks  (Figure  17b).  PBAN 
binder  showed  a  complex  stringy  structure  in  the  recess  (Figure  17c)  sug¬ 
gestive  of  viscous  properties  (it's  not  clear  what  process  caused  formation 
of  the  structure,  but  it  may  be  the  result  of  subsurface  bubbling  during 
burning  or  during  the  quench  depressurization).  It  is  a  little  surprising 
that,  even  with  the  recessed  binder  situation,  all  the  binders  produced 
the  smooth  band  on  the  adjoining  AP  surface  outside  the  recess,  and  the 
point  of  maximum  regression  of  AP  is  always  at  the  outer  edge  of  the  smooth 
band.  Conditions  for  melt  flow  over  the  AP  surface  appear  to  be  very  un¬ 
favorable  in  the  recessed,  thin  binder  situation,  suggesting  that  binder 
effects  on  the  AP  surface  are  not  necessarily  due  only  to  melt  flow. 

3.8  Diffusion  Flame 

High  speed  combustion  photography  shows  luminous  flames  originating 
near  the  oxidizer-binder  interfaces.  These  flames  are  orange  in  color 
and  very  likely  are  plumes  of  radiating  hot  carbon  in  the  diffusion  flame. 

The  nature  of  these  flames  has  been  noted  before  (References  11,  12),  ear¬ 
lier  observations  having  usually  been  made  by  viewing  end-on  to  the  sample 
so  as  to  see  in  between  the  walls  of  AP  that  occur  in  low  pressure  tests. 

The  flame  seen  from  this  orientation  was  referred  to  as  "brushy"  (Reference 

11). 

With  tapered  sandwiches,  there  is  a  motivation  to  view  the  diffusion 
flame  sheet  side-on,  so  as  to  resolve  the  flame  behavior  as  a  function 
of  binder  thickness.  Initial  tests  showed  that  this  arrangement  was  feasible, 
and  that  significant  differences  in  flame  behavior  did  indeed  occur  as 
a  function  of  binder  thickness.  However  the  diffusion  flame  is  obscured 
from  side  viewing  in  low  pressure  tests  by  the  protruding  wall  of  slower- 
burning  oxidizer,  limiting  useful  results  to  tests  above  4.1  MPa  (600  psi) 
or  so. 

Figure  18  shows  a  frame  of  a  side-view  motion  picture  of  an  HTPB  sandwich 
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Singular  features  of  three  binders  at  1.4  MPa. 

(a)  Polysulfide  binder,  thickness  22  um  (b)  HTPB  binder,  thickness 
(c)  PBAN  binder,  thickness  26  y,ra 


Figure  18.  Diffusion  flame lets  on  a  burning  sandwich.  Binder  thickness 
decreases  to  the  right.  HTPB  binder,  test  pressure  6.9  MPa. 
Picture  shows  full  width  of  the  sample,  about  1  cm. 

while  burning  at  6.89  MPa  (1000  psi).  This  corresponds  to  an  exposure 
time  of  about  1/3000  sec,  and  shows  the  diffusion  flame  to  be  a  line  array 
of  separate  flamelets.  These  flamelets  appear  to  be  transitory  in  nature, 
on  a  very  short  time  scale.  Most  are  not  identifiable  in  successive  frames 
of  motion  pictures  taken  at  2500  frames  per  second  (a  few  flamelets  at 
the  thick-binder  end  of  the  interface  are  recognizable  for  2  or  3  frames). 
Many  motion  pictures  were  made  of  sandwiches  (Table  1)  and  examined  for 
flamelet  behavior.  The  transient  behavior  does  not  appear  to  be  an  oscil¬ 
lation  in  position,  but  rather  a  truly  transitory  existence.  The  average 
distance  between  flamelets  is  greater  in  the  thick-binder  end  of  the  inter¬ 
face  as  is  the  size  of  the  flamelet.  With  HTPB  binder  there  were  usually 
two  rows  of  flamelets,  one  for  each  interface.  Indeed,  there  is  some  indi¬ 
cation  that  the  distinctness  of  the  two  rows  of  flamelets  with  HTPB  binder 
resulted  from  a  tendency  for  the  two  interfaces  to  regress  one  ahead  of 
the  other  with  that  particular  binder.  Resolution  in  pictures  does  not 
establish  whether  the  two  rows  remain  identifiably  different  at  the  thin- 
binder  end  of  the  burning  surface,  although  it  seems  likely  that  they  will 
not.  It  should  be  noted  that,  while  the  side-view  pictures  seem  to  show 
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a  line  of  vertical  flamelets,  end  view  pictures  reported  in  the  past  show 
the  flamelets  are  at  various  angles  to  the  mean  burning  surface,  giving 
the  "brushy"  appearance. 

Partial  results  from  the  photographic  test  series  in  Table  1  show 
that  a  transitory  flamelet  array  is  common  to  ail  the  binders  tested,  over 
all  the  pressures  that  yielded  a  view  of  the  relevant  region  of  the  combus¬ 
tion  zone  (in  the  tests  at  1.4  MPa,  200  psi,  the  view  was  obscured  by  the 
wall  of  unburned  AP;  in  4.1  MPa  tests  the  view  was  partially  obscured). 

The  sandwiches  with  PS  binder  gave  such  low  flame  luminosity  that  this 
behavior  could  not  be  determined,  or  included  in  the  description  of  flame- 
let  behavior. 

The  transitory  nature  of  the  diffusion  flamelets  led  to  the  specula¬ 
tion  that  a  continuous  flame  sheet  might  not  be  dynamically  stable,  a  spe¬ 
culation  that  motivated  the  analysis  in  Section  7.  Another  speculation 
about  the  transitory  flamelet  field  stems  from  the  fact  that  the  flamelets 
originate  in  the  region  where  the  AP  surface  processes  appear  to  be  influ¬ 
enced  by  the  binder  (region  of  the  "smooth  band").  The  nature  of  this 
surface  is  unknown,  but  it  is  reasonable  to  speculate  that  the  nature  is 
locally  intermittent,  (an  example  of  intermittent  behavior  would  be  emer¬ 
gence  of  oxidizer  vapors  through  a  thin  binder  melt  overlay  by  bubble  blow¬ 
ing)  . 

Whatever  the  cause  of  the  evident  instability  of  individual  flamelets, 
the  resulting  nature  of  the  flamelet  field  is  a  matter  of  considerable 
practical  importance  in  that  it  may  play  a  critical  role  in  the  response 
of  combustion  to  flow  disturbances  that  lead  to  combustion  instability. 

Of  course  understanding  of  the  cause  of  the  flamelet  response  will  suggest 
means  to  change  it. 

3.9  Discussion  of  Sandwich  Burning  Results 

3.9.1  Results  of  the  present  investigation  are  similar  to  those  of  earlier 
investigations  in  some  respects,  different  in  others.  The  differences 
are  most  conspicuous  in  the  thin-binder  geometry,  which  also  is  the  confi¬ 
guration  least  investigated  by  others  but  most  closely  simulating  propellant 


45 


combustion  conditions.  The  similarities  with  previous  studies  (i.e., 
arising  from  tests  of  similar  geometry  and  test  conditions)  include: 

1.  Consistency  with  the  general  classification  system  in  Figures 
6  and  7. 

2.  A  characteristic  (pressure  dependent)  pattern  of  depressions 
and  evidence  of  a  frothy  melt  on  the  AP  surface. 

3.  Presence  of  a  leading  "edge"  of  the  burning  surface  of  the  AP 
near,  but  not  at,  the  AP-binder  interface. 

4.  Anomalous  smooth  surface  quality  of  the  AP  between  the  leading 
"edge"  and  the  interface. 

5.  Width  of  the  band  of  modified  surface  in  4  comparable  to  the 
particle  diameter  of  AP  in  propellants. 

6.  Evidence  of  melt  flow  over  the  interface  (not  certain  at  low 
pressure  or  with  thin  binder). 

7.  A  tendency  for  protrusion  of  the  binder  lamina  to  disappear  as 
the  binder  thickness  is  reduced. 

8.  Diffusion  flames  that  are  "brushy"  when  viewed  in  line  with  the 
sandwich  lamina. 

Novel  results  in  the  present  investigations  include: 

1.  Minimal  protrusion  of  binder  with  thin  binder  lamina,  recessed 
in  low  pressure  tests  (possibly  all  pressures). 

2.  Persistence  of  the  anomalous  "smooth"  band  of  AP  surface  adjoin¬ 
ing  the  interface,  even  with  low  binder  thicknesses  not  previously 
studied. 

3.  Apparent  retardation  of  AP  regression  rate  in  the  "smooth  band", 
leading  to  the  presence  of  a  leading  edge  of  the  AP  front  outside 
the  smooth  band,  even  with  thin  binder  at  low  pressure,  where 
the  binder  is  recessed  (Figure  12,  14,  15,  17). 

4.  Retardation  of  the  smooth  band  so  extreme  at  high  pressure  that 
the  unaffected  AP  surface  apparently  regresses  laterally  under 
the  retarded  smooth-band  surface  (PBAN  binder,  Figure  16). 

5.  Substantial  difference  in  state  of  binder  surface  (e.g.,  fluid¬ 
ity,  residue)  as  a  function  of  kind  of  binder.  Specifically, 

HTPB  and  PBAN  binders  showed  the  most  evidence  of  molten  state, 
and  CTPB  binder  showed  a  relatively  less  viscous  melt.  HTPB 


and  CTPB  binders  tended  to  show  more  char  than  the  others.  PS 
binder  showed  the  least  evidence  of  melt  state  (these  judgements 
are  somewhat  subjective,  based  on  evidence  of  melt  flow  and/or 
bubbles  from  quenched  samples  and  observations  by  combustion 
photography) . 

6.  Unsteadiness  of  the  diffusion  flame,  that  leads  to  intermittent 
flamelets  with  size,  spacing  and  duration  being  larger  for  thick 
binder,  smaller  for  thin  binder. 

7.  A  low  pressure  limit  for  sustained  burning  of  thin-binder  sand¬ 
wiches  (1.4  MPa  for  30  p,m  binder). 

8.  Several  detailed  differences  in  combustion  behavior  according 
to  binders  used. 

9.  Non-reproducible  combustion  behavior  with  thin-binder  sandwiches 
that  was  apparently  due  to  locally  intermittent  behavior,  and 
sometimes  to  proximity  to  the  deflagration  limit. 

3.9.2  The  "smooth  banrf1  on  the  AP  surface  adjoining  the  interface  appears 
to  be  a  major  factor  in  behavior  on  the  dimensional  scale  of  propellant 
microstructure,  but  the  nature  or  cause  of  this  band  is  still  unknown. 

It  has  been  attributed  in  the  past  to  binder  melt  flow.  This  argument 
was  based  on  results  with  thick-binder  sandwiches,  where  the  flow  pattern 
over  the  interface  is  visible  (References  11,  12)  and  the  binder  film  is 
visible  in  profile  of  sectioned  samples  (Reference  6).  This  argument  is 
less  convincing  with  thin-binder  sandwiches  where  a  "deficiency"  of  binder 
leads  to  recessed  binder  surfaces,  but  does  not  eliminate  the  smooth  band. 
Since  the  leading  edge  in  the  AP  surface  is  at  the  outer  edge  of  the  smooth 
band,  the  diffusion  flame  is  exerting  some  effect  at  least  that  far  out. 
Thus  the  smooth  band  is  probably  on  the  fuel-rich  side  of  the  diffusion 
flame  mean  position.  However,  this  flame  may  actually  consist  of  an  array 
of  transitory  flamelets  issuing  from  the  smooth  band — possibly  at  the  sites 
of  more  or  less  randomly  bursting  bubbles  in  the  smooth  band.  Existing 
observational  methods  have  not  thus  far  provided  adequate  resolution  to 
determine  the  actual  details,  but  the  combustion  photography  establishes 
the  short  duration  of  the  flamelets  (of  order  10  to  10  sec). 


Figure  19.  Successive  burning  surfaces  reflecting  retardation  of  the 
AP  burning  surface  near  binder. 

(a)  Sandwich  (b)  Particle 

While  the  microscopic  and  chemical  details  of  the  smooth  band  and 
diffusion  flamelets  are  still  unresolved,  it  seems  clear  that  the  band 
is  a  region  of  retarded  regression  of  the  AP  (relative  to  the  flame-assis¬ 
ted  leading  edge).  In  the  high  pressure  tests  this  retardation  is  so  great 
that  the  unassisted  AP  deflagration  undercuts  the  surface  of  the  band  (PBAN 
binder).  These  effects  are  on  an  easily  resolvable  scale,  the  same  scale 
as  AP  particle  size  in  propellants.  Thus  the  results  raise  profound  ques¬ 
tions  about  how  AP  particles  burn  in  propellants,  and  what  constitute  valid 
models  for  propellant  combustion.  This  is  suggested  by  the  sketch  of  suc¬ 
cessive  burning  surface  profiles  in  Figure  19,  where  part  a  is  the  inter¬ 
pretation  of  results  of  13.8  MPa  sandwich  burning  tests,  and  part  b  is 
a  corresponding  speculation  about  particle  burning  in  propellants.  Not 
only  does  this  type  of  particle  burning  raise  difficulties  with  existing 
analytical  models  of  propellant  burning  rate,  but  adds  features  to  the 
combustion  zone  structure  that  are  likely  to  be  important  in  the  dynamic 
response  that  drives  combustion  instabilities. 

3.9.3  As  with  any  simplified  model  of  a  complex  system,  the  sandwich  mode l 
mus t  fail  to  simulate  some  important  aspects  of  propellant  combustion. 
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While  the  microscopic  and  chemical  details  of  the  smooth  band  and 
diffusion  flamelets  are  still  unresolved,  it  seems  clear  that  the  band 
is  a  region  of  retarded  regression  of  the  AP  (relative  to  the  flame-assis¬ 
ted  leading  edge).  In  the  high  pressure  tests  this  retardation  is  so  great 
that  the  unassisted  AP  deflagration  undercuts  the  surface  of  the  band  (PBAN 
binder).  These  effects  are  on  an  easily  resolvable  scale,  the  same  scale 
as  AP  particle  size  in  propellants.  Thus  the  results  raise  profound  ques¬ 
tions  about  how  AP  particles  burn  in  propellants,  and  what  constitute  valid 
models  for  propellant  combustion.  This  is  suggested  by  the  sketch  of  suc¬ 
cessive  burning  surface  profiles  in  Figure  19,  where  part  a  is  the  inter¬ 
pretation  of  results  of  13.8  MPa  sandwich  burning  tests,  and  part  b  is 
a  corresponding  speculation  about  particle  burning  in  propellants.  Not 
only  does  this  type  of  particle  burning  raise  difficulties  with  existing 
analytical  models  of  propellant  burning  rate,  but  adds  features  to  the 
combustion  zone  structure  that  are  likely  to  be  important  in  the  dynamic 
response  that  drives  combustion  instabilities. 

3.9.3  As  with  any  simplified  model  of  a  complex  system,  the  sandwich  model 
must  fail  to  simulate  some  important  aspects  of  propellant  combustion, 


and  may  also  produce  some  effects  that  are  of  only  secondary  importance 
with  propellants.  These  are  the  prices  that  must  be  paid  for  the  benefits 
of  controlled  variables,  and  simplicity  of  interpretation  of  test  results. 

Perhaps  of  first  importance  to  relevance  is  the  distance  of  things 
on  the  propellant  burning  surface.  Oxidizer  particles  in  propellant  ordi¬ 
narily  have  an  effective  radius  in  the  range  of  2  -  200  micrometers,  so 
that  distances  from  oxidizer  binder  interfaces  are  in  that  range,  and  the 
flame  environment  occurring  in  that  range  is  the  one  relevant  to  propellant 
surfaces.  A  relevant  theory  or  experiment  is  one  that  correctly  describes 
conditions  or  processes  in  the  200  ^m  range  from  the  interface.  This  im¬ 
plies  that  relevant  experimental  observations  of  sandwich  burning  mus t 
pertain  to  the  200  u,m  range  relative  to  the  interface,  and  the  experiment 
must  simulate  propellant  behavior  in  that  region. 

Past  work  pertaining  to  sandwich  burning,  both  analytical  and  experi¬ 
mental,  has  not  been  too  faithful  to  the  above  test  of  relevance  to  pro¬ 
pellant  burning,  a  point  that  will  not  be  belabored  here  except  to  point 
out  that  increasing  emphasis  has  been  placed  in  the  present  work  on  details 
close  to  the  binder  interface,  and  on  use  of  binder  thicknesses  comparable 
to  binder  structural  elements  in  the  matrix  of  propellants. 

A  second  basic  concern  regarding  relevance  of  sandwich  burning  to 
propellants  is  the  relatively  steady  state  nature  of  burning  in  sandwiches . 
Being  made  primarily  of  particulate  ingredients,  the  combustion  of  hetero¬ 
geneous  propellants  is  locally  intermittent  (on  a  time  scale  corresponding 
to  the  burning  time  of  the  ingredient  particles).  Thus  any  condition  on 
the  burning  surface  of  a  sandwich  that  requires  a  generation  time  longer 
than  the  propellant  particle  burning  time  will  not  be  fully  developed  on 
the  propellant  surface  because  of  the  interruption  at  particle  burnout. 

An  example  of  this  is  the  "needle"  structure  reported  (Reference  9)  on 
the  deflagration  surface  of  AP  at  high  pressure,  which  involves  a  layer 
of  thickness  about  the  same  as  the  diameter  of  AP  particles  in  propellants. 
This  does  not  mean  that  the  phenomena  observed  during  more  nearly  steady 
burning  are  not  relevant,  but  rather  that  those  resulting  from  slow  pro¬ 
cesses  will  not  become  fully  developed  in  particle  burning. 

A  third  basic  concern  regarding  relevance  of  sandwich  burning  has 
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to  do  with  the  " edgewise  burning"  of  the  interface.  In  the  sandwich,  the 
geometry  constrains  the  burning  front  to  proceed  in  the  direction  of  the 
interface.  Because  of  the  granular  nature  of  ingredients  in  propellants, 
the  interfaces  are  not  only  discontinuous,  but  inclined  in  varying  degree 
to  the  mean  burning  surface  regression.  There  may  be  a  significant  portion 
of  the  particle  burning  in  which  edge  burning  is  approximated,  but  it  is 
not  presently  established  what  contribution  each  part  of  the  particle 
burning  history  makes  to  the  overall  burning  rate,  or  to  the  structure 
of  the  burning  surface  and  gas  phase  combustion  zone.  These  questions 
cannot  be  answered  by  sandwich  burning,  but  are  probably  important  issues 
in  determining  steady  and  oscillatory  burning  characteristics  and  proper 
mechanistic  models  thereof. 

3.9.4  From  the  results  of  the  sandwich  burning  tests  and  the  foregoing 
considerations,  it  is  concluded  that  the  principal  results  of  value  to 
propellant  applications  come  from  examination  of  the  surface  behavior 
within  about  200  micrometers  of  the  binder-oxidizer  interface  and  the 
diffusion  flames  issuing  from  that  region.  Further,  it  appears  that  re¬ 
sults  relevant  to  propellants  can  be  obscured  if  thick  binder  laminae  are 
used.  In  addition  unless  caution  in  interpretation  is  exercised,  one  may 
draw  conclusions  from  results  of  sandwich  burning  experiments  that  are 
not  relevant  to  propellants  because  of  the  quasi-steady  nature  of  sandwich 
burning  (as  opposed  to  the  interraittency  induced  by  the  granular  nature 
of  the  propellant  microstructure). 

In  the  present  work,  the  "novel  results"  in  3.9.1  above  are  believed 
to  all  be  relevant  to  propellants,  and  pose  a  rather  different  picture 
of  propellant  combustion  zone  behavior  than  that  generally  described  in 
the  past.  This  is  illustrated  by  the  somewhat  speculative  picture  of  the 
combustion  zone  microstructure  in  Figure  20.  The  accompanying  photomi¬ 
crographs  of  quenched  burning  surfaces  attest  to  the  complexity  of  the 
real  combustion  zone  structure  on  a  microscopic  scale  and  illustrate  de¬ 
tails  similar  to  the  sketch,  as  modified  by  the  expected  trends  in  the 
features  of  the  combustion  zone  as  discussed  in  the  following. 
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The  sketch  in  Figure  20a  is  based  on  behavior  at  about  3.5  MPa  (  ~ 

500  psi),  corresponding  to  the  second  part  of  Figure  20b.  The  sketch  is 
also  based  on  an  AP  particle  diameter  of  about  300-400  |j,m.  The  features 
of  the  combustion  zone  would  be  expected  to  vary  with  particle  size,  pres¬ 
sure,  time  during  burning,  etc.,  in  the  following  manner. 

1.  SMOOTH  BAND,  present  at  all  pressures,  retards  surface  regres¬ 
sion.  May  cover  whole  particle  when  particles  are  small.  May 
have  progressive  effect  with  time  during  burning  of  each  par¬ 
ticle.  May  provide  anchor  sites  of  diffusion  flamelets. 

★ 

2.  LEADING  EDGE  OF  AP  surface,  caused  by  di f f us icjn- flame  augmented 
AP  deflagration.  Contribution  of  the  flame  to  regression  rage 
at  high  pressure  (e.g.,  10-15  MPa)  is  minor,  contribution  at 
low  pressure  (<2  MPa)  is  dominant.  Effect  may  be  washed  out 

by  smooth-band  effect  with  small  particles. 

3.  AP  DEFLAGRATION  dominates  AP  surface  in  the  central  area  of  the 
particle  surface,  especially  with  large  particles  and  at  higher 
pressures.  At  lower  pressures  this  area  is  raised  because  its 
regression  is  not  effectively  assisted  by  the  diffusion  flame 
(large  particles). 

4. -  DIFFUSION  FLAMELETS  occur  in  the  vicinity  of  the  oxidizer-binder 

interface,  and  are  intermittent  and  localized  on  a  time  scale 
of  about  l  millisecond.  The  causes  of  this  behavior  are  not 
determined  by  the  experiments.  However,  it  is  significant  that 
the  flamelets  apparently  originate  over  the  smooth  band  rather 
than  the  interface,  implying  that  their  singular  behavior  may 
be  related  to  the  nature  of  the  smooth  band.  Since  the  "root" 
of  the  diffusion  flame  dominates  the  overall  flame  behavior  and 
its  contribution  to  burning  rate  and  dynamic  behavior,  it  seems 
crucial  that  the  nature  and  role  of  the  smooth  band  be  further 
clari f ied. 


In  this  report  the  "diffusion  flame"  is  used  to  designate  the  AP 
binder  flame,  regardless  of  whether  its  rate  is  "diffusion  limited" 
or  "kinetically 
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TESTS  ON  OTHER  ORDERED  MICROSTRUCTURES 

4.1  Philosophy  of  Tests 

Use  of  simple  microstructures  such  as  sandwiches  is  helpful  in  un¬ 
derstanding  tests  results,  but  any  such  configuration  necessarily  fails 
to  exhibit  faithfully  all  aspects  of  propellant  combustion  behavior.  In 
the  following,  some  examples  of  this  are  described,  along  with  some  ideas 
about  ordered  microstructures  that  might  serve  well  to  demonstrate  aspects 
of  combustion  behavior  not  exhibited  by  sandwiches.  Following  this  dis¬ 
cussion,  some  exploratory  tests  will  be  described,  tests  designed  to  sup¬ 
plement  the  results  obtained  in  sandwich  burning  tests.  These  tests  are 
catalogued  in  Table  2. 

4.2  Combustion  Behavior  and  Test  Design 

In  combustion  of  real  propellants,  there  are  aspects  of  behavior  that 
depend  on  the  smallness  of  distances  in  the  propellant  microstructure, 
the  interaction  of  neighboring  flames,  and  intermittency  of  composition 


* 

Table  2.  Tests  on  Ordered  Structures 


Check  on  burning  rate  of  pressed 
AP  alone 

Beaded  Sandwiches 


Nylon  threads  in  AP 
HTPB  threads  in  AP 

Sandwich  with  600  y,m  spheres 

Stacked  Cylinders 


4.14  MPa 


Quench 


1/2 

0/3 


Movie 


6.89  MPa 


Quench  Movie 


8/9 

5/8 

4/4 

2/5 

l/l 

1/6 

3/3 

5/5 

3/3 

1/3 

5/6 

with 

useful 

tun  - - - - -  —  , 

test  results  (i.e.,  apparatus  function  OK  and  sample  burn  regular) . 
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Figure  21.  Nature  of  flame  sheets. 

(a)  Closed  interface  line  typical  of  oxidizer  particles  in 
propellants 

(b)  Open  interface  line  typical  of  sandwiches 

(c)  Experimental  arrangement  for  studying  behavior  of  a  closed- 
interface  counterpart  of  a  sandwich  (fuel  coated  oxidizer 
rod  in  oxidizer) 

in  any  line  of  progression  of  the  burning  surface  (due  to  heterogeneity). 
These  are  typical  propellant-dependent  differences  from  the  edge-burning 
sandwich  model,  and  these  differences  are  known  to  be  important  to  steady 
state  and/or  oscillatory  combustion  of  propellants. 

The  foregoing  is  illustrated  by  the  following  example.  One  might 
reasonably  question  to  what  extent  the  surface  and  flame  behavior  observed 
in  sandwich  burning  would  occur  if  the  diffusion  flame  and  interface-re¬ 
lated  effects  issued  from  a  small  closed  interface  curve  such  as  the  border 
of  a  partly  burned  AP  particle  in  a  propellant  burning  surface,  instead 
of  the  linear  interface  of  a  sandwich  (Figure  21).  At  what  diameter  of 
a  closed  interface  do  the  adjoining  and  opposite  elements  of  the  flame 
surface  (flamelets)  or  the  "smooth  band"  start  interacting  with  each  other, 
and  what  is  the  effect  on  propellant  burning  rate  or  combustion  stability? 
Presumably  these  issues  would  be  clarified  by  studying  combustion  of  rods 
of  AP  coated  with  binder  and  surrounded  by  more  AP  (Figure  21c).  End¬ 
burning  behavior  as  a  function  of  rod  diameter  would  give  some  clues  to 
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Figure  22.  Arrangements  of  fuel  filaments  in  oxidizer. 

(a)  Two  filaments,  far  enough  apart  for  separate  flames 

(b)  Two  filaments,  interacting  flames 

(c)  Two  filaments,  different  size,  interacting  flames 

the  real  propellant  effects  posed  by  the  initial  question. 

A  second  fundamental  question  of  composite  propellant  burning  is  the 
interaction  of  dif fuion  f lamelets  from  adjoining  oxidizer  particles  (Figure 
22).  Do  flamelets  from  neighboring  particles  tend  to  support  each  other 
(possibly  suppressing  flickering),  or  "repel"  each  other  (by  consuming 
each  other's  local  fuel  supply)?  Does  the  behavior  of  a  big  flame  envelope 
of  a  large  oxidizer  particle  modify  the  smaller  envelopes  of  smaller  neigh¬ 
boring  oxidizer  particles  (or  vice  versa)?  If  adjoining  envelopes  flicker, 
do  they  flicker  together?  Do  such  interaction  effects  possibly  develop 
only  in  oscillatory  environments,  and  if  so  at  what  frequency?  Some  of 
these  questions  reflect  mechanistic  issues  currently  being  considered  in 
computer  modeling  of  the  combustion,  but  the  factual  information  on  real 
behavior  is  lacking.  The  real  combustion  behavior  could  be  studied  in 
a  controlled  laboratory  experiment  in  which  the  oxidizer  particles  are 
replaced  by  oxidizer  rods,  packed  in  varying  degrees  of  proximity.  How¬ 
ever,  such  a  configuration  is  difficult  to  produce,  and  in  the  present 


work,  the  experiment  is  " inverted'  by  using  small  binder  rods  pressed 
in  an  AP  test  sample. 

The  third  fundamental  aspect  of  composite  propellant  combustion  that 
is  not  addressed  by  sandwich  models  (or  analytical  models)  is  the  local 
intermittency  of  the  microstructure,  i.e.,  the  structural  and  chemical 
discontinuities  corresponding  to  particulate  ingredients.  In  an  ideal 
experiment,  it  would  be  desirable  to  start  with  intermittency  while  re¬ 
taining  an  ordered  microstructure,  thus  avoiding  those  complications  in 
characterization  of  combustion  that  are  due  to  the  difficulty  of  charac¬ 
terizing  the  original  sample  microstructure  in  the  first  place.  Some 
efforts  were  made  in  the  present  program  to  achieve  this  goal,  without 
much  success  to  date. 

4.3  Fuel  Strands  -  Beaded  Sandwiches 

In  an  effort  to  produce  diffusion  flame  arrays  that  are  more  localized 
than  conventional  sandwiches  provide,  tests  were  run  on  the  configuration 
shown  in  Figure  23b.  The  configuration  in  part  b  of  the  figure  was  pro¬ 
duced  by  modification  of  the  sandwich  fabrication.  The  die  used  in  press¬ 
ing  the  AP  sheets  was  milled  so  that  its  face  had  a  series  of  small  par¬ 
allel  ridges  (Figure  23a),  which  produced  small  parallel  grooves  in  the 
pressed  surface.  Sandwiches  were  made  by  pairing  these  special  sheets 
with  flat  ones,  and  using  minimal  HTPB  binder  in  the  flat  interface  (20-25 
p,m) .  Excess  binder  filled  the  grooves,  giving  binder  "beads"  or  "fila¬ 
ments"  roughly  160  ym  square  and  80  ym  square  (two  different  pressing  dies). 
Tests  were  made  on  those  samples  by  edge  burning,  with  combustion  photo¬ 
graphy  and  rapid  depressurization  quench  (Table  2).  Tests  on  the  "Beaded 
Sandwiches"  were  not  particularly  rewarding,  but  results  are  worth  record¬ 
ing.  Combustion  photography  showed  a  luminous  flame  from  each  binder 
"bea<f' ,  very  little  flame  elsewhere.  Figure  24  shows  three  typical  se¬ 
quences  of  frames  from  motion  pictures  taken  at  1000  frames  per  second 
(the  first  sequence  is  taken  without  external  illumination).  The  pictures 
show  an  irregular  flame  envelope  above  a  160  u>m  binder  bead,  which  changes 
only  slowly  with  time  compared  to  the  flamelets  seen  in  flat  sandwiches. 
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Figure  23.  Arrangement  for  sandwiches  with  beaded  binder  lamina 


Figure  24.  Picture  sequences  showing  flame  at  the  site  of  the  beads 
(ridges)  on  the  binder  lamina. 

(a)  No  external  illumination 

(b)  and  (c)  Sequences  with  external  illumination  (in  (b)  and 
(c)  the  large  light  area  below  the  flames  is  reflected 
external  light)  Pictures  taken  at  1000  frames  per  second 


This  result  suggests  that  the  closed  flame  envelope  may  be  more  stable 
(or  unstable  at  lower  frequency)  than  a  plane  flame  sheet.  This  is  a 
speculative  conclusion,  because  the  geometry  of  the  burning  surface  was 
too  complex  for  interpretation  of  behavior.  Further,  there  is  usually 
a  projecting  char  of  binder  inside  the  flame  envelope  and  even  with  such 
a  "flame  holder",  the  flame  envelope  was  not  stationary,  but  simply  fluc¬ 
tuated  in  a  slower  manner  than  the  plane  flame  sheets  of  the  conventional 
sandwiches . 

The  general  surface  character  of  the  quenched  beaded  sandwiches  was 
similar  to  other  thin-binder  sandwiches  burned  at  the  same  pressure  (6.9 
MPa,  1000  psi),  except  near  the  sites  of  the  binder  beads.  The  binder 
was  roughly  50  Hm  thick  at  sites  other  than  the  beads.  The  AP  surface 
sloped  upward  towards  the  AP-binder  interface  in  a  manner  similar  to  re¬ 
gular  sandwiches,  and  the  binder  laminae  ranged  from  slightly  recessed 
to  slightly  projecting  relative  to  the  interface  oxidizer  (Figure  25). 

The  smaller  beads  of  binder  (  ~  125  nm)  generally  projected  above  the 
adjoining  AP,  in  a  manner  corresponding  to  behavior  of  conventional  sand¬ 
wiches  with  125  urn  binder  (Figure  25a) .  The  samples  with  250  p.m  binder 
beads  yielded  recesses  at  the  sites  of  the  beads  (Figure  25b),  with  evi¬ 
dence  of  melt  behavior  in  the  recesses.  This  result  is  not  consistent 
with  conventional  sandwich  behavior,  but  occurred  consistently.  There 
may  be  a  problem  of  injection  of  molten  binder  during  quench,  a  classic 
but  never  verified  rationalization  of  unexplained  features  of  dp/dt  quench¬ 
ed  samples  (this  point  was  not  raised  earlier  in  this  report  because  re¬ 
cessed  binders  were  noted  primarily  under  conditions  not  conducive  to 
injection  of  melt  during  quench,  e.g.,  thin  laminae,  low  pressures,  and 
particularly  with  PS  binder  which  shows  minimal  signs  of  melts). 

Binder  melt  flow  over  the  AP  surface  and  "smooth  band"  effects  were 
present  but  not  conspicuous  with  the  beaded  sandwiches  at  the  pressures 
tested.  Samples  were  burned  upside  down  and  quenched,  and  gave  nominally 
similar  quenched  surfaces,  with  some  enhancement  of  the  binder  protrusion 
at  the  beads. 


(Continued)  Surface  of  quenched  beaded-binder  sandwiches  (6.9  MPa),  (b)  Large  beads 
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Figure  26.  Arrangement  for  hydrostatic  pressing  AP  samples  with  fuel 
filaments . 

4.4  Fuel  Strands  -  Fuel  Filaments 

Another  approach  to  fabrication  of  fuel  strands  in  AP  was  to  dry  press 
powdered  AP  samples  containing  threads  of  polymer.  For  this  purpose,  samples 
of  monofilament  nylon  fish  line  were  used.  The  segments  of  thread  were 
simply  laid  in  the  AP  powder  in  the  pressing  die,  covered  by  more  AP,  and 
pressed  into  wavers  in  the  usual  manner  in  the  arbor  press.  In  an  alternate 
approach,  samples  were  pressed  in  a  hydrostatic  press.  These  samples  were 
prepared  by  filling  a  rubber  tube  (Figure  26)  with  AP  powder,  using  rubber 
stoppers  to  seal  the  samples  in  the  tube.  Monofilament  threads  were  sup¬ 
ported  through  the  sample  by  threading  them  through  the  stopper.  This 
sample  assemblage  was  subjected  to  207  MPa  hydrostatic  pressure  for  60 
minutes.  The  pressed  sample  appeared  to  be  well  consolidated,  but  lacked 
the  nearly  translucent  appearance  of  the  die-pressed  wafers.  In  a  further 
adaptation  of  the  method,  samples  were  made  by  both  the  hydrostatic  press 
and  die-press  methods  with  the  polymer  filament  replaced  by  fine  steel 
wires.  After  pressing,  the  steel  wires  were  withdrawn  and  the  holes  were 
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filled  with  uncured  propellant  binder.  This  approach  had  the  advantage 
of  assuring  a  uniform  polymer  filament  (not  deformed  by  AP  particle  in¬ 
trusion  during  pressing)  and  permitting  use  of  conventional  propellant 
binders.  As  noted  below,  test  results  on  samples  made  by  all  methods 
suggested  that  voids  and/or  a  gap  between  the  binder  and  AP  in  the  hole 
may  have  existed  in  many  samples  tested.  However,  further  work  seems  to 
hold  promise. 

Actual  tests  on  the  fuel  strand  are  listed  in  Table  2.  Hydrostati¬ 
cally  pressed  samples  tended  to  burn  rapidly  down  the  interface  between 
the  fuel  strand  and  the  AP.  This  was  true  of  both  HTPB-filled  types  and 
pressed- in-place  nylon  filament  samples.  Burning  rate  tests  on  samples 
prepared  without  the  fuel  strands  were  found  to  be  much  higher  than  with 
die-pressed  samples  of  AP.  These  "suspicious"  results  led  to  abandonment 
of  the  hydrostatic  pressing  method.  The  behavior  was  not  satisfactorily 
explained,  but  poor  AP  consolidation  was  suspected 

Fuel  strand  samples  made  by  the  more  familiar  die-pressing  method 
also  led  to  more  rapid  burning  than  expected  down  the  interface,  but  the 
burning  was  even,  and  yielded  quenched  samples  with  axially  symmetrical 
profiles  about  the  axis  of  the  fuel  strand.  An  example  of  a  quenched  two- 
strand  HTPB-filled  sample  quenched  from  6.9  MPa  (1000  psi)  is  shown  in 
Figure  27.  The  high  magnification  part  of  the  figure  (part  b)  shows  that 
the  sample  burned  with  the  now- familiar  "smooth  ban<f'  in  the  AP  near  the 
binder  interface.  However  combustion  has  proceeded  much  further  down  the 
AP  surface  there  than  in  the  rest  of  the  AP.  Similar  behavior  was  observed 
at  4.15  MPa  (600  psi),  as  shown  in  Figure  28.  With  nylon  strands,  the 
region  of  the  fuel-oxidizer  interface  was  similar  to  that  with  HTPB,  and 
burning  proceeded  far  ahead  of  the  main  AP  surface  as  with  HTPB  strands. 
Results  with  nylon  were  inconclusive  because  the  one  die-pressed  sample 
that  was  quench-tested  was  quenched  so  early  in  burning  (Figure  29)  that 
penetration  of  burning  down  the  interface  region  may  not  have  had  suffi¬ 
cient  time  to  yield  a  fully  developed  profile. 

Although  the  number  of  tests  on  strands  and  beaded  sandwiches  was 
limited,  the  results  pose  a  striking  contrast  (Figure  30).  The  beaded 
sandwiches  yielded  retarded  regression  of  the  AP  in  the  region  near  the 
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Figure  29.  Surface  of  quenched  AP  sample  having  nylon  filament 
die-pressed  in  the  sample.  Test  at  6.9  MPa. 
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Figure  30.  Comparison  of  quench  profiles  for  HTPB  strands  in  AP  and  HTPB 
beads  (ridges)  on  sandwich  laminae. 

(a)  Profile  of  cast  HTPB  strand  sample 

(b)  Profile  of  beaded  HTPB  sandwich  (section  through  bead, 
perpendicular  to  laminae),  250  p.m  bead 

(c)  Profile  of  beaded  HTPB  sandwich,  125  lim  bead 

fuel  beads,  while  the  samples  with  pre-forraed  threads  and  with  filled  holes 
showed  accelerated  AP  regression  in  the  same  region.  However  both  showed 
the  typical  smooth  band  with  retardation  of  rate  immediately  adjacent  to 
the  fuel.  Both  results  are  suspect  because  experience  with  sample  prepara¬ 
tion  is  limited  and  quality  is  suspect. 


4.5  Ordered  Arrays  of  Spherical  Oxidizer  Particles 

In  order  to  study  those  aspects  of  propellant  combustion  that  relate 
to  the  heterogeneity  ui  the  direction  of  surface  regression,  a  microstruc¬ 
ture  consisting  of  ordered  arrays  of  oxidizer  particles  is  needed.  Am¬ 
monium  perchlorate  spheres  of  uniform  size  would  be  desirable,  because 
they  can  be  arranged  in  a  "lattice-like"  array  to  give  a  determinate  micro¬ 
structure.  However  the  nearest  approximation  to  spheres  that  could  be 
obtained  was  a  commercial  grade  consisting  of  rounded  particles  of  nominal 
600  y,m  size.  These  were  not  sufficiently  uniform  in  shape  and  size  to 


b  7 


•00  ftm  AP  9PHCRU 


Figure  31.  Idealized  sandwich  with  ordered  array  of  oxidizer  spheres  in 
the  binder  lamina. 

yield  a  lattice-array  packing.  It  was  decided  to  test  them  in  a  conventional 
sandwich  configuration  with  the  particles  in  a  binder  layer  of  appropriate 
thickness  to  accomodate  one  layer  of  the  AP  particles  (Figure  31)  Four 
of  these  samples  were  prepared  and  tested,  with  inconclusive  results. 

Motion  pictures  were  obtained  on  one  test,  which  showed  the  burning  of 
the  filled  binder  layer  proceeding  rapidly  and  unevenly  ahead  of  the  AP 
sheets.  One  quenched  sample  was  obtained,  which  showed  irregular  burning 
had  occurred  with  local  penetration  to  the  interior  of  the  sample.  These 
difficulties  can  no  doubt  be  overcome,  and  further  use  of  this  configura¬ 
tion  is  anticipated. 

4.6  Stacked  Cylinder  Configuration 

Propagation  of  the  burning  surface  through  a  heterogeneous  propellant 
involves  not  only  deflagration  of  oxidizer  surfaces  and  propagation  down 
oxidizer-binder  interfaces,  but  also  burning  through  binder  layers  such 
as  those  separating  oxidizer  particles.  In  Section  6,  this  latter  aspect 


of  interparticle  burn- through  is  addressed  as  an  integral  part  of  a  burning 
rate  model.  However  it  is  recognized  that  very  little  is  known  about  this 
interparticle  burn-through  process. 

The  process  of  interparticle  burn-through  was  explored  experimentally 
by  preparing  and  burning  samples  (Table  2)  of  the  form  illustrated  in  part 
a  of  Figure  32.  These  samples  were  made  by  pressing  cylindrical  pellets 
of  AP,  6.35  mm  in  diameter,  and  cementing  them  in  the  indicated  array  using 
HTPB  binder.  The  samples  were  ignited  at  the  top,  and  motion  pictures 
were  taken  to  observe  how  the  deflagration  proceeded  through  the  sample. 
Figure  33  shows  three  frames  from  one  of  the  motion  pictures.  Part  b  of 
Figure  32  shows  the  sequences  of  burn-through  between  cylinders  from  five 
tests.  The  numbers  indicate  the  time  sequence  of  burn- through.  In  4  out 
of  5  cases,  the  bottom  cylinder  was  ignited  by  a  sequence  of  "diagonal" 
burn-throughs  before  a  direct  vertical  burn-through  occurred. 

4.7  Discussion 

The  investigations  in  this  section  were  largely  exploratory,  and  the 
results  should  be  regarded  accordingly.  Of  the  methods  used,  some  appeared 
to  be  amenable  to  control,  some  were  not.  Fuel  strands  by  the  beaded  sand¬ 
wich  method  resulted  in  regular  burning,  in  a  manner  compatible  with  sand¬ 
wich  burning  results.  However  the  presence  of  a  char  on  the  binder  strand 
invalidates  the  method  as  a  test  of  flamelet  stability.  Further,  it  would 
be  desirable  to  establish  whether  binder  was  ejected  during  quench  of  the 
sandwich  with  coarse  binder  beads.  The  objective  of  current  testing  with 
this  type  of  sample  is  to  determine  the  trend  of  behavior  with  pressure. 

The  results  with  nylon  strands  pressed  in  place  suggest  that  dimen¬ 
sional  changes  after  pressing  are  leaving  a  gap  between  oxidizer  and  fuel, 
and  the  method  has  been  dropped.  Results  with  samples  prepared  by  filling 
holes  in  the  oxidizer  with  binder  after  pressing  are  inconclusive,  and 
the  samples  are  difficult  to  prepare  reproducibly .  However  the  method 
holds  promise. 

The  sandwiches  with  spherical  oxidizer  particles  in  the  binder  layer 
did  not  yield  any  reproducible  results,  but  it  is  anticipated  that  useful 


GMTOR  PASTE 


•)  BEFORE  IGNITION 


»)  RESULTS  OF  FIVE  TESTS 


Figure  32.  Stacked  array  of  AP  cylinders. 

(a)  Initial  configuration  with  binder  interfaces 

(b)  Diagram  showing  order  of  burn-through  between  cylinders 
for  five  tests 
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Figure  33.  Photographic  sequence  of  burn-through  between  stacked  cylinders 


Figure  34.  Sketch  of  propagation  of  burning  surface  through  propellant 
microstructure.  Example  argues  the  least  time  path  via  par¬ 
ticles  1  to  2  to  3,  interpreting  results  of  Figures  32  and 
33. 

At  time  t^,  interface  between  particle  1  and  2  is  exposed 
to  flame  of  particle  1,  pyrolyzes  binder  and  ignites  particle 
2  at  nearest  point. 

At  time  t  particle  1  is  burning  out  without  igniting 
particle  3.  Burning  front  in  particle  2  has  proceeded  from 
the  ignition  point. 

At  time  t  ,  the  deflagrators  front  in  particle  2  has  pro¬ 
gressed  past  the  nearest  point  to  particle  3  and  has  burned 
through  the  binder  to  ignite  particles  3  at  that  point.  In 
the  meantime,  particle  3  may  have  ignited  at  the  "top”,  but 
the  results  indicate  that  the  delay  causes  the  deflagration 
front  from  that  point  to  lag  the  one  originating  at  the  particle 
2-particle  3  point. 

results  concerning  AP  particle  burning  can  be  obtained  by  this  method, 
and  further  tests  are  planned. 

The  tests  with  stacked  oxidizer  cylinders  appeared  to  answer  a  fairly 
critical  question  concerning  the  least  time  path  for  propagation  of  the 
combustion  front  through  the  propellant  microstructure.  Specifically, 
burn-through  of  the  binder  layer  between  oxidizer  particles  appears  to 
be  a  significant  factor  in  propagation  of  the  combustion  front.  Further, 
burn-through  is  apparently  more  rapid  between  particles  adjoining  each 
other  at  lateral  surfaces  than  adjoining  particles  oriented  in  the  direction 
of  propagation  of  the  macroscopic  burning  surface  (Figure  34).  This  result 


suggests  that  burn- through  on  lateral  interfaces  is  supported  by  the  sus¬ 
tained  convection  of  hot  oxidizer,  which  is  absent  when  the  interfacial 
binder  is  reached  just  as  the  oxidizer  particle  burns  out.  However  this 
result  is  speculative  because  of  the  relatively  large  size  of  the  oxidizer 
cylinders  in  the  tests  (compared  to  oxidizer  particles  in  propellants). 

Collectively,  the  various  tests  on  ordered  structures  were  not  par¬ 
ticularly  rewarding,  and  relatively  difficult  to  carry  out.  The  results 
may  be  more  revealing  than  they  seem  to  be,  but  more  testing  is  needed 
to  distinguish  between  orderly  combustion  behavior  and  artifacts  of  the 
experimental  methods. 


SURFACE  CHARACTERISTICS  OF  QUENCHED  PROPELLANTS 


5.1  Philosophy  of  Tests 

In  keeping  with  the  philosophy  of  progressing  from  geometrically 
simple  systems  like  the  sandwich  samples  (Section  3)  to  the  laterally 
intermittent  structures  (sought  in  Section  4),  one  seeks  to  progress  to 
the  isotropic,  "semirandom"  heterogeneous  situation  in  real  propellants. 
Intermittency  in  the  burning  direction  was  discussed  in  Section  4.5  and 
4.6,  and  efforts  to  achieve  such  intermittency  using  ordered  structures 
were  described.  Results  were  only  marginally  successful,  and  it  was  de¬ 
cided  to  do  some  studies  using  a  set  of  available  propellants.  These  tests 
were  run  primarily  as  a  reference  point  for  testing  relevance  of  work  on 
model  systems  such  as  sandwiches.  Further  tests  are  planned  using  for¬ 
mulations  and  ingredients  tailored  to  evaluate  specific  aspects  of  com¬ 
bustion  behavior. 

5.2  Test  Design 

These  tests  were  run  on  a  set  of  three  propellants  provided  by  the 
Thiokol  Chemical  Corporation  (Huntsville)  in  connection  with  another  program. 
Their  primary  attributes  were  relatively  narrow  oxidizer  size  distributions 
and  a  contemporary  binder  (HTPB).  The  stated  average  particle  sizes  were 
48,  81  and  273  nm  for  the  three  samples.  The  nominal  formulations  were 
83  percent  AP  and  17  percent  binder,  but  were  varied  because  of  difficul¬ 
ties  in  processing.  The  details  are  shown  in  Table  3,  along  with  a  cata¬ 
loging  of  tests  that  were  run.  Also  shown  in  the  table  are  tests  run  on 


Table  3.  Summary  of  Propellants  Tested 


Propellant 

Designation 

Composition 

AP/Binder/AL 

AP  Mass  Avg. 
Particle  Size 
(micro  meters) 

Binder 

Type 

T-48 

70.0/30.0/0 

48 

HTPB 

T-81 

83.1/16.9/0 

81 

HTPB 

T-273 

85.0/15.0/0 

273 

HTPB 

GIT 

75.0/25.0/0 

600 

HTPB 

UTP3001* 

73.0/12.0/15.0 

Blend* 

PBAN 

♦Nominal  formation,  a  commercial  propellant,  see  JANNAF  Propellants  Manual. 
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a  propellant  made  at  GIT  using  the  rounded  600  Wn  oxidizer  referred  to 
in  Section  4.5.  Both  quench  tests  (dp/dt),  and  combustion  photography 
tests  were  carried  out,  with  the  photographic  procedures  designed  primarily 
to  observe  flame  behavior.  Quenched  samples  were  examined  by  optical  micro¬ 
scopy,  with  selected  samples  examined  by  scanning  electron  microscope. 

5.3  Results 

The  propellants  burned  in  a  normal  manner,  with  visible  orange  flames, 
with  some  smoke,  particularly  at  high  pressures  with  coarse  oxidizer. 

In  the  2-7  MPa  range  the  luminous  region  of  the  flame  was  2-5  mm  thick. 

For  the  propellants  with  finer  AP  particle  size ,  the  flame  layer  was  made 
up  of  a  large  number  of  plumes  which  extended  outward  from  the  surface 
in  straight  parallel  lines  (Figure  35).  Within  the  capability  of  the 
photography  to  resolve  it,  the  plumes  appear  to  be  transitory  on  a  time 
scale  of  one  or  two  milliseconds,  in  a  manner  similar  to  the  flickering 
observed  in  the  sandwich  burning  tests.  The  lateral  dimensions  and  spacing 
of  the  plumes  are  of  the  same  order  as  the  oxidizer  particle  size  (finer 
detail  would  not  be  resolved  in  the  pictures).  For  the  propellants  with 
coarser  AP  particle  size,  the  thickness  of  the  luminous  region  was  the 
same  (roughly)  as  with  fine  AP,  but  the  filaments  were  larger  and  more 
persistent,  with  the  greatest  width  near  the  base  being  of  the  same  order 
as  the  oxidizer  particle  size  (Figure  35).  It  is  not  valid  to  assume  that 
the  luminous  "flamelets"  are  a  faithful  representation  of  the  distribution 
of  reaction,  because  they  probably  reflect  the  distribution  of  moving  hot 
carbon. 

Figure  36  shows  the  surfaces  quenched  from  two  different  pressures. 

Parts  36a  and  b,  for  T-81  propellant, show  irregularly  shaped  oxidizer  par¬ 
ticles  with  irregular  quality  of  the  central  areas  of  the  particle  surfaces. 
Although  the  pictures  do  not  reveal  the  "height"  of  the  surface  details, 
visual  inspection  showed  the  oxidizer  surfaces  to  be  substantially  raised 
in  2.1  MPa  tests,  and  recessed  in  6.9  MPa  tests.  In  both  cases,  the  irregular 
porous  and  "bubble"  quality  of  the  central  areas  persisted,  and  these  cen¬ 
tral  areas  were  raised  relative  to  the  adjoining  "leading  edge"of  the 
AP  particle  surface  at  both  pressures.  Results  with  T-48  propellant  (not 
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Figure  35.  Frames  from  motion  pictures  of  combustion  of  AP-HTPB  propellants. 

Pictures  at  1000  frames  per  second;  frame  height  corresponds  to  7  mm 
a)  T-81  Propellant  at  2.1  MPa;  b)  T-273  Propellant  at  0.7  MPa 
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Quenched  burning  surfaces  of  AP-HTPB  propellants  quenched  by 
rapid  depression.  Pictures  by  scanning  electron  microscopy. 
T-81,  quenched  from  a)  2.1  MPa;  b)  6.9  MPa 
T-273,  quenched  from  c)  2.1  MPa;  d)  6.9  MPa 


shown  in  the  Figure)  were  similar  to  those  with  T-81  propellant,  except 
for  a  tendency  to  surface  flooding  with  molten  binder  or  its  decomposition 
intermediates . 

The  quench  tests  with  the  coarse-oxidizer  propellant  (T-273,  Figure 
36c,  d)  showed  the  oxidizer  particles  to  be  round,  and  the  quenched  par¬ 
ticles  exhibited  porosity  in  the  central  area  of  their  burned  surfaces, 
as  with  the  finer  oxidizer.  This  central  area  was  less  raised  in  the  6.9 
MPa  tests.  The  samples  showed  what  appeared  to  be  a  profusion  of  binder, 
suggestive  of  accumulation  and  melt  flow,  although  there  was  no  unequivocal 
test  used  to  map  the  binder-covered  areas  of  the  surface. 

In  general,  the  surface  profiles  were  consistent  with  results  of 
sandwich  burning  experiments.  Thus  at  all  pressures,  the  leading  "edge" 
of  the  regression  in  the  oxidizer  is  at  some  distance  (10  -  250  u,m)  in 
from  the  edge  of  the  particle  burning  surface  (most  evident  in  Figure  36c, 
d),  with  the  AP  surface  then  rising  again  further  from  the  interface  at 
lower  pressure  (36c),  but  simply  leveling  off  in  6.9  MPa  tests  (36d). 

In  the  coarse-oxidizer  propellant,  the  binder  accumulation  and  apparent 
flow  is  analogous  to  thick-binder  sandwich  combustion.  The  most  distinc¬ 
tive  feature  of  the  quenched  particles  (compared  to  sandwiches)  is  the 
relatively  greater  porosity  of  the  AP  surface  (far  from  the  interface) 
in  the  case  of  particles.  The  "smooth  band"  observed  in  the  sandwich¬ 
burning  tests  was  not  conspicuous  in  most  of  these  scanning  electron  mi¬ 
croscope  pictures  of  quenched  propellants.  This  is  apparently  an  artifact 
of  microscope  technique,  as  illustrated  by  Figure  37,  a  later  picture  of 
the  T-81  sample  in  Figure  36b.  In  this  picture,  the  smooth  band  is  readily 
observable  (the  conductive  film  used  on  the  surface  to  make  SEM  possible 
has  developed  flows,  especially  on  the  lower  left  smooth  band). 

Figure  38  shows  scanning  electron  microscope  pictures  of  the  quenched 
surface  of  the  GIT  propellant  (Table  3)  with  600  ^m  round  AP  particles. 

The  appearance  of  binder  melt  flow  is  particularly  convincing  in  this  propel] an: 
(25%  HTPB  binder).  Also  evident  are  the  porous  regions  in  the  central 
area  of  the  quenched  particle  surfaces.  The  rectangular  array  of  cracks 
on  some  particle  surfaces  have  been  observed  previously  on  the  quenched 
surfaces  of  single  AP  crystals  (Reference  8),  and  are  probably  produced 
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Figure  37.  Illustration  of  the  "smooth  band"  on  a  quenched  AP  particle. 

T-81  propellant  quenched  from  6.9  MPa.  Microscope  technique  for 
this  picture  shows  the  smooth  band  better  than  in  previous 
quenched  propellant  figures. 
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Figure  38.  Quenched  burning  surfaces  of  AP-HTPB  propellant  with  600pm 
AP,  25%  binder.  Samples  quenched  from  6.9  MPa. 


in  the  surface  layer  due  to  contraction  upon  cooling  back  through  the  cubic- 
orthorhombic  phase  transition  after  quench.  Due  to  the  pervasiveness  of 
the  binder,  no  generalizations  about  the  surface  profiles  near  the  peri¬ 
phery  of  the  AP  surface  are  ventured  here.  However,  it  is  notable  that 
the  nature  and  exposure  of  the  AP  surfaces  is  extremely  varied.  This  appears 
to  be  due  to  the  different  extent  of  burning  of  the  various  particles, 
and  suggests  a  complex  and  time-dependent  role  of  the  binder  during  the 
history  of  burning  of  each  particle. 

5.4  Discussion 

The  results  of  this  series  of  tests  are  similar  to  results  reported 
by  other  investigators,  and  are  useful  here  primarily  in  the  context  of 
the  results  of  tests  on  ordered  structures.  As  examples,  the  following 
analogies  and  contrasts  help  to  determine  what  aspects  of  tests  on  ordered 
structures  relate  to  conventional  propellants,  and  what  aspects  do  not: 

1.  The  "smooth  band",  with  apparent  retardation  of  the  AP  burning 
rate  is  common  to  both  propellants  and  ordered  configurations.  Whether 
the  smooth  surface  and  retarded  rate  are  causally  connected  is  not  estab¬ 
lished.  Further,  it  is  not  established  whether  the  burning  rate  is  ac¬ 
tually  retarded,  or  whether  it  simply  adjoins  a  region  of  the  AP  surface 
that  has  abnormally  high  rate  because  of  proximity  to  the  diffusion  flame. 

2.  Excess  binder  leads  to  accumulation  and  flow  of  binder  melts, 
which  can  dominate  the  propellant  surface.  In  sandwiches,  this  corresponds 
to  thick  binder  laminae,  which  result  in  melt  flow  over  the  oxidizer. 

The  "smooth  band"  was  initially  thought  to  be  a  result  of  binder  melt  flow 
over  the  AP  surface,  but  the  band  has  been  evident  in  many  tests  under 
conditions  where  such  flow  is  improbable,  e.g.,  low  pressure-thin  binder 
sandwiches.  Quench  tests  of  sandwiches  burned  in  an  upside-down  position 
have  not  shown  conspicuous  differences  in  either  melt  flow  or  smooth  bands, 
indicating  that  intermolecular  forces  are  dominant,  i.e.,  the  behavior 
is  not  fluid  flow  except  possibly  in  the  case  of  excess  binder  situations. 

3.  There  is  persistent  evidence  that  the  deflagration  of  AP  oxidizer 
involves  a  very  complex  surface  degradation.  This  is  even  more  conspicuous 
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in  propellants  than  with  sandwiches,  as  evidenced  by  the  porous  raised 
areas  near  the  center  of  the  burning  surfaces  of  particles  in  quenched 
propellants . 

4.  Given  the  complex  interaction  of  binder  and  oxidizer,  it  seems 
likely  that  the  role  of  the  binder  changes  with  time  during  the  AP  particle 
burning,  and  hence  that  the  AP  burning  itself  changes  with  time  during 
particle  burning.  This  is  an  unwelcome  development  from  the  standpoint 
of  analytical  modeling  of  propellant  combustion,  and  merits  further  exam¬ 
ination  to  estimate  the  range  of  behavior  involved  and  the  effect  on  burn¬ 
ing  rate  and  dynamic  response  of  the  combustion  zone. 
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6.0  LEAST  TIME  PATH  MODEL  FOR  BURNING  RATE 

The  current  generation  of  statistical  models  of  heterogeneous  propel¬ 
lant  combustion  stems  from  the  original  analyses  of  Hermance  (Reference 
14)  and  of  Beckstead,  Derr  and  Price  (Reference  15).  These  models  are 
"light  on  statistics"  and  heavy  on  mechanistic  modeling  of  the  combustion 
zone.  The  later  developments  involving  more  complete  statistical  descrip¬ 
tion  of  the  combustion  were  initiated  by  Glick  (References  16,  17),  and 
necessarily  embodied  relatively  little  mechanistic  modeling  (but  used  an 
ensemble  assumption  that  drew  from  mechanistic  models  or  experimental 
results  for  mechanistic  inputs).  These  approaches  avoid  direct  confron¬ 
tation  of  the  basic  burning  rate  questions:  "What  is  the  quickest  path 
through  the  propellant  microstructure  for  burning  to  proceed  to  a  given 
point  that  is  distant  from  the  burning  surface  (distant  on  the  scale  of 
ingredient  particle  size).  How  long  does  it  take?"  Without  confronting 
these  questions,  one  has  not  really  directly  confronted  the  processes  that 
determine  burning  rate. 

It  is  not  clear  that  a  direct  confrontation  of  these  questions  would 
be  successful,  but  at  least  it  would  help  bring  out  the  areas  of  ignorance 
that  block  progress  on  the  burning  rate  problem.  It  was  with  this  reali¬ 
zation  that  efforts  were  undertaken  to  develop  a  "least  time  path  model" 
for  describing  burning  rate.  To  minimize  the  statistical  complications 
in  the  first  analysis,  the  model  was  based  on  a  hypothetical  propellant 
made  up  of  oxidizer  spheres  arranged  in  a  lattice  array.  The  statistical 
aspect  of  the  problem  was  introduced  by  allowing  a  limited  range  of  varia¬ 
tion  of  diameter  of  the  spheres,  which  in  turn  resulted  in  a  substantial 
percentage  variation  in  thickness  of  the  binder  interface  between  particles. 
An  analysis  and  computed  results  are  presented  in  the  journal  reprint  ir. 
Appendix  A  (Reference  18).  An  end  result  of  the  analysis  was  a  realization 
that  one  must  know  more  about  the  three-dimensional  aspects  of  heat  flow 
and  surface  regression  in  individual  particles,  before  relevance  of  any 
model  can  be  assured.  Aside  from  these  uncertainties,  the  model  gave 
plausible  burning  rates. 


STABILITY  ANALYSIS  OF  DIFFUSION  FLAMES* 

7.1  Introduction 

During  experimental  work  on  this  program  it  was  observed  that  the 
binder-oxidizer  diffusion  flame  was  rarely  if  ever  steady.  High  speed 
photography  of  combustion  of  sandwiches  showed  the  flame  sheet  to  be  more 
of  a  "brush"  (as  noted  earlier  by  Boggs,  et  al,  References  11,  12).  In 
most  cases  the  individual  flamelets  in  the  brush  were  of  such  short  dura¬ 
tion  that  they  could  not  be  identified  in  successive  frames  of  movies  at 
several  thousand  frames  per  second.  While  a  variety  of  explanations  were 
contrived  for  this  unsteadiness  relating  it  to  unsteady  behavior  of  binder 
melt  or  char  accumulation,  these  explanations  seemed  to  be  of  questionable 
relevance  when  applied  to  the  time  scale  of  unsteadiness  observed  in  the 
movies.  However  the  question  arose  as  to  the  stability  of  the  diffusion 
flame  sheet  itself.  It  was  noted  that  steady  state  solutions  to  the  ana¬ 
lytical  problem  had  been  used  in  the  past  without  verification  that  the 
solutions  were  stable,  and  it  was  decided  to  examine  this  question  further. 
The  following  is  a  brief  summary  of  that  investigation. 

7.2  Posing  the  Analytical  Problem 

Diffusion  flames  can  arise  when  unmixed  fuel  and  oxidizer  flows  come 
into  contact. 

Analytical  modeling  of  steady  state  diffusion  flames  (e.g.,  the  Burke- 
Schumann  Problem,  Reference  19)  usually  include  the  following  assumptions: 

1.  Diffusion  processes  along  the  convective  direction  are  ignored 
(boundary  layer  type  assumption). 

2.  The  region  of  chemical  reaction  is  restricted  to  a  geometrical 
surface. 

3.  The  geometry  of  the  problem  is  simple. 

4.  The  flow  is  considered  to  be  of  low  speed  and  assumptions  are 

*  A  more  detailed  description  of  this  work  is  in  preparation  as  a  thesis 
by  T.  S.  Sheshadri. 
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made  on  the  velocity  and  mass  flux  to  either  automatically  satis  fy 
or  solve  the  mass  and  momentum  conservation  equations. 

The  problem  is  then  solved  by  solving  the  energy  and  species  conservation 
equations . 

However,  there  arises  the  question  of  whether  these  flames  are  dyna¬ 
mically  stable  or  not.  That  is  to  say,  if  disturbed  from  the  steady  state 
position,  will  the  flame  return  to  its  steady  state  position  or  not.  This, 
of  course,  depends  in  part  on  the  nature  of  the  disturbance.  Since  any 
arbitrary  disturbance  can  be  expanded  in  a  Fourier  series,  it  is  sufficient 
to  consider  the  stability  of  the  flame  with  respect  to  a  single  general 
terra  in  such  a  series.  If  such  a  general  term  grows,  decays  or  remains 
constant  with  time  the  flame  is  unstable,  stable  or  neutral  respectively 
with  respect  to  that  general  term.  If  a  flame  is  to  be  stable  for  some 
arbitrary  disturbance  it  must  be  stable  for  all  the  Fourier  component  terms 
that  make  up  the  disturbance.  Thus,  if  the  flame  is  shown  to  be  unstable 
for  a  single  component  term  it  will  be  unstable  for  the  disturbance  as 
a  whole. 

7.3  Analysis 

In  order  to  answer  the  question  of  stability  of  the  steady  state, 
the  unsteady  mass,  momentum,  energy  and  species  conservation  equations 
are  considered.  All  dependent  variables  are  considered  to  be  the  sum  of 
a  steady  state  part  (which  is  known  from  the  steady  state  solution)  and 
a  time  dependent  perturbation  about  it.  The  conservation  equations  are 
then  linearized  in  the  perturbation  quantities.  This  linearization  is 
essential  if  all  arbitrary  disturbances  are  to  be  analyzed  in  a  general 
way.  In  addition,  there  is  the  obvious  advantage  of  simplification  in 
the  mathematics.  As  to  the  relevance  of  the  linear  theory,  the  following 
can  be  said.  If  the  linear  theory  indicates  instability,  so  will  the  non¬ 
linear  theory.  On  the  other  hand  if  the  linear  theory  indicates  stability, 
there  can  still  be  finite  amplitude  unstable  oscillations  as  a  result  of 
finite  initiating  disturbances. 

The  perturbation  quantities  are  taken  as  a  general  Fourier  component 
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term  of  the  arbitrary  disturbance.  Their  behavior  in  time  can  be  obtained 
by  solving  the  linearized  conservation  equations  and  this  answers  (within 
the  scope  of  the  linear  theory),  the  question  of  the  dynamic  stability 
of  the  flame. 

The  above  procedure  was  carried  out  for  a  diffusion  flame  problem 
whose  geometry  differs  somewhat  from  that  of  the  Burke-Schumann  problem, 
and  is  illustrated  below: 

There  are  two  concentric  tubes,  both  of  which  terminate  at  z  =  0. 
Oxidizer  flows  through  the  inner  tube  and  outside  the  outer  tube.  Fuel 
flows  between  the  inner  and  outer  tubes.  Diffusion  flames  can  exist  in 
the  region  z  5  0. 

7.4  Results 

The  analysis  shows  that  if  perturbation  solutions  about  the  steady 
state  exist  at  all,  they  must  exist  in  conjugate  pairs — one  corresponding 
to  an  unstable  solution,  the  other  to  a  stable  one.  Both  of  these  solu¬ 
tions  are  found  to  exist  by  computer  analysis  of  the  problem,  thus  estab¬ 
lishing  the  instability  of  this  class  of  diffusion  flames. 

The  choice  of  geometry  was  motivated  by  the  possibility  of  application 
to  the  combustion  of  composite  solid  propellants,  where  diffusion  flames 
play  a  major  role  over  certain  pressure  ranges.  The  results  of  the  ana¬ 
lysis  suggest  that  the  diffusion  flame  "sheets"  from  the  individual  oxi¬ 
dizer  particles  in  a  propellant  may  be  inherently  unstable.  It  is  too 
much  to  infer  from  the  analysis  what  the  diffusion  flamelets  would  be  like, 
but  the  indicated  inherent  unsteadiness  could  be  very  important  in  low- 
pressure  burning  and  in  oscillatory  burning.  The  results  are  consistent 
with  the  observation  that  the  flame  sheet  in  sandwich-burning  tests  is 
spatially  discontinuous  and  temporally  intermittent.  More  generally,  the 
results  indicate  that  the  classical  solution  to  a  particular  simple  steady 
state  flame  problem  represents  a  physically  unattainable  solution. 


STATISTICAL  MODELING  OF  THE  COMBUSTION  ZONE 


8.1  Introduction 

Description  of  combustion  zone  behavior  is  very  difficult  when  the 
combustion  zone  is  a  chaotic  array  of  surfaces  and  flamelets.  Since  a 
detailed  description  over  a  representative  surface  area  is  virtually  im¬ 
possible,  one  resorts  to  statistical  means  of  describing  what  is  happening. 
Even  in  statistical  descriptions,  there  is  a  need  for  severe  assumptions 
in  the  interest  of  mathematical  tractability  (indeed,  a  statistical  de¬ 
scription  of  the  propellant  microstructure  alone  is  a  formidable  task). 
Early  approaches  made  severe  assumptions  about  statistics  of  the  surface 
microstructure  in  order  to  preserve  some  capacity  to  represent  the  flame 
structure  (References  14,  15).  This  was  followed  by  approaches  that  more 
fully  describe  the  chaotic  array  of  the  burning  surface,  but  contained 
no  modeling  of  the  flame  structure  (References  16,  17).  Influence  of  flame 
structure  was  introduced  with  the  aid  of  results  from  the  earlier  flame- 
structure-oriented  models.  In  the  course  of  the  present  studies,  three 
efforts  were  carried  out  that  pertained  to  statistical  modeling.  One  was 
a  review  of  current  models,  which  was  reported  in  Reference  5.  A  critique 
of  models,  based  on  Reference  5,  is  presented  in  the  following  Section 
(8.2).  This  is  followed  (Section  8.3)  by  a  detailed  description  of  the 

original  surface-statistics-oriented  analysis,  carried  out  primarily  to 

★ 

clarify  the  assumptions  used.  Finally,  Section  6.0  presents  a  new  sta¬ 
tistical  approach  (published  as  Reference  18),  which  is  designed  to  force 
consideration  of  mechanisms  dominant  in  determining  burning  rate,  even 
in  the  process  of  setting  up  the  statistical  approach  (tractability  is 
preserved  in  this  approach  by  limiting  the  disorder  in  the  propellant 
microstructure) . 


*  Since  Section  8  was  written,  significant  extension  of  statistical 

models  has  occurred,  and  the  resulting  models  are  reviewed  in  Refer¬ 
ence  20.  However,  the  following  discussion  remains  relevant,  and 
instructive  to  many  readers.  Accordingly,  it  is  presented  for  "tu¬ 
torial"  purposes  and  to  aid  in  perspective. 
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8.2  Critique  of  Statistical  Models 

For  practical  purposes,  a  description  of  combustion  zone  structure 
is  of  no  interest,  except  to  the  extent  that  it  is  needed  to  describe  or 
control  macroscopic  burning  of  the  propellant.  In  problems  of  steady-state 
burning,  the  detailed  behavior  can  be  averaged  in  both  space  and  time 
(i.e.,  we  are  concerned  with  regression  of  surfaces  that  are  large  compared 
to  the  scale  of  the  propellant  microstructure,  and  with  times  that  are 
long  compared  to  the  burning  time  of  elements  of  the  microstructure). 

In  transient  problems  we  are  concerned  with  the  transient  responses  of 
large  surfaces,  but  often  on  time  scales  comparable  to  (or  shorter  than) 
the  burning  time  of  elements  of  the  microstructure.  Thus  it  is  natural 
that  modelers  would  seek  to  describe  the  microscopic  details  of  the  com¬ 
bustion  zone  in  statistical  terms  that  allow  surface-wise  averaging  of 
surface  and  flame  elements.  Of  course  the  goal  is  to  achieve  surface-wise 
averaging  of  the  behavior  of  elements.  This  usually  starts  with  statis¬ 
tical  descriptions  of  some  classes  of  elements  of  the  combustion  zone  that 
are  considered  to  be  important  to  macroscopic  burning  behavior.  This 
presumes  an  a  priori  knowledge  of  what  determines  the  surface  average 
behavior.  This  last  point  is  extraordinari ly  important,  because  a  sub¬ 
stantial  body  of  analysis  and  "supporting"  experiments  has  developed  on 
this  intuitive  foundation.  It  is  important  that  this  intuitive  foundation 
be  strengthened  by  appropriate  mechanistic  studies;  if  this  is  not  done 
the  evolution  of  a  relevant  statistical  model  will  be  arrested. 

The  foregoing  arguments  are  well  illustrated  by  past  work.  The  early 
work  of  Hermance  (Reference  14)  was  done  at  a  time  when  the  idea  of  attack 
of  oxidizer  vapor  on  binder  surface  or  fuel  vapor  on  oxidizer  surface  was 
strongly  advocated  by  some  investigators.  In  the  Hermance  burning  model, 
such  reactions  were  assumed  to  cause  burning  to  proceed  down  the  interface 
between  oxidizer  and  binder.  This  particular  element  of  the  combustion 
zone  was  found  in  the  analysis  to  be  the  most  important  source  of  heat 
feedback  to  the  solid,  and  hence  the  dominant  factor  in  overall  burning 
rate.  However,  concurrent  experimental  studies  of  burning  of  propellants 
and  oxidizer-binder  "sandwiches"  (Reference  6)  revealed  that  interfacial 
burning  did  not  occur  with  conventional  binders  and  oxidizers,  thus  casting 
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Figure  39.  Combustion  zone  structure  assumed  in  original  Beckstead- 
Derr-Price  model  (Reference  15). 

in  doubt  the  foundation  of  the  "heterogeneous"  theory. 

The  statistical  approach  of  the  Hermance  model  was  combined  with  a 
more  relevant  representation  of  the  combustion  zone  structure  by  Beckstead 
Derr  and  Price  (Reference  15)  and  this  "BDP"  model  has  been  moderately 
successful  in  describing  the  macroscopic  steady  state  combustion  behavior 
of  nona lumini zed  ammonium  perchlorate  propellants  with  hydrocarbon  bindo.  ■ 
This  model  describes  the  combustion  (Figure  39)  in  terms  of  four  kinds 
of  flame  elements, 

Genera  1  Flame  Structure  Assumptions 

1.  Exothermic  decomposition  of  AP  at  its  surface 

2.  Exothermic  gas  phase  flame  of  the  AP 

3.  Primary  diffusion  flame  between  binder  pyrolysis  products  and 
intermediate  decomposition  products  of  AP 

Secondary  diffusion  flame  between  binder  pyrolysis  products 


A. 


In  the  formulation  of  the  original  model,  all  flame  elements  were  tied 
to  the  exposed  surface  size  of  the  oxidizer  particle  by  partly  assumed, 
partly  calculated  flame  structure.  Thus  the  burning  surface  became  a 
statistical  array  of  oxidizer  surface  elements,  each  with  its  own  indepen¬ 
dent  oxidizer  and  diffusion  flame  structure.  Further,  it  was  assumed  that 
the  statistical  array  of  oxidizer  surface  elements  could  be  replaced  by 
an  equal  number  of  elements  of  a  single,  "average"  size,  in  the  process 
replacing  also  the  variety  of  flame  structure  by  an  array  of  identical 
flames.  These  assumptions  are  tabulated  here  because  they  will  come  up 
again. 

Some  Impl ici t  Flame  Structure  Assumptions 

1.  A  unique  flame  structure  is  determined  for  each  size  of  exposed 
oxidizer  surface  (assumed  circular). 

2.  The  flame  structure  in  1  is  indepedent  of  the  time  in  the  burning 
of  the  particle  (except  as  the  size  of  the  exposed  surface  changes 

3.  The  flame  structure  in  1  is  independent  of  the  nature  of  surround¬ 
ing  (or  preceding  or  underlying)  microstructure. 

4.  The  differences  in  flames  of  different  sizes  (insofar  as  their 
effect  on  macroscopic  combustion  behavior  is  concerned)  averages 
in  the  same  way  as  does  oxidizer  surface  area. 

These  were  perfectly  appropriate  strategies  to  use  in  the  interest  of  math¬ 
ematical  tractabi lity,  but  they  represent  further  compromises  with  reality, 
the  impact  of  which  can  only  be  assessed  in  retrospect,  and  then  only  if 
appropriate  mechanistic  studies  are  made  to  evaluate  them  and  motivate 
further  improvement  in  modeling.  Efforts  to  correlate  experimental  b;  in. 
rate  trends  with  predictions  of  the  model  have  suggested  that  the  flames 
associated  with  small  oxidizer  surfaces  must  be  "weighted"  in  the  averaging 
process  because  they  are  more  effective  on  a  per  unit  area  basis  than  large 
oxidizer  surfaces,  i.c.,  to  an  extent  not  balanced  by  the  opposite  effect 
for  surfaces  that  are  above  the  average  size.  This  is  a  plausible  approach 
within  the  constraints  of  present  modeling,  but  the  experimental  trends 
the  revised  model  seeks  to  accomodate  may  in  fact  be  due  to  processes  not 
encompassed  in  the  present  models  at  all. 

A  more  general  treatment  of  the  statistical  aspects  of  the  combustion 
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zone  was  proposed  by  Glick  (Reference  16)  and  developed  further  by  Glick 
and  Condon  (Reference  17),  and  Cohen,  et  al  (Reference  21).  This  treatment, 
as  originally  proposed,  was  concerned  primarily  with  statistical  descrip¬ 
tion  of  the  array  of  surfaces  of  oxidizer  particles,  so  as  to  avoid  settling 
for  an  average  size.  In  this  more  general  statistical  description  the 
implicit  assumptions  1-3  above  are  still  accepted,  but  the  fourth  is  not. 
Instead  it  is  assumed  that  the  surface  regression  rate  of  the  oxidizer 
particle  is  a  unique  function  of  the  exposed  particle  surface  area,  and 
that  the  net  oxidizer  mass  flow  is  the  sum  of  that  from  the  individual 
particle  surfaces,  each  regressing  according  to  its  size.  In  a  limited 
sense,  this  constitutes  the  "weighing"  process  called  for  in  the  last  para¬ 
graph.  In  order  to  find  the  functional  dependence  of  the  oxidizer  surface 
regression  rate  on  the  surface  area,  it  is  assumed  that  the  rate  for  a 
given  size  oxidizer  surface  element  is  the  same  as  the  rate  of  a  propellant 
containing  only  that  size  oxidizer  particle  areas,  as  calculated  by  a  mech¬ 
anistic  model  such  as  the  BDP  model.  This  scheme  provides  a  computational 
procedure  to  calculate  burning  rate  of  the  propellant  as  a  function  of 
oxidizer  particle  size  distribution.  Its  validity  is  unknown,  as  it  con¬ 
tains  a  logical  inconsistency  that  prevents  evaluation.  The  inconsistency 
is  that  small  particles  are  supplied  to  the  propellant  burning  surface 
only  in  the  number  present  in  the  propellant.  If  they  burn  faster  than 
the  larger  particles  (i.e,,  higher  surface  regression  rate),  their  popula¬ 
tion  on  the  surface  will  be  depleted  until  they  contribute  to  the  net 
oxidizer  flow  in  the  proportion  corresponding  to  their  mass  proportion 
of  total  oxidizer  in  the  propellant.  Thus  they  cannot  increase  propellant 
burning  rate  simply  by  burning  faster  individually.  They  must  in  the 
process  increase  the  rate  of  combus tion  of  larger  particles .  At  present 
the  statistical  models  do  not  provide  for  this,  as  reflected  in  Implicit 
Assumption  Number  3.  The  mechanistic  models  contain  a  potential  for  such 
effects,  although  it  is  by  an  indirect  means  relying  on  an  aspect  of  the 
models  discussed  in  the  following. 

In  the  mechanistic  models,  a  second  kind  of  averaging  takes  place, 
disguised  in  the  use  of  "one  dimensional"  heat  transfer  to  the  propellant 
surface.  The  model  seeks  to  determine  the  amount  and  average  site  of  heat 
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release  in  the  various  reaction  steps  noted  above.  The  "site"  of  heat 
release  refers  to  distance  from  the  burning  surface,  and  an  average  dis¬ 
tance  is  calculated  in  order  to  make  a  one  dimensional  heat  transfer  cal¬ 
culation  to  establish  the  effectiveness  of  that  particular  reaction  step 
in  supplying  heat  back  to  the  surface  to  drive  the  combustion  wave.  In 
this  kind  of  calculation,  the  concentrated  sources  are  presumed  to  supply 
heat  impartially  to  all  of  the  surface.  While  current  studies  on  this 
are  seeking  to  deal  with  contributions  from  more  than  one  oxidizer  particle 
size  (thus  averaging  contributions  over  the  particle  size  range),  they 
are  not  examining  the  significance  of  the  three  dimensionally  local  nature 
of  reaction  sites.  Such  an  examination  is  necessary  before  a  rational 
basis  can  be  established  for  how  the  burning  of  small  particles  affects 
the  burning  of  large  particles,  i.e.,  flame  interaction  (elimination  of 
Implicit  Assumption  Number  3). 

Some  modifications  of  the  original  BDP  mechanistic  model  have  been 
made  recently  to  accomodate  for  the  inclusion  of  nitramine  (HMX,  RDX) 
particles  (References  22  -  25)  and  exothermic  binders  (e.g.,  NG/NC)  (Re¬ 
ference  26).  These  modifications  have  not  received  objective  evaluation, 
so  will  not  be  examined  here.  Their  nature  is  chosen  to  deal  with  dif¬ 
ferences  in  deflagration  of  the  "new"  particulate  materials,  and  corres¬ 
ponding  modification  of  flame  structure  resulting  from  the  monopropellant 
nature  of  the  binders,  the  melting  properties  of  ingredients  (and  surface 
flow),  and  the  different  stoichiometry  of  nitramines  and  of  NG/NC  binders. 
Many  of  the  difficulties  with  earlier  models  have  been  confronted  in  the 
more  recent  work  by  Beckstead  (Reference  26),  but  a  final  report  is  not 
currently  available.  In  general,  the  new  work  increases  the  flexibility 
of  the  model  to  deal  with  differences  in  ingredient  characteristics,  but 
of  necessity  deals  superficially  with  detailed  combustion  steps  and  the 
statistics  of  "chaotic"  surfaces. 

As  a  matter  of  perspective,  it  should  be  understood  that  proper  model¬ 
ing  of  the  transient  behavior  of  composite  propellants  is  an  order  of 
magnitude  more  difficult.  Not  only  are  the  Implicit  Flame  Structure  Assump¬ 
tions  of  doubtful  relevance,  but  time-wise  averaging  is  no  longer  applicable. 


The  entire  array  of  flame  elements  and  local  structure  changes  with  time 
during  a  transient,  in  a  manner  that  is  largely  unexplored.  Recent  efforts 
to  adapt  the  statistical  approach  to  description  of  transient  response 
have  not  argued  these  points.  Instead  they  have,  in  effect,  determined 
the  surface-averaged  steady  state  burning  and  then  used  classical  one 
dimensional  theory  to  calculate  transient  response  of  a  homogeneous  pro¬ 
pellant  having  those  steady  state  burning  characteristics.  To  take  care 
of  the  fact  of  a  spread  of  sizes  of  oxidizer  areas  in  the  surface,  the 
one  dimensional  approach  is  applied  over  the  size  range,  and  a  net  tran¬ 
sient  response  is  calculated  as  the  sum  of  the  independent  contributions 
of  each  particle  size.  It  is  extremely  unlikely  that  this  approach  will 
give  valid  results,  as  none  of  the  Implicit  Assumptions  are  realistic, 
and  any  one-dimensionalization  of  heat  feedback  is  also  invalid  when  ap¬ 
plied  on  a  time-resolved  basis.  Further,  there  is  not  enough  known  about 
response  of  the  individual  flame  structural  elements  to  environmental 
disturbances  (e.g.  pressure  oscillations)  to  know  what  to  average  in  order 
to  obtain  the  time  resolved  macroscopic  combustion  zone  response. 

Beyond  these  problems  for  experimentalists  and  modelers  are  the  very 
real  problems  of  how  combustion  characteristics  are  modified  by  certain 
low-concentration  propellant  additives,  and  how  to  model  the  role  of  in¬ 
gredients  such  as  aluminum,  HMX,  or  nitroguanadine,  ingredients  whose 
behavior  is  either  unknown,  or  known  to  be  very  complex. 

8.3  Surface-Statistics-Model 

As  noted  above,  one  line  of  endeavor  in  recent  years  is  aimed  at  a 
more  complete  representation  of  the  statistics  of  the  propellant  hetero¬ 
geneity  in  the  combustion  models.  The  original  approach  by  Glick  is  pre¬ 
sented  here  with  some  elaboration  and  clarification  of  the  assumptions 
used,  or  implied  by  the  approach. 

8.3.1  The  burning  surface  of  a  propellant  is  revealed  by  examination  of 
quenched  samples,  and  shows  a  myriad  of  different  size  oxidizer  surfaces, 
the  different  sizes  reflecting  different  size  oxidizer  particles  and  different 
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times  during  burning  of  individual  particles  (Figure  20).  For  the  present 
purposes,  it  is  assumed  that  each  exposed  oxidizer  particle  can  be  charac¬ 
terized  by  a  surface  area,  sQx  (i.e.,  there  is  a  single  well-defined  con¬ 
tinuous  closed  boundary  between  oxidizer  and  binder  for  each  exposed  oxi¬ 
dizer  particle).  Consequently  there  is  a  size  distribution  function  Fqx, 
for  surface  area  elements  of  oxidizer,  defined  by 
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In  this  expression,  N  represents  the  number  of  all  exposed  oxidizer  surface 
elements  in  a  large  typical  surface  population  of  elements;  dN  is  the  number 
of  those  elements  in  the  surface  area  range  sqx  to  sQx  +  dsQx;  and  Fqx 
is  that  function  of  sqx  that  satisfies  the  equation.  As  a  matter  of  con¬ 
venience  i£  Ls  assumed  that  these  variables  (N  and  F  x)  refer  to  a  unit 
of  planar  propellant  surface  area.  This  normalization  and  the  use  of  dif¬ 
ferential  notation  imply  that  Fqx  is  continuous,  i.e.,  based  on  an  infinite 
population  of  surface  elements,  or  derived  from  a  "smoothed"  fit  to  a 
finite  population.  By  the  same  token,  it  is  then  necessarily  assumed  that 
there  are  no  significant  variations  in  Fqx  on  a  macroscopic  scale  (if  there 
were,  a  characteristic  burning  rate  would  not  exist,  a  problem  that  some¬ 
times  occurs  in  practice  when  propellants  are  not  properly  mixed).  It 
is  appropriate  to  note  here  that  F  can  be  interpreted  as  a  probability 
function,  i.e.,  F  ds  for  a  particular  value  of  s  is  the  probability 
that  a  randomly  chosen  oxidizer  surface  element  will  fall  in  the  size  in¬ 
terval  s  to  s  +  ds  . 
ox  ox  ox 

Dealing  with  the  oxidizer  surface  elements  in  statistical  terms  comes 
more  naturally  than  with  binder  (fuel)  surface  elements,  because  the  oxidizoi 
elements  are  the  more  active  contributors  to  combustion,  are  relatively 
defineable,  and  are  related  to  the  statistics  of  the  original  oxidizer 
particles  used  in  the  propellant.  However,  the  binder  plays  a  role  in 
the  burning  rate  of  the  propellant,  and  thus  needs  to  be  considered  in 
the  statistical  modeling.  Glick  (Reference  16)  suggests  that  the  web  of 
exposed  binder  (fuel)  on  the  burning  surface  be  viewed  as  an  array  of  fuel 
surface  elements,  each  associated  with  a  particular  oxidizer  particle, 
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and  having  a  particular  surface  area,  s ^ .  At  the  moment,  one  might  defer 
yie  exact  procedure  for  identifying  the  fuel  surface  element  or  the  pairing 
rule  relating  s^'s  to  80x's>  and  simply  hypothesize  a  distribution  (pro¬ 
bability)  function  for  fuel  elements  similar  to  that  for  oxidizer. 
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Since  the  fuel  and  oxidizer  surface  elements  are  paired,  N  is  the  same 
for  both  (every  choice  of  an  oxidizer  particle  brings  a  fuel  particle). 
The  distribution  functions  have  the  property  that 
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when  the  integration  is  over  the  full  element  size  range.  This  implies 
simply  that  all  N  particles  are  in  that  range,  or  that  the  probability 
of  a  particle  being  in  that  range  is  one.  Ultimately,  one  is  interested 
in  the  fuel-oxidizer  surface  elements  as  pairs.  The  statistical  distri¬ 
bution  of  size  of  pairs  involves  two  size  variables,  s  and  sc .  If  one 
„  ox  r 

defines  d  N  to  mean  the  number  of  pairs  (out  of  N  total)  with  sqx  between 
sqx  and  sQx  +  dsQx  and  s^  between  s^  and  s^  +  ds^,  then  a  frequency  dis¬ 
tribution  for  pairs  can  be  defined  as 

^77^  -  3ct*.„d<Of  (2) 

where  the  script  form  of  3  denotes  a  distribution  function  for  pairs.  This 
is  depicted  in  Figure  40.  In  this  definition,  it  is  assumed  that  N  and 
3  refer  to  the  population  per  unit  planar  surface  area,  analogous  with 
the  definitions  of  FQX  and  F^.  In  the  figure,  d  N/N  is  the  volume  under 
the  surface  element  dsQxds£,  the  probability  that  a  randomly  chosen  pair 
will  be  in  the  particular  interval  dsQxdSj.  Double  integration  over  the 
entire  size  range  present  will  give  all  the  particle  pairs,  i.e.,  the 
integral  of  d  N/N  equals  one;  this  corresponds  to  a  volume  of  one  under 
the  surface  (probability  of  one  that  a  randomly  selected  pair  from  the 
population  will  be  a  member  of  the  population). 
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Figure  40.  Combined  distribution  function  for  oxidizer  and  fuel  elements 
(Equation  2). 

The  function «7  is  related  to  the  functions  F  and  Fr.  Thus  if  one 

ox  f 

considers  a  specific  value  of  s  .it  relates  to  a  curve  on  the  (s  ,  s,) 
r  ox’  ox  t 

surface,  parallel  to  the  sQx  =  0  coordinates  plane  (Figure  40).  The 
volume  under  a  differential  slice  of  the  &  surface  based  on  the  constant 
sqx  curve  is  the  proportion  of  oxidizer  particles,  dN/N  in  the  interval 
ds^  .  In  other  words 
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So7ds^  is  the  probability  that  a  random  sample  will  have  a  value  of  s£ 
between  sqx  +  dsQx  regardless  of  value  of  s^  (and  analogous  statement  for 
JHs  ). 
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The  analytical  form  of  the  function  7  is  determined  by  the  distribution 
functions  of  sqx  and  s^;  but  any  specific  distribution  function, J7,  depends 
on  a  pairing  rule  for  fuel  and  oxidizer  surface  elements.  Do  large  oxidi¬ 
zer  surface  elements  tend  to  occur  in  combination  (pair)  with  correspond¬ 
ingly  large  fuel  surface  elements,  and  small  sqx  with  small  s^?;  large 
8ox  with  small  s^?;  random  combinations?  To  answer  this  question  ration- 
ally,  one  would  have  to  go  back  and  define  fuel  surface  elements,  presumably 
in  a  manner  based  on  the  real  issue,  characterization  of  an  element  of 
the  statistical  flame  population  (References  16  and  26).  For  most  propel¬ 
lants,  the  particulate  content  is  so  high  that  random  packing  of  particles 
in  the  solid  is  impossible,  so  a  valid  rule  for  pairing  of  oxidizer  and 
fuel  surface  area  elements  would  have  to  reflect  the  fact  that  small  par¬ 
ticles  pack  in  the  "voids”  among  large  particles.  Thus  fuel  surface  ele¬ 
ments  would  have  to  reflect  this  fact  in  their  definitions  and  in  the 
resulting  rules  regarding  how  the  elements  pair  with  the  oxidizer  surface 
elements  to  give  the  distribution  functions  Fqx,  Ff  and  3.  This  aspect 
of  the  statistical  models  has  been  circumvented  almost  without  comment 
in  most  modeling  work,  simply  by  making  unevaluated  assumptions  (see  com¬ 
ments  in  Reference  26). 

Continuing  on  the  approach  used  by  Click  (Reference  15),  the  pairing 
problem  was  "solved"  by  assuming  that  the  pairing  was  a  matter  of  pure 
chance;  then  for  a  particular  increment,  dsQx,  of  the  FQX  population,  the 
number  of  those  oxidizer  elements  that  are  paired  with  a  particular  incre¬ 
ment,  dsj,  of  fuel  elements  is  determined  by  the  proportion,  F^,  of  the 
fuel  elements  in  that  particular  ds^  increment. 

Analytically  the  foregoing  restrictive  assumption,  i.e.,  random  pair- 
ing,  means  that  all  constant-sQx  curves  on  the  distribution  surface  are 
functionally  alike  and  ail  the  constant-s^  curves  are  functionally  alike, 
which  implies  that  (*ox»  s^)  has  the  form  (Reference  27 ) 
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(in  the  following  the  functional  notation  will  be  dropped,  leaving^  with 
the  subscript  "f"  referring  to  a  function  of  s^,  or f  with  the  subscript 
"ox"  referring  to  a  function  of  sqx).  Thus,  for  any  "random  pairing" 
distribution  surface,  <7,  the  function  is  a  simple  function  of  s^  only, 
and fQX  is  a  simple  function  of  sqx  only. 

For  a  distribution  function  of  the  form  above  (Equation  4),  the  dis¬ 
tribution  fQX  and  are  simply  related  to  FQx  and  Fx-  Thus 
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and  following  the  form  of  Equation  3, 
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This  implies  that 
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and  substituting  in  Equations  4  and  5 
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Recalling  that  double  integration  of  Equation  5  yields  a  volume  of  1  under 
the  surface  7  (probability  of  1),  Equation  5  means  that 
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Then  the  denominator  in  the  bracket  in  Equation  9  has  a  value  of  one,  and 

—  -  J  ct*f  =  /,  do,  4*.^  ~  (ID 

In  other  words,  for  the  special  form  of  the  distribution  function  in  Equa¬ 
tion  4, 

J  ~  /-t  foy  -  /=f  /Zy  (12) 

This  is  the  form  of  the  distribution  function  used  originally  by  Glick 
(Reference  15).  When  the  realities  of  defining  and  pairing  of  fuel  and 
oxidizer  surface  elements  are  more  thoroughly  examined,  and  applied  to 
situations  such  as  bimodal  oxidizer  particle  size  distributions,  this 
assumption  of  random  pairing  will  need  review  and  may  have  to  be  modified. 
However,  it  will  be  used  for  the  present,  and  attention  will  be  turned 
to  other  considerations. 

8.3.2  Mass  flow  from  the  burning  surface  is  the  objective  of  most  ana¬ 


lytical  models  of  solid  propellant  combustion,  Indeed,  the  singular  nature 
of  the  surface  (as  a  boundary  separating  domains  of  physical  and  chemical 
states  and  of  transitions  of  physical  and  chemical  processes)  traditionally 
leads  to  analytical  formulation  of  conservation  laws  for  the  surface. 

These  are  used  as  boundary  conditions  for  joining  models  of  processes  in 
the  solid  with  processes  in  the  gas.  This  convention  is  obviously  jeopar¬ 
dized  by  presence  of  a  chemically  heterogeneous,  geometrically  complex 
surface,  and  the  principal  role  of  the  statistical  model  is  to  represent 
the  real  process  in  some  average  form.  Early  approaches  (References  14, 

15)  sought  to  achieve  this  by  replacing  the  distribution  of  sizes  of  oxi¬ 
dizer  surface  elements  by  that  resulting  from  some  "effective"  single  size 
spherical  particles.  Since  this  still  gave  a  distribution  of  oxidizer 
surface  element  areas  (different  times  during  burning  of  the  "monodispersc" 
particles),  a  further  approximation  was  made,  i.e.,  all  surface  elements 
were  the  same  size,  that  size  being  the  time-average  area  for  burning  of 
a  particle.  This  purely  hypothetical  surface  was  then  used  to  write  mass 
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conservation  equations.  The  mass  burning  rate  (per  unit  area)  on  the 
surface  elements  was  considered  to  be  a  function  of  surface  element  size, 
but  only  after  arriving  at  the  hypothetical  "equal  circle"  representation 
of  the  surface  (i.e.,  only  for  different  propellants).  As  noted  by  Glick, 
this  does  not  take  into  account  the  fact  that  the  mass  burning  rate  per 
unit  area  actually  is  not  expected  to  be  the  same  for  all  surface  elements 
because  of  actual  size  difference.  If  an  average  procedure  is  to  be  used 
in  writing  the  conservation  laws,  the  average  must  be  an  average  of  the 
quantity  being  conserved;  in  this  case,  mass  flux,  not  area.  Glick  went 
to  some  lengths  to  stress  this.  However,  the  situation  is  complex,  and 
successive  studies  will  no  doubt  continue  to  reveal  shortcomings  in  prior 
analyses.  For  example,  the  original  analysis  by  Glick  did  not  consider 

1.  The  dependence  of  mass  flux  (flow  per  unit  area)  from  surface 
elements  of  the  same  size  at  different  stages  in  the  particle 
burning  (®ox  is  assumed  to  be  a  function  of  sox,  with  reference 
to  size  of  particle  or  stage  of  its  burning). 

2.  The  nonuni formi ty  of  mass  flux  as  a  function  of  location  on  the 
element  surface. 

3.  The  imprecise  definition  of  the  binder  surface  elements  and 
corresponding  difficulty  of  deciding  anything  about  mass  flux 
from  them. 

Regardless  of  these  unresolved  problems,  it  is  no  doubt  worthwhile 

to  pursue  the  formulation  of  mass  conservation  equations  within  this  limited 

context,  and  generalize  from  there  as  familiarity  with  the  analysis  and 

advances  in  mechanistic  knowledge  subsequently  permit.  Thus,  following 

the  general  approach  of  Glick,  the  mass  flow  rate  per  unit  area  of  oxidizer 
2 

produced  from  the  d  N  surface  elements  in  the  size  range  s  to  s  +  ds 

ox  ox  ox 

and  s^  to  s^  ♦  ds^  is 

d2/mox  =  '*>\xCt4ltllc/*/>/-/>r?0%‘<kxfr/SG{49Xc/4f')  (13) 

Thus  contributions  to  flux  from  pairs  of  surface  elements  in  which  the 

oxidizer  surface  elements  are  in  the  range  s _  to  s  ♦  ds  (but  all  s, 
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sizes)  is 
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and  the  total  oxidizer  flux  per  unit  planar  surface  area  is 
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In  these  integrations  it  should  be  noted  that  m  remains  under  the 

ox 


inner  integral,  reflecting  its  possible  dependence  not  only  on  oxidizer 
surface  element  size,  but  also  on  fuel  surface  element  size.  In  Equation 
15,  no  use  has  been  made  of  the  random  pairing  assumption  used  by  Glick. 
If  the  assumption  is  used,  Equations  12  and  15  give  Glick’ s  equation  7 
(Reference  15) 


i0  x  “  A  J^OX 


^xf^ox 


Further,  if  the  dependence  of  m^  on  fuel  surface  element  size  in  the  pair 


is  neglected,  this  simplifies  further  to 

””*x  --  N  J"»ox-6-.x  Fox  frf  dxifdxi' 

\x  “  A/ /Tr  ^  V 


*x 


(17) 


In  a  practical  sense,  this  last  assumption  may  not  represent  a  serious 
further  departure  from  reality,  as  the  dependence  of  oxidizer  product  flux 
on  fuel  element  size  cannot  be  realistically  accounted  for  in  Equation 
16  anyway,  because  the  oxidizer  surface  elements  have  not  been  realistical 
paired  with  the  fuel  surface  elements  by  the  random  pairing  assumption 
used.  All  three  of  the  forms  15  -  17  may  merit  future  consideration, 
depending  on  the  particular  situation  and  application. 

It  should  be  noted  that  a  parallel  development  can  be  made  to  obtain 
symmetrical  equations  for  fuel  product  flux.  If  needed,  those  equations 
can  be  written  from  Equations  14-17  by  interchanging  subscripts  ox  and 
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f.  As  noted  before,  calculations  based  on  fuel  product  flux  are  usually 
not  used  because  of  poor  definition  of  the  fuel  surface  elements  and  rela¬ 
tively  passive  role  of  the  fuel  surface.  This  may  not  always  be  so  in 
future  developments.  However,  it  ijs  necessary  that  the  oxidizer  and  fuel 
product  mass  flows  be  in  a  ratio  dictated  by  the  mixture  ratio  of  binder 
and  oxidizer  in  the  propellant.  Thus  if 


sPrX  j* 

£  aor  ’  -  ur  / 

oce„  +  ocf  =  / 

(18) 

<r+\T 

<*ox  ~  OCf 

/-xoX 

(19) 

The  quantity  m^,  is  the  net  mass  flow  rate  per  unit  planar  burning  surface 
area,  the  a's  are  a  result  of  the  propellant  formulation,  and  suitable 
representation  of  mQx  (or  m^)  can  be  obtained  from  the  foregoing  (e.g.. 
Equations  15-17). 

8.3.3  Application  of  Equations  17  and  18  provides  a  scheme  for  calculating 
mass  burning  rate  per  unit  planar  surface  area,  and  hence  linear  surface 
regression  rate.  The  actual  calculation  is  dependent  on  knowing  the  vari¬ 
ables  on  the  right  in  Equations  15-17.  "Knowing  the  variables"  means 
knowing  the  functional  forms  necessary  to  carry  out  the  integrations. 

Using  the  results  depends,  of  course,  on  knowing  the  dependence  mQx  on 
other  variables  of  interest,  variables  not  explicitly  shown  in  the  equations, 
variables  such  as  fuel-oxidizer  pairing,  propellant  temperature,  environ¬ 
mental  pressure,  and  deflagration  characteristics  of  the  oxidizer.  In 
short,  Equations  15-17  and  19  constitute  a  statistical  statement  of  the 
mass  conservation  at  the  burning  surface.  All  other  important  aspects 
of  the  combustion  process  must  be  independently  modeled  in  a  form  giving 
®ox  as  a  ^uncti°n  of  all  relevant  variables,  including  specifically  those 
variables  (sqx,  s^)  needed  to  carry  out  the  indicated  integrations. 

In  addition  to  knowing  mQx  (sQX,  sf)  to  evaluate  mQx  and  n^,  the 


distribution  function  of  surface  element  pairs  must  be  determined.  This 
is  ordinarily  thought  of  as  deriving  from  the  oxidizer  particle  size  dis¬ 
tribution,  and  has  been  based  entirely  on  that  in  some  analytical  models 
(e.g.,  References  20-23).  As  evidenced  by  the  foregoing  (and  the  analysis 
of  Glick),  a  more  rigorous  statistical  model  of  the  surface  is  not  affected 
by  the  exact  meaning  of  sqx  and  s^,  at  least  up  to  Equation  19  inclusive. 

The  "day  of  reckoning"  comes  when  these  vagaries  focus  into  a  determination 
of  mQx  vs  sqx>  s ^  pairs,  a  determination  that  cannot  be  made  rigorously 
without  reference  to  the  complex  micro-scale  processes  in  the  combustion 
zone  that  determine  the  surface  details  discussed  in  the  present  report. 

These  pitfalls  may  be  "boldly"  circumvented  by  making  direct  assumptions 
about  the  surface  geometry  (all  models  to  date  do);  and  assumptions  about 
dependence  of  mQx  on  particle  size,  about  pairing  of  sqx  and  s^,  and  about 
other  variables  such  as  pressure.  It  must  be  anticipated  that  such  reck¬ 
less  appearing  assumptions  will  be  tried,  in  the  interest  of  mathematical 
tractabi lity .  In  a  sense,  they  are  a  logical  step  in  the  search  for  un¬ 
derstanding  and  useful  results. 

Those  investigators  who  have  use  of  the  foregoing  statistical  approach 
have  obtained  the  dependence  of  mQx  on  pressure  and  temperature  from  the 
results  of  earlier  published  models  (References  14,  15).  This  evades 
direct  consideration  of  energy  conservation  and  combustion  processes  in 
general,  bringing  instead  a  more  visible  treatment  of  the  statistics  of 
the  burning  surface. 

In  the  interest  of  early  progress  towards  realistic  modeling,  it  is 
helpful  to  have  clearly  in  mind  the  nature  of  the  assumptions  in  present 
models.  The  following  summarizes  those  assumptions  explicitly.  There 
is  a  growing  awareness  that  some  of  these  assumptions,  e.g.,  assumptions 
5-9  are  unacceptable,  and  recent  work  has  begun  to  provide  more  realistic 
representation  of  the  detailed  combustion  process  involved. 

8.3.4  Assumptions  in  the  Surface  Statistics  Model 
1.  It  is  assumed  that  exposed  oxidizer  particles  each  exhibit  a  define- 

able,  connected  burning  surface  area.  No  assumption  of  shape  of  surface 
is  made.  It  is  implied  that  the  surface  elements  have  areas  related 
to  particle  size. 
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2.  It  is  assumed  that  the  exposed  binder  (fuel)  surface  area  can  be 
allocated  in  some  manner  into  surface  elements  assignable  one-on-one 
to  the  oxidizer  surface  elements  (although  no  unique  procedure  is 
specified  to  do  this). 

3.  It  is  assumed  that  the  microscopic  heterogeneity  is  macroscopically 
uniform,  so  that  a  regression  rate  of  the  surface  (burning  rate)  can 
be  assigned  that  is  uniform  on  a  macroscopic  scale  (i.e.,  when  aver¬ 
aged  over  many  particle  diameters),  and  characteristic  of  the  pro¬ 
pellant. 

4.  Based  on  3,  it  is  assumed  that  the  burning  surface  can  be  character¬ 
ized  by  an  average  number,  N,  of  oxidizer  surface  elements  per  unit 
burning  surface  area  (the  same  number  of  fuel  surface  elements). 

In  this  assumption,  the  surface  elements  may  be  non-flat,  but  the 
number  N  is  the  average  number  of  such  elements  per  unit  planar  sur¬ 
face  area,  i.e.,  area  projected  on  a  plane  perpendicular  to  the  di¬ 
rection  of  the  mean  burning  rate. 

5.  In  certain  instances,  it  is  assumed  that  the  pairing  of  oxidizer 
surface  elements  and  fuel  surface  elements  is  uncorrelated  with  ele¬ 
ment  sizes;  i.e.,  the  probability  of  a  given  size  fuel  surface  element 
pairing  with  any  particular  oxidizer  surface  element  depends  only 

on  the  relative  frequency  of  occurrence  of  that  size  in  the  oxidizer 
population. 

6.  The  mass  flow  rate  per  unit  real  surface  area  (flux)  of  an  oxidizer 
surface  element  is  single  valued  over  the  surface  of  the  element  (at 
any  moment) . 

7.  The  mass  flux  of  oxidizer  products  may  be  different  for  different 
size  oxidizer  surface  elements. 

8.  The  mass  flux  of  oxidizer  products  from  a  given  size  oxidizer  surface 
element  may  depend  on  the  size  of  the  companion  fuel  surface  element 
(Equation  16),  or  may  not  (Equation  17). 

9.  The  mass  flux  of  oxidizer  products  from  a  given  size  oxidizer  surface 
element  is  not  dependent  on  the  circumstances  leading  to  the  element 
size  (i.e.,  no  distinction  is  made  between  particle  size  contribution 
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to  element  size  and  time-during- burning  contribution  to  element  size) 
This  implies  that  no  consideration  is  given  to  the  local  microstruc- 
tural  details  except  as  they  are  reflected  in  surface  element  area. 

10.  It  is  assumed  that  the  dependence  of  burning  rate  on  pressure,  pro¬ 
pellant  temperature,  and  other  variables  not  explicitly  contained 

in  the  statistical  model  can  be  brought  in  as  rate  dependence  of  the 
oxidizer  surface  element. 

11.  It  is  assumed  that  a  surface  element  size  distribution,  ,  can  be 
obtained  for  the  oxidizer-fuel  surface  elements  pairs,  either  by  a 
more  detailed  theory,  or  by  direct  observation  of  burning  surfaces, 
or  by  further  simplifying  assumptions. 

Recent  developments  in  modeling  have  sought  to  remedy  (or  evaluate) 
the  limitations  of  assumptions,  and  these  efforts  were  discussed  in  a 
JANNAF  Workshop  and  reported  in  an  excellent  review  by  the  chairman  of 
the  workshop  (Reference  20).  A  very  significant  effort  has  been  put  into 
design  of  better  physical  and  statistical  models,  with  heavy  reliance  on 
rather  extensive  computer  programs  to  extract  burning  rate  predictions. 

The  reader  is  referred  to  Reference  20  for  the  status  of  results.  It  is 
notable  that  "improvements"  in  the  models  (i.e.,  to  conform  to  reality) 
are  usually  conceived  on  the  basis  of  recognized  deficiencies  in  repre¬ 
sentation  of  detailed  physio-chemical  processes  operating  at  the  micro¬ 
scopic  level,  but  the  adequacy  of  the  improvements  in  the  model  are  judged 
by  the  improved  ability  of  the  model  to  correlate  measured  average  burning 
rates.  This  is  not  a  particularly  stringent  test  of  the  validity  of  the 
modified  models,  as  they  usually  have  greater  flexibility  in  functional 
form  and  "free"  parameters.  As  a  result,  a  variety  of  models  have  emerged, 
but  relatively  little  work  has  been  done  at  the  microscopic  level  to  eval¬ 
uate  or  guide  the  modeling. 


Notation 


average  number  of  oxidizer  surface  elements  (per  unit  planar 
surface  area)  having  area  in  the  interval  s  to  s  +  ds 
(see  Eq.  6)  ox  ox  ox 

average  number  of  element  pairs  sqx,  s^  (per  unit  planar  surface 

area)  having  oxidizer  element  area  between  s  and  s  ♦  ds 
and  fuel  element  area  between  s^  and  s^  +  ds^X(see  Eq*  5)  °X 

normalized  frequency  distribution  of  one  variable 

normalized  frequency  distribution  of  oxidizer  surface  elements 

normalized  frequency  distribution  of  s^  (see  Eq.  16) 

normalized  frequency  distribution  of  two  variables,  s  and 
sf  (see  Eq.  2)  OX 

a  function  of  sQX,  proportional  to  FQx  (see  Eq.  4) 

a  function  of  sf,  proportional  to  (see  Eq.  4) 

mass  flow  rate  per  unit  real  area,  mass  flux 

mass  flux  averaged  over  a  range  of  surface  element  sizes; 
per  unit  planar  surface  area 

mass  flow  rate  per  unit  real  area,  of  oxidizer  products  from 
oxidizer  surface  area  (dependent  on  s  ,  sf,  possibly  other 
variables  such  as  pressure)  OX 

mass  flux  from  fuel  surface 

mass  flux  from  propellant  surface  mass  flow  per  unit  planar 
surface  area 

N  average  number  of  surface  element  pairs  per  unit  planar  surface 

area 

Or  ratio  of  mass  of  ingredients  in  the  solid  propellant 

0/  ratio  of  oxidizer  mass  to  total  mass  in  the  solid  propellant 

Of  ratio  of  fuel  mass  to  total  mass  in  the  solid  propellant 


*  +  i 

m 
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APPENDIX  A 

Some  Statistical  Considerations  in  the  Burning  of 
Composite  Solid  Propellants 


Warren  C.  Strahle* 

Georgia  Institute  of  Technology,  Atlanta,  Ga. 


Propagation  of  ■  deflagration  wave  through  a  near!)  ordered  composite  propellant  is  considered.  The  intent  is 
to  see  some  effects  of  the  particle  packing  statistics  on  the  propagation  rate.  A  highly  oversimplified  model  of 
the  deflagration  physics  is  considered,  and  only  the  statistics  are  emphasized.  It  is  found  that,  for  propellants  of 
usual  packing  density,  the  least-lime  path  of  burning  to  a  given  point  depends  primarily  on  the  packing  statistics 
in  a  line  parallel  to  the  burn  rate  vector;  there  is  little  effect  of  particles  only  a  few  particle  diameters  to  the  side 
of  this  line.  In  a  simplified  model  that  may  be  solved  analytically,  there  arc  clearly  seen  particle-size,  pressure, 
and  packing  density  effects  on  bum  rale  which  are  due  to  the  statistics  alone  and  not  the  deflagration  physics. 


Introduction 

HE  heterogeneity  of  a  composite  solid  propellant  gives 
severe  problems  when  it  comes  to  modeling  the 
combustion  behavior.  If  one  is  positioned  on  a  line  parallel  to 
the  mean  burn  rate  vector  and  follows  the  surface  as  it 
regresses,  a  time-dependent  process  is  seen  by  the  observer, 
involving  alternate  binder  and  oxidizer  burnthrough  (ex¬ 
cluding  metalized  propellants)  and  other  unsteady  processes 
such  as  bubbling  melts  and  melt  flows.  Alternatively,  at  any 
instant  of  time,  if  the  observer  looks  over  the  surface,  a 
spacewisc  heterogeneous  process  is  involved.  Precision  in 
modeling,  therefore,  requires  consideration  of  both  time 
dependence  and  space  heterogeneity. 

Glick 1  and  Glick  and  Condon2  have  tackled  the  space 
heterogeneity  problem  by  application  of  a  statistical  method 
to  the  propellant  structure.  In  their  method,  a  combustion 
model  is  coupled  with  a  statistical  description  of  various  fuel- 
oxidizer  pair  sizes  to  yield  an  average  burn  rate.  The  com¬ 
bustion  model  favored  is  the  BDP3  model,  which  is  essentially 
a  steady-state  model  and  ignores  the  transients  that  would  be 
found  in  an  actual  propellant  burnthrough  problem.  There 
are  also  difficulties  in  the  statistical  description  of  the 
propellant  surface.  An  early  assumption  in  Ref.  1  is  that  for 
any  particle  of  oxidizer  the  probability  of  a  particular  "size” 
of  binder  pocket  being  adjacent  to  the  oxidizer  is  independent 
of  the  oxidizer  size.  This  assumption  is  tantamount  to 
assuming  complete  disorder  to  the  propellant  array.  Fur¬ 
thermore,  it  is  assumed  that  the  existence  of  a  particular 
oxidizer  size  at  one  point  does  not  influence  the  allowable 
particle  size  at  a  neighboring  point;  this  again  is  an  assump¬ 
tion  of  complete  disorder  to  the  propellant  structure. 

These  two  difficulties,  unsteadiness  and  a  statistical 
description  of  the  propellant  surface,  actually  are  coupled. 
The  deflagration  physics  will,  in  fact,  affect  the  surface 
configuration.  It  therefore  would  seem  reasonable  to  actually 
track  the  deflagration  front  in  time,  using  an  appropriate 
deflagration  model,  to  see  how  the  surface  details  unfold  with 
time.  There  are  obvious  difficulties  with  such  an  approach, 
however,  because  one  still  would  need  a  combustion  model 
and  a  statistical  description  of  the  propellant  packing. 

Concerning  the  propellant  packing  statistics,  there  is  a  real 
issue  concerning  the  degree  of  disorder  to  the  oxidizer  par¬ 
ticles  imbedded  in  a  binder  matrix.  For  nonspherical  particles 
with  multimodal  particle  size  distributions  and  wide  cuts 
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about  a  nominal  particle  size,  the  observation  of  Ref.  2  that  a 
propellant  is  a  quite  disordered  structure  is  probably  correct, 
even  if  the  propellant  is  near  the  maximum  possible  solids 
loading.  On  the  other  hand,  it  is  possible  to  conceive  of  a 
nearly  ordered  propellant  structure  if  single-particle-size, 
spherical  particles  are  packed  at  nearly  the  maximum  solids 
loading.  This  becomes  a  completely  ordered  structure  in  the 
limit  that  the  oxidizer  loading  is  a  maximum  and  the  particle- 
size  cut  is  of  zero  variance  about  the  mean  particle  size  (a 
delta-function  size  distribution).  Such  a  packing  is  shown  in 
Fig.  I.  The  lattice  points  of  the  array  form  an  oblique 
parallelpiped.  A  perfectly  ordered  structure  also  would  result 
with  spherical  particles  in  a  multimodal  distribution  at  the 
maximum  solids  loading  if  the  particles  had  narrow  cuts 
about  each  size  in  the  distribution.  The  propellant  of  Fig.  I  is, 
in  principle,  possible  to  manufacture,  although  it  would  not 
be  a  practical  propellant.  If  one  now  backs  off  slightly  from 
the  maximum  packing  density,  a  nearly  ordered  structure 
would  result  with  particles  near  but  randomly  displaced  from 
their  ordered  lattice  points.  This  is  the  kind  of  propellant  to  be 
considered  in  this  paper. 

The  motivation  for  considering  the  nearly  ordered  structure 
is  twofold.  First,  it  is  a  simple  structure  through  which  to 
track  an  unsteady  deflagration  wave.  Second,  the  statistics  of 
the  packing  are  sufficiently  simple  to  investigate  the  question 
of  what  is  the  lateral  extent  of  influence  of  one  point  upon  the 
other.  Stated  otherwise,  given  the  deflagration  arrival  at  a 
point,  how  was  the  time  of  arrival  affected  by  particles  above, 
but  transversely  separated  from,  the  vertical  axis  running 
through  the  particle  of  interest? 
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Although  the  propellant  of  Fig.  I  is  clearly  not  a  practical 
propellant,  it  is  intended  for  use  to  see  some  statistical  effects 
that  are  surely  present  in  real  propellants.  The  model  to  be 
employed  for  the  deflagration  process  is  also  highly  over¬ 
simplified  and  is  selected  not  for  its  deflagration  physics 
credibility  but  its  demonstration  of  statistical  effects. 

This  work  should  be  considered  complementary  to  the  work 
of  Refs.  I  and  2;  this  is  intended  to  answer  several  questions 
that  cannot  be  answered  at  the  present  time  by  their 
procedure.  Although  the  analogy  is  not  precise,  the  work  here 
may  be  considered  a  Lagrangian  approach  to  the  problem 
where  the  deflagration  front  is  tracked  in  time.  The  approach 
of  Refs.  I  and  2  may  be  likened  to  an  Eulerian  view  where  a 
time-stationary  average  surface  is  under  observation. 

Analysis 


Lattice  Statistics 

The  lattice  structure  of  Fig.  I  is  the  starting  point  and 
represents  the  maximum  packing  density  for  spherical  par¬ 
ticles.  If  V0  is  any  volume  containing  a  targe  number  of 
particles,  the  volume  V  of  particulate  material  is  given  from 
elementary  geometry  as  V/V0  = 0.6981,  which  is  the 
maximum  packing  density.  Consider  keeping  the  same  lattice 
geometry,  but  shrinking  each  particle  in  size  by  an  equal 
amount.  Then  each  particle  is  in  an  ordered  array  with 
VI  Vu  <0.6981,  and  the  particles  are  separated  from  each 
other  by  an  equal  amount  at  their  points  of  closest  approach. 
This  perfectly  ordered  propellant  actually  would  have  a  mild 
anistropic  burning  behavior.  In  any  principal  lattice  plane  and 
along  a  direction  parallel  to  a  principal  lattice  axis,  everything 
looks  the  same.  But,  by  making  planar  cuts  in  an  arbitrary 
plane  and  moving  along  an  arbitrary  line,  the  configuration 
changes  to  an  observer,  as  compared  with  observations  on  the 
principal  planes  and  directions.  As  an  assumption  of  the 
analysis,  observation  will  only  be  made  on  a  principal  plane 
and  along  a  principal  direction. 

Now  consider  Fig.  2a,  where  a  typical  three-particle  element 
is  viewed,  and  consider  small  random  displacements  of  the 
particles  from  their  ordered  lattice  points.  The  particle  radius 
is  r,  the  distance  between  ordered  lattice  points  is  d,  the 
displacement  of  the  center  of  any  particle  from  its  ordered 
position  is  the  vector  b  and  the  distance  between  two  particles 
along  the  line  of  centers  is  x.  The  statistics  to  be  employed 
consist  of  assigning  a  probability  distribution  to  the  i’s  and 
then  following  the  burn  through  the  matrix. 

To  simplify  things  further,  the  “one-dimensional”  ap¬ 
proximation  will  be  made  that  4  has  only  a  vertical  component 
6  and  that  all  planes  connecting  particle  centers  which  are 
parallel  to  the  paper  in  Fig.  2  have  the  same  set  of  Vs. 
Consequently,  only  events  in  the  plane  of  the  paper  are  of 
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Fig.  2  Configuration  of  a 
nearly  ordered  array  of 
spherical  part  teles. 
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interest  (there  is  no  interaction  between  adjacent  planes).  One 
now  has  the  situation  of  Fig.  lb.  From  elementary  geometry, 
the  distance  A -A',  which  is  x  for  particles  1  and  2,  is 

x-=  |  (3d2/4)  +  \(d/2)  -bt -b,)2),l2-2r  (1) 

An  approximation  to  Eq.  (1),  which  will  be  used  and  is  valid  if 
(d—2r)/d<  1,  is 

*=d-2/-+[(6/-6,)/21»A  +  [(6(-6i)/2|  '  (2) 

Equation  (2)  is  valid  for  well-packed  propellants  if  the  Vs  are 
also  restricted  so  that  6,  /d<  I . 

There  is,  in  principle,  no  restriction  on  the  magnitude  of  6. 
However,  if  one  draws  a  sphere  (circle)  of  influence  about  the 
ordered  lattice  center  of  r+6,„,  where  6,„  is  the  maximum 
expected  6,  it  is  clear  that  these  spheres  will  overlap  if 
Sm  >  (d-2r)/2.  Consequently,  if  6m>(d-2r)/2,  one 
particle  will  interfere  with  the  allowable  positions  of  an  ad¬ 
jacent  particlefs).  This  again  would  complicate  the  statistics, 
so  it  is  required  that  6OT  s  (rf-2r)/2.  This  may  be  viewed 
physically  as  a  “mixedness”  assumption.  That  is,  if  the 
propellant  is  well  mixed,  no  large  voids  of  particles  would  be 
expected  and  each  particle  would  be  near  its  ordered  lattice 
point. 

For  illustrative  purposes,  the  probability  distribution  for 
the  Vs  will  be  the  uniform  distribution,  whereby  the  dif¬ 
ferential  probability  dPof  Finding  6  between  6  and  6  +  d6  is 

dP=db/2bm  (3) 

Obviously  the  particle  must  be  somewhere  between  -  bm  and 
+  6,„,  so  that 


P[-b„  s6s6m) 


d  P=l 


(4) 


Finally,  two  cases  will  be  considered.  The  first, 
corresponding  to  Fig.  2b,  will  be  periodic  about  two  column 
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widths  and  will  be  symmetric  about  the  left-hand  centerline  on 
Fig.  2b.  This  means  that  only  the  region  between  the  two 
centerlines  need  be  considered  for  burning  calculations.  The 
second  case  is  shown  in  Fig.  3a.  This  case  considers  symmetry 
about  the  left-hand  centerline  and  periodicity  about  three 
column  widths.  Thus,  in  case  1,  only  two  columns  can  in¬ 
teract.  In  case  2.  three  columns  can  interact.  In  principle,  one 
could  continue  (his  escalation  to  an  infinite  number  of  column 
widths  which  would  be  the  limit  of  a  truly  heterogeneous, 
random  propellant.  It  will  be  found  that  this  is  not  necessary. 

Mechanistic  Model 

It  is  assumed  first  that  the  oxidizer  particles  are 
monopropellants  that  have  a  planar  deflagration  rate  of  r0  at 
pressure  p  and  cold  temperature  T.  Second,  it  is  assumed  that 
the  particles  are  ignited  at  the  point  of  closest  approach  to  a 
preceding  adjacent  particle.  (The  binder  burns  through  at  the 
thinest  point,  and  ignition  of  the  oxidizer  does  not  occur  until 
the  binder  burnthrough  has  occurred.)  Third,  the  con¬ 
sumption  of  the  oxidizer  is  by  a  spherical  outgoing 
deflagration  wave  at  rate  r0  which  emanates  from  the  ignition 
point  or  points.  Fourth,  after  the  deflagration  wave  has 
reached  the  point  of  closest  approach  to  an  adjacent  oxidizer 
particle,  there  is  a  binder  burnthrough  time  calculated  by 


tb=Axe"'  ■  (5) 

This  law  will  be  discussed  later.  Finally,  it  is  assumed  that 
oxidizer  ignition  is  instantaneous  after  binder  burnthrough. 

Researchers  and  practitioners  in  this  field  will  find  this 
model  unacceptable  from  the  standpoint  of  reality,  ft  is  not 
presented,  however,  for  its  realism  in  deflagration  physics.  It 
has  the  necessary  elements  to  test  the  statistics  of  the  oxidizer 
matrix  and  to  yield  information  of  the  effect  of  the  statistics 
of  the  net  propagation  rate. 

Equation  (5)  has  the  property  that  i*<* x  for  low  a,  but,  for 
general  a.  th  is  a  nonlinear  function  of  x.  For  large  *,  the 
binder  burnthrough  time  increases  at  an  increasingly  fast  pace 
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with  the  burnthrough  distance.  This  behavior  is  expected  to 
have  its  counterpart  in  reality.  The  self-deflagration  wave  in 
the  oxidizer  contains  an  advancing  thermal  wave  in  the  solid 
phase  which  will  precede  the  gas-phase  reactions  in  arrival  at 
the  binder  burnthrough  point.  This  mechanism  alone  will 
provide  energy  for  binder  gasification.  Following  the  gas- 
phase  arrival,  there  is  a  convective  flow  of  hot  gases  from  the 
receding  oxidizer  which  can  provide  more  heat  transfer  to  the 
binder.  Moreover,  a  diffusion  flame  between  the  binder  gases 
and  oxidizer  gases  may  form  and  will  add  further  heat  trans¬ 
fer  for  gasification.  However,  as  the  oxidizer  surface  recedes 
from  the  binder  burnthrough  point,  there  is  a  longer  gas  flow 
path,  and  a  thicker  heat-transfer  "boundary  layer”  will  exist 
at  the  burnthrough  point.  This  will  reduce  the  heat-transfer 
rate,  making  it  more  difficult  to  burn  through  the  binder. 
Consequently,  Eq.  (S)  embodies  this  expectation,  although  the 
exact  analytical  form  probably  is  not  correct. 

The  binder  could,  in  principle,  burn  through  to  a  particle 
adjacent  and  to  the  side  or  below  a  given  oxidizer  particle. 
However,  the  only  heat  source  to  a  binder  directly  below  the 
particle  of  interest  is  by  an  advancing  thermal  wave.  It  was 
shown  in  Ref.  4  that,  although  the  energy  in  this  thermal  wave 
is  often  sufficient  to  pyrolyzc  the  binder,  the  rate  is  often 
much  too  slow.  This  is  so  because  the  pyrolysis  rate  drops 
rapidly  because  the  surface  temperature  drops  rapidly  as  the 
initial  pyrolysis  carries  away  the  initial  portion  of  the  binder. 
Binder  burnthrough  points  to  the  side  of  a  particle  have 
additional  convective,  hot,  and,  perhaps,  reacting  gases  to  aid 
the  burnthrough  process.  Consequently,  it  will  be  assumed  in 
the  analysis  that  burnthrough  is  only  allowed  sideways  and 
not  directly  downward.  This  necessitates  consideration  or  at 
least  a  two-particle  column  width  for  the  analysis. 


Case  I:  Two-Column  Model 

It  is  first  necessary  to  demonstrate  that  the  ignition  details 
will  not  alter  the  results  to  be  presented.  Consider  Fig.  4, 
which  is  a  two-column  case  at  the  maximum  packing  density 
In  this  case,  burnthrough  at  the  bottom  of  particles  can  also 
occur.  Equal  time  contours  have  been  drawn  on  Fig.  4 
assuming  ignition  in  two  different  situations  at  lime  =  -  3 
units.  On  the  left  is  ignition  on  the  tops  of  two  particles,  and 
on  the  right  ignition  is  assumed  on  a  plane  cut  through  the 
propellant.  By  graphical  construction,  it  is  seen  that  the 
ignition  transient  disappears  after  only  two  particles  in  the 
columns  have  been  consumed.  After  that,  a  surface  profile 
that  is  periodic  in  time  emerges.  This  indicates  that  the 
ignition  will  have  no  lasting  effect  or  induce  any  wild 
oscillatory  behavior  in  the  sequence,  and  it  is  sufficient  to 
consider  an  arbitrary  set  of  particles  in  the  interior  of  the 
propellant  for  analysis. 

The  two-column  model  is  as  shown  in  Fig.  3a  but  with  the 
dotted  line  on  the  right  moved  one  column  width  to  the  left,  to 
exclude  the  third  column.  The  propagation  path  for  a  typcial 
element  is  ignition  at  point  o,  self-deflagration  to  b.  binder 
burnthrough  to  point  c,  and  then  the  process  repeats.  Con¬ 
sidering  a  long  column  of  length  D containing  N  particles,  the 
total  time  to  burn  through  distance  D  is 

•  =  1  AT,  +  +  'AP,  +'/>„  +  fAPj  +  ■  I  by  .  i  y 

Since  the  individual  binder  burnthrough  times  depend  upon 
the  i„ 


M$/» 

and  the  S,  are  random  variables.  The  expected  value  of  r  is. 
using  Eq.  (3), 
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(f*P  i  +  ,b/2) 


66,46  2 


(6) 


For  large  N,  D  =  N  rf/2.  Moreover,  the  burn  rate  is 
r=D/E(t),  so  that  r=d/2t,  where,  from  Eq.  (6), 


JX  JX 

The  distance  between  a  and  b  on  Fig.  3a  is  also  slightly 
affected  by  the  6  values,  but,  under  the  previous  ap¬ 
proximation  of  close  packing,  it  is  only  weakly  dependent 
upon  the  6  values.  Approximately,  rAP/  =  r/r0,  and  /  becomes 

-  r  1  t*m  r  *m  r 

r  +  /»„(jr)dM*,-r+ffc  (7) 

r0  J  J  "  '0 


The  expression  for  rAP,  is  where  the  assumption  of  the  dose- 
packed  hexagonal  array  enters.  The  final  relation  needed  is 
/„„(*).  Using  Eq.  (5)  for  the  burnthrough  law,  a  rather 
important  observation  may  be  made.  If  ihl} oex,  then  t 
becomes  the  time  taken  for  the  AP  to  burn  plus  the  time 
required  for  cookthrough  to  take  place  through  the  average 
binder  thickness.  That  is,  the  statistics  add  nothing  to  the 
problem.  Only  when  r„l:  is  a  nonlinear  function  of  x  does  the 
statistical  treatment  give  interesting  results.  The  reason  for 
this  behavior,  of  course,  is  that  if  ib/}<xx,  the  burnthrough 
rate  is  a  constant  and  the  slowness  of  cookthrough  of  thick 
layers  is  canceled  by  cookthrough  of  an  equal  number  of  thin 
layers.  When  the  function  is  nonlinear,  however,  the  statistics 
will  weight  the  slower  burnthrough  layers  in  a  heavier 
manner.  Placing  Eq.  (5)  in  Eq.  (7)  and  carrying  out  the  in¬ 
tegration  and  forming  the  overall  burn  rate  expression,  there 
results 


(2r/d)  +  {ih/d)Rr0 

(8a) 

rb=KAe-i 

(8b) 

R  =  2fcosba6„  -  / )  /  ( a6m ) 3 

(8c) 

In  Eqs.  (8),  ih  is  the  time  to  burnthrough  of  the  average  binder 
thickness.  #?a  I,  depending  upon  the  product  of  ai„;  this  is 
the  nonlinearity  effect  of  the  binder  burnthrough  law.  If  l'b  is 
short  enough  (if  the  binder  on  average  burns  at  a  faster  linear 
rate  than  the  oxidizer),  r/r0  may  be  slightly  greater  than 
unity,  according  to  this  model,  since  2r/dsl,  depending 
solely  on  the  oxidizer  loading  level. 

Comparing  with  experiment,  for  example,  with  the 
polysulfide-AP  work  of  Bastress,’  it  is  known  that  the  burn 
rate  ratio  of  Eqs.  (8)  will  increase  with  a  decrease  in  AP 
particle  size,  decrease  in  pressure  (the  propellant  has  a  lower 
exponent  than  does  pure  AP),  and  an  increase  in  oxidizer 
loading.  In  Eqs.  (8),  the  first  term  in  the  denominator  of  the 
burn  rale  expression  depends  solely  on  the  oxidizer  loading 
level,  whereas  the  second  term  depends  upon  all  parameters. 
The  factor  tb  Rid  will  decrease  with  AP  particle  size  at  a  rale 
depending  on  the  magnitude  of  aA  and  a6m;  consequently, 
the  particle-size  effect  is  qualitatively  predicted.  If  the  binder 
burnthrough  law  is  pressure-independent,  as  it  is  expected  to 
nearly  be,  the  pressure  behavior  is  also  qualitatively  predicted 
properly,  since  r0  increases  with  pressure  in  the  second  term 
of  the  denominator.  Furthermore,  the  oxidizer  loading  effect 
is  also  properly  predicted  because  Rid  will  decrease  with  an 
increase  in  oxidizer  loading  much  faster  than  Ir/d  will 
decrease. 

It  should  be  mentioned  that  there  is  a  serious  limitation 
when  compared  with  experiment.  In  the  experiments  cited, 
t/rt  was  significantly  greater  than  unity,  except  at  suf¬ 


ficiently  large  particle  size.  The  current  model  cannot  produce 
such  results  unless  the  pure  AP  rate  is  augmented,  perhaps  by 
consideration  of  an  alternate  heat  source  from  oxidi/er-bindcr 
reactions.  On  the  other  hand,  there  are  several  AP-bindcr 
systems  that  do  exhibit  burn  rates  less  than  that  of  pure  AP 
over  a  significant  pressure,  particle-size,  and  oxidi/er  loading 
range. 6 

The  major  point  is  that  some  of  the  trends  apparent  in 
experiment  are  contained  in  the  statistical  treatment.  The 
deflagration  model  is  so  naive  that  it  was  not  expected  to 
show  the  full  picture  with  regard  to  various  variables. 


Case  2:  Three-Column  Model 

Returning  to  Fig.  3,  consider  the  three-column  model.  The 
question  now  is  whether  or  not  there  is  significant  interaction 
of  the  third  column  with  the  left  two  columns.  In  Fig.  3a,  a 
perfectly  ordered  array  is  shown.  Evidently  the  propagation 
path  abedef  takes  exactly  the  same  lime  as  path  ijklmn.  In  this 
case,  the  same  burn  rate  is  obtained  by  considering  either 
path,  and  the  third  column  adds  nothing  to  the  problem. 
Consider  next  Fig.  3b,  which  is  a  somewhat  disordered  array. 
Here  x^  —xtm  and  xdf  -x)k.  Again  path  abedef  takes  exactly 
the  same  time  as  path  ijklmn,  and  there  is  no  interaction 
between  the  third  column  and  the  left-hand  two  columns.  This 
kind  of  event  will  be  called  one  of  "type  A.” 

Figure  3c,  on  the  other  hand,  shows  a  different  type  of 
event.  Here,  if  the  nonlinearity  in  the  binder  burnthrough  law 
is  strong  enough,  tjk  ►  r^,  or  idc.  Consequently,  path  abedef  or 
o  is  faster  than  path  ijkl.  The  bottom  particle  has  been 
consumed  solely  by  the  left-hand  path, and  the  right-hand 
path  is  slopped  effectively  at  point  I.  This  event  will  be  called 
one  of  “type  B,"  and  the  identical  event  with  the  left  and  right 
columns  interchanged  will  be  called  one  of  “type  C.” 

Neither  type  B  or  C  events,  if  followed  by  the  same  type  or 
one  of  type  A,  will  alter  a  burn  rate  calculation  as  compared 
with  a  two-column  calculation.  All  that  is  happening  is  that 
one  column  is  being  dragged  along  by  the  other.  Since  the 
burn  rate  is  defined  by  the  deepest  penetration  divided  by  the 
time,  the  left-hand  column  in  the  case  of  type  B  or  A  in¬ 
teraction  is  the  sole  determinant  of  the  burn  rate.  Even  a  type 
B  followed  by  type  A  followed  by  type  C  event  does  not  in¬ 
crease  the  burn  rate  statistically.  This  is  so  because,  im¬ 
mediately  after  the  type  A  events,  everything  starts  out 
"fresh”  with  the  center  column  particle. 

Consider  now,  however.  Fig.  3d,  where  two  opposite 
unusual  events  follow  each  other;  in  Fig.  3d,  this  is  shown  as 
type  B  followed  by  type  C.  Here,  everything  would  start  fresh 
at  points  /  and  o,  but  now  th  is  so  long  that  the  left-hand 
columns,  t!t  would  have  had  to  be  considered,  and  it  would 
have  decreased  the  expected  burn  rate.  Here,  however,  path 
opqrs  allows  that  long  time  interval  to  be  bypassed.  This  is  a 
clear  interaction  of  all  three  columns  which  augments  the 
burn  rate.  The  question  is,  how  frequent  are  such  events? 
First  of  all,  the  event  of  Fig.  3d  cannot  happen  unless  the 
binder  burnthrough  law  is  nonlinear.  Since  the  distance  from 
a  to  e  and  i  to  m  is  always  the  same,  the  total  time  to  burn¬ 
through  the  binders  would  be  the  same  if  the  law  were  linear; 
all  events  would  be  those  of  type  A.  Consequently,  the 
nonlinearity  is  essential  if  the  statistics  are  to  give  interesting 
effects.  Second,  the  answer  may  be  provided  analytically,  as 
in  the  two-column  model,  but  the  problem  becomes  so 
algebraically  complex  that  it  was  done  by  computer. 

A  three-column  array  was  set  up  on  the  computer  with  all 
4,’s  chosen  by  a  standard  FORTRAN  random  number 
generator.  Again  the  uniform  distribution  of  Eq.  (3)  should 
be  reproduced  if  enough  samples  were  taken.  A  baseline  case 
where  t\  -0.5  r/r0  was  chosen.  An  AP-binder  propellant  was 
assumed  at  75 wt  of  AP  with  a  binder  specific  gravity  of 
1.3.  This  yields  d/2r  ■  1 .01  $ .  If  the  mean  binder  burnthrough 
time  were  equal  to  the  binder  burnthrough  time  over  the  mean 
distance  t\ ,  this  would  yield  r/r»- 0.677.  The  nonlinearity  in 
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about  25%  at  ram  275,  where  the  burn  rate  is  about  half  the 
zero-a  value.  This  suggests  that  the  two-column  analytical 
result  may  be  adequate  for  real  propellants,  since  burn  rate 
depressions  from  the  monopropellant  oxidizer  rale  often  are 
not  severe.  That  is,  the  effective  a  is  probably  not  very  large. 
The  major  conclusion,  therefore,  is  that  a  two-column  model 
probably  is  not  too  bad;  there  is  little  influence  of  neighboring 
columns  on  the  two-column  situation.  Another  way  of  saying 
this  is  that  the  horizontal  correlation  length  scale  is  only  of  the 
order  of  a  couple  of  particle  widths.  From  a  practical  stand¬ 
point,  this  says  that,  if  the  transient  problem  of  combustion 
through  a  column  that  is  about  as  wide  as  two  of  the  largest 
oxidizer  particles  can  be  solved,  an  important  breakthrough  in 
burn  rate  modeling  will  have  been  achieved. 

Considering  more  than  three  columns  would  elevate  the 
result  in  Fig.  5  toward  the  t'h  =  th  line,  but  not  very  much  at 
moderate  a.  There  is  a  point  of  diminishing  returns  which  is 
reached  rapidly  as  the  column  width  is  expanded.  The 
probability  of  a  compound  helpful  event  from  extra  columns 
becomes  too  small. 

_  As  is  reduced  toward  zero,  all  curves  collapse  to  the 
=  /ft  line.  If  th  is  depressed  by  depressing  K  in  Eq.  (5),  the 
curves  all  tend  toward  r/ru  =  1 .0.  The  other  effects  of  oxidizer 
loading,  particle  size,  and  oxidizer  size  are  the  same  with  the 
three-column  model  as  with  the  two-column  case. 


Fig.  A  Illustration  of  number  of  “helpful"  events  as  a  function  of 
the  binder  burnihrouxh  nonlinearity  parameter. 


the  binder  hurnthrough  law  will,  in  general,  reduce  this  value 
through  Eqs.  (8),  but  the  three-column  model  will  mitigate 
this  drop. 

The  results  are  shown  in  Fig.  5.  Shown  is  the  analytical  two- 
column  result  from  Eqs.  (8).  This  is  a  lower  bound.  An  upper 
bound  is  the  curve  (constant)  that  would  result  if  fh  =  fh  were 
assumed.  The  three-column  result  lies  in  between  and  was 
obtained  with  a  120-particle  column  depth.  With  fewer 
particles,  some  severe  scatter  can  occur,  and  approximate 
bounds  on  the  results  obtained  with  only  a  40-particle  depth 
column  are  shown.  The  scatter  becomes  worse  at  high  a 
(strong  nonlinearity)  because  any  deviation  from  perfection  in 
the  random  number  generator  strongly  weights  the  long 
binder  burnthrough  events.  This  is  also  the  reason  why  the 
exact  (120-particle)  curve  lies  close  to  the  upper  bound  of  the 
40-particle  curve;  deviation  from  perfect  sampling  heavily 
weights  the  long-time  events. 

The  breakaway  of  the  three-column  result  from  the  two- 
column  result  begins  precisely  with  the  onset  of  events  shown 
in  Fig.  3d.  Figure  6  makes  a  count  of  the  number  of  “helpful” 
events  vs  a.  Even  at  high  a,  the  number  of  helpful  events  is 
relatively  small  compared  with  the  number  of  particles  in  the 
column,  but  the  effects  on  rate  are  striking  because  an  ex¬ 
tremely  long-time  burnthrough  has  been  avoided  by  the 
helpful  event. 

The  three-column  result  makes  a  100%  augmentation  in  the 
two-column  result  for  ram 400.  The  augmentation  is  shown  to 


Conclusions 

1)  Using  a  simplified  deflagration  model  and  a  nearly 
ordered  array  of  spherical  oxidizer  panicles  in  a  binder 
matrix,  a  simple  model  predicts  known  effects  of  oxidizer 
size,  oxidizer  loading,  and  pressure  level  through  the  statistics 
of  the  particle  packing. 

2)  A  nonlinearity  in  the  binder  burnthrough  time  vs 
distance  law  is  required  to  obtain  interesting  results  from  the 
statistical  treatment. 

3)  The  transverse  correlation  length  scale  for  burning 
through  the  matrix  is  relatively  small.  Hence,  combustion 
modeling  of  a  column  of  propellant  only  a  couple  of  particle 
diameters  wide  would  be  a  useful  undertaking.  This  is 
probably  the  most  important  conclusion  of  the  work. 
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